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ProiesChet® THE THIRD EDITION 


WE regard this book as an introductory text, so we did not 
feel any need for major modifications for this new edition. 
However, the broadened use of semiclassical techniques in 
quantum mechanics, in both atomic and nuclear physics, 
prompted us to add Appendix VII which introduces the 
reader to asymptotic, or ‘large-j’, approximations for many of 
the quantities discussed here. These relations provide some 
insights into the transition from quantal to classical 
descriptions of the angular momentum aspects of various 
phenomena, as well as often providing ways of simplifying 
computations that involve angular momenta that are large 
compared to Planck’s constant fh. Such large-j expressions 
can be quite accurate, even if the j-values are not very large, 
and their use can provide classical insights into the workings 
of complex phenomena. In other cases, large-j relations are 
valid only if some averaging over quantal oscillations is 
performed, emphasizing that the transition from quantal to 
classical descriptions need not converge uniformly. 

We also took the opportunity to correct the few misprints 
remaining in the second edition. 

A number of other, more detailed, books pertaining to this 
subject have appeared since our second edition was published 
in 1968. A nice treatment by R. N. Zare, Angular momentum: 
understanding spatial aspects of chemistry and physics (Wiley, 
New York, 1988) provides many examples of applications in 
atomic and molecular physics. Another text by I. Lindgren 
and J. Morrison, Atomic many-body theory (Springer, Berlin, 
1982) makes use of angular momentum theory in detailed 
applications of the theory of atomic structure. In particular, 
they use the graphical techniques discussed here in Chapter 
VII, but they find it physically more transparent to reverse 
our convention for the direction of the arrow on each leg of the 
graphs. We might also mention Angular momentum in 
quantum physics and The Racah—Wigner algebra in quantum 
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theory, both by L. C. Biedenharn and J. D. Louck (Addison- 
Wesley, Reading, Massachusetts, 1981), and Jrreducible 
tensor methods by B. L. Silver (Academic Press, New York, 
1976). Last, but not least, an English translation of the very 
extensive compendium of formulae and_ relationships 
Quantum theory of angular momentum by D. A. Varshalovich, 
A. N. Moskalev, and V. K. Khersonskii has been published 
(World Scientific, Singapore, 1988). Other useful results and 
applications are given in refs. [88-95]. 


May 1998 D.M.B. 
Italy and Oak Ridge G.R.S. 


PREFACE TO'THE SECOND EDITION 


SINCE publication of the first edition, a number of books have 
appeared that treat various aspects of the quantum theory of 
angular momentum. One, by Yutsis, Levinson, and Vanagas 
(Israel Programme for Scientific Translations, Jerusalem, 1962), 
develops the graphical methods introduced by Levinson. We 
have added a new chapter which gives a brief introduction to 
graphical techniques. We might also mention Quantum Theory 
of angular momentum, edited by C. L. Biedenharn and H. Van 
Dam (Academic Press, New York, 1965) which reprints a 
number of important papers in the field and includes an 
extensive bibliography of basic theory and applications up to 
1965. In particular, detailed applications to the shell theory of 
nuclear structure have been given by A. de Shalit and I. Talmi, 
Nuclear shell theory (Academic Press, New York, 1963), while 
B. R. Judd has described Operator techniques in atomic spectro- 
scopy (McGraw-Hill, New York, 1963). Additional tables of 
numerical values for the various coefficients have now appeared; 
we might mention the Tables of Racah coefficients by A. F. 
Nikiforov, V. B. Uvarov, and Yu. L. Levitan (Macmillan, New 
York, 1965), and A table of Clebsch-Gordan coefficients by B. E. 
Chi (Report prepared by Rensselaer Polytechnic Institute, 
Troy, New York, 1962). 

Errors and misprints found in the first edition have been 
‘corrected. There are only two changes of notation: The Fano 
X-coefficient has been replaced by the Wigner 9-7 symbol and 
the definition of certain quantities in section 6.1.2 has been 
changed to correct and clarify equation (6.21). The Wigner 6-7 
symbol has come to be used more frequently and formulae in- 
volving this function have been collected in Appendixes II and 
VI. A few other useful formulae have been added to the appen- 
dixes. 

We would like to thank many people who have brought 
various errors and misprints to our attention. We are especially 
indebted to Dr. K. T. R. Davies, Dr. H. J. Rose, Dr. O. Hausser, 
and Dr. J. Lopes for checking the text and many of the formulae. 
: D.M.B. 
March 1967 G.R.S. 


PREFACE TO THE FIRST EDITION 


DuRrine recent years important technical advances have been 
made in the quantum theory of angular momentum and its 
application to physical problems, both in nuclear and atomic 
physics. Our intention is to present these new techniques and 
to explain their physical significance without undue reference 
to their highly formal group-theoretic origins. Since this book 
was started other texts have been published by Rose [54], 
Edmonds [22] and Messiah [46]. Some overlap with their work 
is inevitable, but we feel our approach is sufficiently different to 
make a useful contribution. We have, in particular, attempted 
to emphasize the physical applications and to provide a source 
of formulae for workers in this field. Much of the underlying 
formal theory was developed by Wigner as early as 1937 [77], 
and is discussed in the book by Fano and Racah [31] and the 
well-known text of Wigner recently translated into English, 
[78]. 

References made to the literature on the physical applica- 
tions of this theory are necessarily selected somewhat arbi- 
trarily. We can hope to do no more than provide a starting 
point for wider reading. Where possible reference has been 
made to review articles rather than individual papers. We 
beg the indulgence of any authors whose work does not seem 
to be given proper recognition. 

Since there is a bewildering variety of notations and phase 
conventions it may be a help to the reader to have an indica- 
tion of some of those adopted in this book. Others can be 
found in the relevant sections. 

We use J and L to denote angular momentum vectors. 
Generally L refers to orbital angular momentum, but there is 
no fixed rule. In order to simplify formulae angular momentum 
is measured for the most part in natural units corresponding 
to h = 1. For spherical harmonics Y,, we adopt the usual 
phase convention of Condon and Shortley [17] and make 
frequent use of the ‘renormalized’ spherical harmonics 
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Crug = [(2k + 1)/47]*Y,, which help avoid unnecessary fac- 
tors of 47 in many formulae. There are many notations for 
Clebsch-Gordan coefficients (cf. Appendix I) but fortunately 
most definitions agree. We often use the more symmetric 
Wigner 3-7 symbol. There are two main definitions for 
rotation matrices and ours is explained in sections 1.4 and 
2.4 and in Appendix V. 

D.M.B. 
November 1961 G.R.S. 
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CHAPTER I 


SYMMETRY IN PHYSICAL LAWS 


1.1. Introduction 

A THEORETICAL investigation of a physical problem may come 
upon difficulties of two kinds. The exact physical laws govern- 
ing the behaviour of a system may not be known, making it 
impossible to arrive at a complete theoretical description of the 
system. Problems of fundamental particle structure and 
reactions present this first difficulty in an acute form. More 
often the situation is analogous to that encountered in prob- 
lems of atomic and molecular structure, where the inter- 
actions (Coulomb forces) are well-known, but where the struc- 
ture in a particular problem may be so complicated that no 
exact solution can be found. 

Fortunately the basic interactions in most physical problems 
have symmetry properties which affect the structure of a 
composite system in a way independent of the details of the 
interactions. Symmetries of physical laws may lead to con- 
servation laws, and the laws of conservation of energy, 
momentum, angular momentum and isotopic spin arise in this 
way. Again, in the theory of molecular structure the sym- 
metry of the configuration of nuclei in a molecule produces a 
symmetry in the electronic structure. An understanding of 
the effects of symmetry often enables one to distinguish 
between properties of a physical system which are conse- 
quences of conservation laws, and properties depending upon 
details of structure and interaction. For example, the angular 
distribution between two radiations emitted successively 
from a nucleus depends partly on symmetry properties, i.e. 
on the angular momenta of the states involved and partly on 
the detailed structure of the states. An understanding of the 
dependence on symmetry of the interactions enables one to 
obtain information about the detailed structure from experi- 
ments. A more familiar classical example is given by the 
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motion of a particle in a central field. The symmetry of the 
interaction leads to a plane orbit and conservation of angular 
momentum while the exact shape of the orbit depends upon 
the detailed form of the central interaction. 


1.2. The Symmetry of Physical Laws 

We develop the ideas suggested in the introduction by defining 
a symmetry transformation of a dynamical system as a trans- 
formation which applied to any dynamically possible state of 
motion carries it into another possible state of motion. As an 
example consider the solution to a classical problem involving 
two particles expressed by giving the positions of the particles 
as functions of the time. If the two particles are identical in 
their intrinsic properties and interactions then interchang- 
ing the trajectories of the particles yields a second solution 
to the problem. In the sense of the above definition the 
operation of interchanging two particles is a symmetry trans- 
formation. 

In quantum mechanics a symmetry transformation has a 
very simple representation. A state of motion of a quantum 
system is determined by giving the wave function y(t) as a 
function of time. This wave function represents a possible 
state of motion if the Schrédinger equation is satisfied; 


, ay 

Asymmetry transformation is represented by a linear operator 
U acting on the wave function with the property that Uy(t) 
satisfies the Schrdédinger equation whenever it is satisfied by 
p(t), i.e. if p(t) represents a possible state of motion then so 
does Uy(t). Symmetry transformations preserve the orthog- 
onality of wave functions and, in general, can be represented 
by unitary operatorst, so that U-1 = U+. If the symmetry 

} A unitary operator satisfies the condition UtU = UU+ = 1, where U+ 
is the adjoint of the operator, i.e. the adjoint operator U+ is also the inverse to 
U, Aunitary transformation preserves the orthogonality and normalization of 


wave-functions. Some symmetry transformations are represented by anti- 
unitary operators, cf. section 1.7, 
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transformation U is time independent then 


pel” _ any 
en eel 
ay 
hence th 7 alos 
and UH = AU. 


Thus the operator represeuting a time independent symmetry 
transformation commutes with the Hamiltonian of the system. 


1.3. The Symmetry Group of a Dynamical System 

An abstract group is characterised by the following prop- 
erties. It is a set of elements a, b... with a multiplication 
law defined satisfying the conditions that 


(1) itis associative a(bc) = (ab)ec, 

(2) there is a unit element 1 such that la = al =a, 

(3) every element a has an inverse a! with the property 
that aa-! = aa = 1, and a—! is itself an element of the 
group. 

The group multiplication law is often non-commutative so 
that ab + ba. A sub-group of a given group is a sub-set of the 
group elements which itself satisfies the conditions of a group, 
in particular that the product of any two elements of the sub- 
group must lie within the sub-group, and the inverse of an 
element of the sub-group must lie in the sub-group. 

The set of all non-singular square matrices of order n with 
the matrix multiplication law provides a particular example of 
a group. Conditions (1) to (3) are satisfied since matrix 
multiplication is associative, there is a unit matrix of order n, 
and every non-singular square matrix has an inverse. Also 
the product of two non-singular matrices is non-singular. 
This group is called the full linear group of order n. The uni- 
tary group of order 7 is the set of all m x m matrices A with 
the property that A~+ = A+ where A?* is the adjoint matrix to 
A. It forms a sub-group of the full linear group of order n. 
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Returning to symmetry transformations it is clear that the 
successive application of two symmetry transformations pro- 
duces a third symmetry transformation, and it can be shown 
that the set.of all symmetry transformations with this law 
of combination forms a group, the symmetry group of the 
system. For a quantum system the symmetry transformations 
are represented by unitary operators, and the successive 
application of two symmetry transformations S and 7’ is 
represented by the operator product 7'S. Froin section 1.2 we 
see the set of symmetry transformations of a quantum system 
is contained in the set of all unitary operators which commute 
with the Hamiltonian operator of the system. 


1.4. Geometrical Symmetries 

An important sub-group of symmetry transformations of 
many dynamical systems has a geometrical origin. The space 
of classical physics is described by a Euclidean geometry, 
implying that all points of space and all directions are equiva- 
lent and only statements relating to relative position and 
relative orientation have a meaning independent of the 
coordinate system. If the Euclidean character of space is 
reflected in physical laws no point in space nor any directions 
should be distinguishable absolutely by the performance of 
experiments. The operations of translation and rotation 
applied to a physical system should therefore belong to the 
symmetry group of the system. Thus the geometrical charac- 
ter of space determines symmetries of physical laws, or, per- 
haps more correctly, the physical symmetries determine the 
geometry of space. Lorentz transformations and space 
reflection come also into the category of symmetries with a 
geometrical origin. 

In practical problems often a distinction is made between a 
physical system and its surroundings and the effect of the 
surroundings is approximated by boundary conditions or by a 
set of known fields or forces applied from the outside. This 
approximate treatment may destroy the geometrical sym- 
metries partially or completely. For example, the equations of 
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motion of an atom in an external magnetic field are not 
invariant for arbitrary rotations unless the external field is 
also rotated. 

The geometrical symmetry operations of translation and 
rotation can be looked at in two ways, called active and passive. 
In the active sense introduced above they transform one state 
of a system into another. A rotation actually rotates the 
system from one position to another. The passive approach 


(b) ie 


Fig. 1. (a) Represents a rotation of the system. A vector OP is carried into 
OP’ by the rotation. 

(b) Represents the equivalent rotation of axes. The vector OP 
remains fixed, but its coordinates are changed. The position of OP 
relative to the new axes in (b) is the same as that of OP’ relative to 
the axes in (a). 


interprets the symmetry operations as coordinate transforma- 
tions, and the symmetry of the dynamical laws is expressed by 
stating that the equations of motion are invariant under 
translation or rotation of axes. Any translation or rotation of a 
system may, however, be induced by a coordinate transforma- 
tion and the two views are equivalent. This fact can be seen 
most easily by considering a simple example, the rotation of a 
one-particle quantum system through an angle « about the 
z-axis. (See Fig. 1.) Rotating the system transforms its wave 
function to a new one so that the value of the new wave 
function at the point (6+«) is the same as that of the old one 
at the point ¢ where ¢ is the azimuthal angle. Let the original 
wave function be y(r,0,¢) and the rotated wave function 
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y' = D,y then this argument shows that 
y'(r,6, 6+) = y(r,9,¢), 


or 
y' (r,9,0) = p(r,9, p—«). 


Thus the rotated wave function can be obtained from the 
original wave function ie making the coordinate transfor- 
mation (r,0,¢) — (7,0,¢—«), i.e. by a rotation of the coordinate 
axes through an angle —« Poe the z-axis. Both these views 
of the rotation are important and will be used with the con- 
vention that a positive (counter clockwise) rotation refers to a 
positive rotation of the system, hence a negative rotation of 
axes. 


1.5. Conservation Laws 

There is a relation between the geometrical symmetries of a 
physical system and the momentum and angular momentum 
conservation laws, having its origin in the fact that the linear 
and angular momenta are respectively the generators of 
translations and rotations of the system. As an example 
consider the canonical transformation generated by L,, the 
z-component of the angular momentum of a classical particle: 

If F is any function of the coordinates q; and the momenta 
p, of a dynamical system and « is small then the transforma- 
tion§ 

Pp, = pital{p.F}, 


fl 
q =a: ta{q,f}, oe 


is an infinitesimal canonical transformation in the sense that 
Hamilton’s equations of motion retain the same form when 
written in terms of p; and g;. F is called the generating func- 
tion of the transformation and it can be shown that if G(p,, q;) 

{ This convention agrees with Rose [54], Messiah [46]. Other authors use 
the opposite convention, a positive rotation referring to a positive rotation 


of the axes; Wigner [78], Edimonds [22]. Cf. section 2.4 and Appendix V. 
§ The Poisson bracket {G, F} of two functions of p,; and gq, is defined as 


aG OF 3G oF 
(|e 
=2 ag, Op, Ope O9r 
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is any function of the coordinates and momenta then the same 
function of the transformed coordinates is 


G(p;,.4%) = G(pi.d:) +a{G,F}. (1.2) 


The orbital angular momentum of a particle is the vector 
product L = rap and in particular L, = xp,—yp,. If we 
replace p, and q; in equation (1.1) by the rectangular coordi- 
nates and momenta of the particle and F by L, then simple 
calculation shows that the infinitesimal transformation 
induced is 


ae’ = 2—a4y, PD, = Pe Py; 
y =ytan, fy =o, (1.3) 
2! — a Pp. = Pp. 


It corresponds to a rotation of the particle position and 
momentum through an infinitesimal angle « about the z-axis. 
The angular momentum component L, generates infinitesimal 
rotations about the z-axis. 

In classical dynamics the rate of change of any function 
F of coordinates and momenta with time is given by the 
equation P 
— = {F,H} (1.4) 


with H(p,q) the Hamiltonian of the system. Invariance of the 
equations of motion under rotations implies that the Hamil- 
tonian should be unchanged by a rotation 


H(r’,p’) = A(r,p). (1.5) 
Transformation theory (equation (1.2)) requires, however, that 
H(r',p’) = H(,p)+a{H,L,} 
implying that debe =) 
aL, 
dt 


and from equation (1.4) 


Thus the motion is such that L, is constant in time, so we see 
that invariance under rotations about an axis (equation (1.5)) 
implies conservation of angular momentum about that axis. 


2 
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Similar considerations establish the connection between 
translation invariance of the equations of motion and con- 
servation of linear momentum. The laws of conservation of 
momentum and angular momentum are closely connected 
with the geometrical symmetries of translation and rotation, 
hence in a sense they are geometrical in origin. 

The above discussion has been given in terms of a simple 
classical system, but the results are general and hold also in 
quantum mechanics. The derivation follows an identical pat- 
tern in view of the close analogy between classical dynamics 
expressed in terms of Poisson brackets and quantum mecha- 
nics in the Heisenberg represcniation, reflected in the corre- 
spondence between the Poisson bracket and the quantum 
commutator 

VG —GF =(|F GG) Fai Gi. (1.6) 


F and G are functions of the coordinates and momenta in the 
classical case and of the corresponding operators in quantum 
mechanics. 

A transformation D applied to a quantum system can be 
interpreted according either to the Schrédinger or the Heisen- 
berg representation. If A is an observable then results of 
observations correspond to matrix elements of A,(1|A|2) 
between states of the system. The corresponding matrix ele- 
ment in the transformed system is (1|D+A D]2) and the trans- 
formation can be interpreted either as a transformation of the 
state vectors |1)-» D|1) and |2) — D|2) or as an operator 
transformation A — A’ = D+AD leaving the state vectors 
unchanged. The second interpretation is more appropriate to 
the Heisenberg representation and we shall use it to derive an 
expression for the operator D, which rotates a quantum system, 
from the classical results of equations (1.2) and (1.3) and the 
correspondence principle expressed in the Poisson bracket 
relation (1.6). When a system is rotated through an infinitesimal 
angle « about the z-axis an operator A transforms to A’ accord- 
ing to equation (1.2) 


A’ = A+a{A,L,} 
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to first order in «. The correspondence relation (1.6) implies 


ed 


A! = A——(AL,—L, A), 
yee 


Thus the operator D, for an infinitesimal rotation about the 
z-axis is ; 
Daa==n) 
a h aire 
By integrating the operator for infinitesimal rotations it can 
be shown that the operator for rotation through a finite angle 


a about the z-axis has the explicit form 


D, = exp(—ta«L,/h). (lez) 


ed 


In the following, to simplify formulae, we suppose that angular 
momentum is measured in natural units so that A = 1. 

Some quantum systems have non-classical internal degrees 
of freedom (‘spin’ degrees of freedom) in addition to the 
classical ones. The orbital angular momentum operator 
L = ) 7; 4 p; generates rotations of the classical variables, 
while the spin angular momentum operator S = },S, rotates 
the internal degrees of freedom. The generator for rotations of 
the system as a whole is the total angular momentum J = 
L-+S and J rather than L or S separately is conserved as a 
result of invariance under rotations. If L and S happen to be 
conserved separately, as is approximately the case in some 
atoms, it is a specifically physical rather than geometrical 


property. 


1.6. Commutation Rules for J 

Let D, = 1—10Jz be the operator which rotates a system 
through an infinitesimal angle « about the z-axis. This rotation 
applied to a vector operator A rotates it through an angle « 
about the z-axis. The components of the rotated vector A’ are 


Al =A,—«A,, A, =A,+oAy A, =A, (1.8) 
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From the discussion of section 1.5, however, Az transforms to 
A DE BL, Ul 
= (l+taJd,)A,(1—t20J,), (1.9) 
= A,+ta(J, A,—A, J,). 
Comparison of equations (1.8) and (1.9) yields 
eed =, ul hy, 
and similarly (J y-A.| eee 


ie ay Ad = 

Similar commutation relations are obtained with J, and J, by 
rotation about the x- and y-axes respectively. If A is the 
angular momentum vector J itself we get 


Id See 


J Jee le (1.10) 
Pdi lh 6s 
or JAI =i, 


for the commutation relations of the components of J. These 
commutation relations also follow directly for the orbital 
angular momentum L = r A p from the commutation relations 
ofr with p. 


1.7. Parity 

Reflection through the origin x ~ —2, y > —y, z > —z is 
a third symmetry operation of a geometrical and therefore 
‘intuitive’ nature. It differs from the operations of trans- 
lation and rotation in that it is discontinuous. Classically this 
implies that invariance under reflections leads to no conser- 
vation law, in the way that invariance under rotations leads 
to conservation of angular momentum. This is not the case in 
quantum mechanics. If P is the operator reflecting a system 
through the origin and P is a symmetry operation then by 
section 1.5 in the Heisenberg representation 


_éP 
ih =(P, H] = PH—HP = 0, 
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and the operator P is constant in time, or ‘parity’ is con- 
served. It has been found, however, that the interaction 
Hamiltonian responsible for 8-decay is not invariant under 
coordinate inversion and so does not commute with P. This 
may be true also of some other interactions so that parity 
would not be conserved for those interactions. 


1.8. Time Reversal 
Another symmetry operation often occurring in conjunction 


with rotational symmetry is the time reversal transformation 
[Wigner 78]. 


t— —t, 
Time reversal] transforms other dynamical variables as follows: 
ror, p—- —p;J— —J. 


A stationary state of a quantum system has a simple time 
dependence proportional to exp(—zHt/h), where H is the 
energy of the state. The time reversal transformation changes 
the sign of ¢ in this exponential factor but has no effect on the 
main part of the wave function. For this reason invariance 
under time reversal gives no conservation law and no addi- 
tional quantum numbers. On the other hand time reversal does 
say something about non-stationary processes leading, for 
example, to the law of detailed balance for nuclear reactions, 
[Blatt and Weisskopf 9]. 

As a further introduction to the ideas of time reversal con- 
sider a particle moving in one dimension. If the Hamiltonian 
of the particle is a function of the coordinate x and the momen- 
tum p then invariance under time reversal is expressed by the 
equation Hey i). 
Alternatively H is a real operator when expressed in terms of x 
and p = a If d(z,t) is a solution of the Schrédinger 


equation for the system then d(x, —#) is a solution of the time 
reversed Schrédinger equation, but since H is real the complex 
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conjugate $*(zx, t) is also a solution. If the phase of ¢ is chosen 
so that ¢ is real when t = 0 these two solutions have the same 
initial value at t = 0, and it follows that the solutions are 
identical for all times 


d(x, —t) = p*(z, f). (1.11) 
If ¢ is a stationary solution we can write 
d(at) = 6,(z)e 2". 


Equation (1.11) shows that the phases can be chosen so that 
$() is real, hence stationary states can be represented by real 
wave functions in the coordinate representation. It can be 
shown further that a complete set of wave functions (with 
definite phases) can always be found so that all matrixelements 
of operators invariant under time reversal are real. If we 
introduce a time reversal operator 6 for the above system by 


Od(z, t) = p* (a, t) 
then 6(¢,+¢,) = o3 +45 = 0¢,+64, (1.12) 
and (ao) =a" d= — a -0o, 


if ais a complex number. An operator which satisfies equations 
(1.12) is called anti-linear. A real wave function is invariant 
under this transformation. 

The results found in this special example can be generalized 
to apply to an arbitrary quantum system invariant under time 
reversal. A time reversal operator can be defined which is 
always anti-linear, but it cannot always be represented by 
simple complex conjugation. Suppose one can find a complete 
set of wave functions {|m)} invariant under the time reversal 
operation (thus with definite phases). If |x) is an arbitrary 
vector which is invariant under time reversal then 


B|a) = fa); 
but |x) and 6|«) can be expanded in the complete (invariant) 


set (72), J) = Y|m)(mla) (1.13) 
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and @|x) = > O(|m)(m|a)) 
= 5 (0|r2))(may* (antilinear 
° property of 7) 
= > |m)(mla)*. (|m) is invariant) 
i (1.14) 


Comparing equations (1.13) and (1.14) we see that if |«) is 
invariant the expansion coefficients (m|«) are real. A similar 
calculation shows that the matrix elements of any operator 
invariant under time reversal, between states of the set |m) are 
real. Thus only one real parameter instead of two is required 
to specify the value of a matrix element. Special problems 
arise in systems with angular momentum since the angular 
momentum operator is not invariant; but changes sign on 
time reversal. A further discussion of this point is given in 
section 4.9. 


CHAPTER II 


REPRESENTATIONS OF THE 
ROTATION GROUP 


2.1. Group Representations in Quantum Mechanics 

AN important part of the theory of groups is that concerned 
with the representation of their elements by matrices. If G 
is an abstract group then a representation of G with dimension 
n is a correspondence between the elements of G and a subset 
of the matrices of order n such that to each element a of the 
group G there is an » x » matrix F(a) with the property that 
if a and b are group elements then 


R(a)R(b) = R(ab), 


i.e. group multiplication corresponds to matrix multiplication. 
In a quantum mechanical formalism the elements of the 
symmetry group G of a system with Hamiltonian H are 
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represented by unitary operators in the Hilbert space of state 
vectors or wave functions, all of which commute with H 
(section 1.2). If H is an n-fold (n-finite) degenerate eigenvalue 
of H, there exists an n-dimensional manifold-@ in the Hilbert 
space such that all the state vectors in.W are eigenvectors of 
H with eigenvalue E. Let ¢ be any vector in-@ and S any 
symmetry transformation then the operator S commutes with 


ea HSd = SH¢ = S(E¢) = ES. 


Thus S¢ is also an eigenvector of H with eigenvalue H, and 
S¢ is also in the manifold .#: S transforms -# into itself. 
Choosing |1)...|n) as a set of orthonormal state vectors 
spanning @, then the transformation of W by the operator S is 
represented by the matrix S,, = (i|S|j). For if |a) is any state 
in-@ then |a) and S|a) can be expanded in the orthonormal set 
as |a) = > a,|t) and S|a) = >6,|7); but 
a 4 


Sla) = YaS|t) = ¥ a|9)G|S|*). 


Therefore > Sia, 


These matrices form a representation of the symmetry group 
G of the quantum system, for if S and 7 are in G then 


GISTIp) = E GSI ATs) 


and the matrices have the same multiplication law as the 
group elements they represent. Hence follows the important 
result that to every n fold degenerate eigenvalue of the 
Hamiltonian of a system there corresponds an n-dimensional 
representation of the symmetry group of the system. The 
representations of the symmetry group can be used to classify 
degenerate states of the Hamiltonian. 

There are an infinite number of ways of choosing the basis 
of the manifold.@ in the above discussion, and for each choice 
of basis the symmetry transformations are represented by 
different matrices. These representations are simply related 
and one can pass from one to another by a change of basis or 
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by a unitary transformation. Representations differing only 
by a unitary transformation are considered the same. 


2.1.1. Reduction of a Representation 

Consider a representation A of a group G on a manifold_4. 
Suppose there exist sub-manifolds.#, and.4, of W such that 
M+, =M and any matrix of the representation A trans- 
forms a state in.#, to a state in./, and a state in. /, to a state 
in.@,. If such a decomposition is possible then the represen- 
tation A is said to be reducible. Otherwise the manifold 4 is 
irreducible under the operations of the group and the represen- 
tation is irreducible. If a basis is chosen for the manifold.4 so 
that the vectors |1)... |r) span.@, and the vectors |r+1)... 
|v) span.@, then any matrix 7';; of the representation A takes 
the partially diagonalized form 


Wg | te te 


ao iat, = 900 LO 


i, eT 


and the sub-matrices of 7',; correspond to representations of 
the group G of dimension r and n—r respectively. The repre- 
sentation A has been reduced to a sum of two representations, 
A, of dimension r and A, of dimension n —r. This reduction is 
written symbolically as 

Ay= A,+A, 


Rn 


2.2. The Irreducible Representations of the Rotation 
Group 

In the following we confine our attention to systems with 

rotational symmetry and consider only the rotational sub- 

group of the symmetry group. Consider a finite manifold @ 
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associated with an irreducible representation A of the rotation 
group. Any rotation can be produced by a succession of 
infinitesimal rotations, thus a necessary and sufficient condi- 
tion for the irreducibility of A is that the manifold-W should be 
irreducible with respect to the Benetators J,, Jy, J, of infini- 
tesimal rotations. 

To find the irreducible representations we introduce the 
operators J, defined by 


J, =JI,tid,. (2.1) 


The operators J,, J,, J, obey the angular momentum com- 
mutation laws (1.10), and it follows that 


I See 


(2.2) 


Let |j) be the eigenvector of J, with the largest eigenvalue ). 
Then as 2.2 gives 


BE) = dep ysl) Sp ly) 


and J_|j) is an eigenvector of J, with ane j—1. Similarly 
(J_)*|7) is an eigenvector of J, with eigenvalue j —2 and so on. 
Let |j), |7—1), ... |g—r) with |j—r) = (J_)"|7) be a sequence of 
eigenvectors of J, generated by successive application of J_. 

Again J,J,|j) =(j+1)J,|j); but since j is already the 
largest eigenvalue of J, in the manifold.@, J,|j) must vanish. 
The square of the total angular momentum has the expression 


P= Ji+Ji +37 
=), eee a 
= J_J,4374J,. (2.3) 
It follows that 
Jj) = (J_J, 43245) |9) =9G+D/N), 


so that |j) is an eigenvector of J? with eigenvalue j(j+1). J? 
commutes with J_, hence |j—r) is also an eigenvector of J? 
with the same eigenvalue j(j +1). 

Because @ is a finite manifold the sequence |j) ...|j—r) of 
orthogonal eigenvectors of J, generated from |j) must 
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terminate, say, with r =n, i.e. J_[j—n) = 0. Now from 
equation (2.3) 
J[j—n) = (J, J_+0?7—J,)]j —n), 
= (F,—JJJ—n), 
= (9 eo 
but it was shown that J?]j—n) = j(j+1)|j—n), 


hence DG SIL) = (pa) eh 
ee) 
or Di 


The number 7 is a positive integer, hence 7 must be a positive 
integer or a positive integer plus one-half. The operators 
J, J_, J, transform the vectors |7) ...|7—) amongst them- 
selves, and since the manifold .@ is irreducible these vectors 
must span.@, so the representation is of dimension n+1 = 
27+1. The basis vectors |j) ...|7—n) are eigenvectors of J, 
with eigenvalues ranging in integer steps from +7 to —j. 
After normalization we denote these basis vectors by |j, m) 
where J,|j,m) = m|j,m) and —j < m < j. They are all eigen- 
vectors of J? with eigenvalues j(j7-+1) and the non-vanishing 
matrix elements of J are given by 


(jm +1 J, | jm) = {(Gam+1)(jFm)}H. (2.4) 


The phases of the off diagonal matrix elements are arbitrary, 
but normally they are chosen as above [Condon and Shortley 
17], thus fixing the relative phases of |jm) and |jm’). 

The matrix elements of the infinitesimal rotation operators 
are determined by the dimension (2j +1) of the representation 
once the z-axis has been chosen; thus the representation of 
dimension (2j+1) is unique. All other representations of the 
same dimension can differ only by a unitary transformation. 
This representation of dimension (27+ 1) is usually denoted by 
Q, and corresponds to an eigenvalue 7(j7+1) of J? with j 
integral or half-odd integral. In particular the basis |jm) of the 
representation Z, depends upon the particular choice of the 
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z-axis. The transformations of the basis corresponding to a 
change of axes will be discussed in section 2.6. 


2.3. Integral Representations and Spherical Harmonics 

The representations with 7 integral arise in problems con- 
cerning the orbital angular momentum of a single particle. 
The state vectors are functions of the particle coordinates, and 
the angular momentum has the explicit form 


0 
L=r1rap, withep, — —7— jcc. 
Ox 


or in polar coordinates, 


L, =L,1+iL,, 
alee (55 +7 cot 6 =) (2.5) 
06 0d)’ 
i) 7 
2 ad 
Also JL -|- ane ails sin =) -- z = (2.6) 
sin 6 06 06/ — sin? 6 dg? 


If we take eigenfunctions of L? and L, as basis vectors of the 
irreducible representations then these eigenfunctions are the 
spherical or surface harmonics, 


Yim(9$) = Om(I)Pn($); 


l l—m)"]} 
with @,,,(0) = (-1[5 a oI P™(6), if m> 0 


= oe i 0) 
©,,(g) = (2m) tems, 
P7"(9)(m > 0) is the associated Legendre polynomial [Jahnke 
and Emde 40]. This definition of the spherical harmonics 


involves an arbitrary choice of phase and we follow Condon 
and Shortley [17]. With this choice 


(Fim(9))* = (—1)™ ¥,_,n(99). (2.7) 
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An alternative phase differing by <' from the above is used by 
some authors [7] to give the spherical harmonics a more 
convenient transformation under time reversal. 


When m = 0 
Oa) 


+1\ 
Yi9(8, Q) = C-) P,(cos 0), (2.8) 


where P; (cos @) is a Legendre polynomial [Jahnke and Emde 
40}. 

Reflection through the origin replaces (6, ¢) by (7—6, 7+4) 
hence cos 6 —- —cos 6 and it follows from the property of 
Legendre polynomials 


PM(n—6) = (—/-™ PP), emer = (—ymeint 


that the spherical harmonics have a definite parity (—1)! 
where / is the order of the spherical harmonic. Spherical 
harmonics are normalized and orthogonal over the unit sphere 


[Xrm(9, $)*¥ in(8, $) EQ = 4(Wl’) 6mm’), 


where dQ = sin 6 d@d¢ is the infinitesimal element of solid 
angle, and 6(a db) is unity if a = 6 and is zero if a # b. They 
form a complete set for expanding bounded functions of 6 
and ¢. 

In some problems it is more convenient to use modified 
spherical harmonics with a different normalization 


47 \} 
Orml8, 8) =(5575) Vins) (2.9) 
With this normalization 
Ce o)y— ee, (cos u), (2.10) 
and Opa) = oy. (2.11) 


2.4. Explicit Representation of the Rotation Matrices 
In the previous sections we have found the possible irre- 
ducible representations of the rotation group from the commu- 
tation properties of the angular momentum operators J,, J,, 
J,. It remains to discuss the matrices representing finite rota- 
tions and for this purpose it is necessary to introduce a set of 
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parameters to specify arbitrary rotations. These parameters are 
the Euler angles («, 8, y). If a set of orthogonal axes (z, y, z) 
are rotated to new positions 2’, y’, and z’, the Euler angles are 
defined as follows (Fig. 2) by making the rotation in three 


Fie. 2. Illustration of the rotation defined by Euler angles (a,8,y). 


steps. First transform the axes to new positions (z,, y,, z) by 
rotating through an angle « about the z-axis. Next rotate 
through an angle 6 about the y,-axis to positions (xj, y,, z’) 
and finally make a rotation through an angle y about the z’- 
axis. Positive rotations are defined by the right hand screw 
sense. 
From equation (1.7) the rotation operator corresponding to 

this rotation has the explicit form 

D(a, B, y) = exp(—tyJ,,)exp(—iBJ, Jexp(—taJ,). (2.12) 
A little thought or alternatively, a direct reduction of equation 
(2.12) expressing Jy, and Jz’ in terms of Jy and Jz shows that 
the same rotation (a, f, y) will be produced by making first a 
rotation through an angle y about the original z-axis, then a 
rotation through an angle # about the original y-axis and 
eventually a rotation through an angle « about the original 
z-axis. Again from equation (1.7) the rotation operator D(a, B, y) 
is 

D(a, B, y) = exp(—iad,Jexp(—éBJ,)exp(—iyJ,). (2.13) 


| 
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Often we shall use R = (a, f, y) as an abreviation for the 
Euler angles specifying a rotation. If R, and R, refer to two 
successive rotations R,R, is the combined rotation. The 
inverse rotation to R is R-! = (—y, —8,—«). With this nota- 
tion we have 

D(R,)D( Ry) = D(k,f,), 


(2.14) 
D-\(R) = D(R-). 


In an irreducible representation of the rotation group of 
dimension 2/+1 corresponding to an angular momentum J 
the rotation (« 8 y) is represented by the matrixt 


IM"| D(x. B y)IM) = Dry ye(o By). (2.15) 
The operator D+ is the adjoint of D, hence its matrix elements 
are related to the matrix elements of D by 
(IM|D+|IM’) = IM'|D\IM)* 
— (ee. 
Also because D is a unitary operator: 


D*(a By) = D(a By) = D(—y —B —a), 
hence 
(Digs (a By ))* = Diy s¢( 1 A, 
The operator D is unitary 
DD — DD. = |. 
hence the rotation matrices 2.15 are unitary matrices, and 


> (Dir x(B))* Pieu( Bk) = 6(M, N), 


7 


> Dy R) Dyy(B))* = 6M, N). (2.16) 


t We use the convention of Rose [54], Messiah [46], that D(a B y) rotates the 
system through Euler angles (a 8 y), while others, (Wigner [78], Fano and Racah 
[31], Edmonds [22] and Rose [53}) use the opposite convention, i.e. that D(« B y) 
rotates the system through angles (—a —f —y). The relations between our 
rotation matrices and theirs are thus 


D(a By) corresponds to D(—« —f —y) 
Dy yx By) corresponds to 94,,(—a —B —y) = (—1)¥-¥ (QL (a B y))* 
diy(B) corresponds to d4,y(—B) = (—1)¥~* d4yy(8). 
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Because the basis vectors of the representation are chosen as 
eigenfunctions of J, and D has the form of equation (2.13) the 
matrices simplify as follows: 


(a By) = UMlexp(—taJ,)exp(—iBJ,)exp(—tyJ,)|IN), 
= e "eM +N) (T Mlexp(—ipJ,)|IN), 
= eNO diev(B)- (2.17) 


Phases of the rotation matrices depend upon the convention 
adopted for the Euler angles and on the choice of phases of the 
matrix elements of J. With the Condon and Shortley choice 
of phases (equation (2.4)) the reduced rotation matrices d4,y are 
real and can be expressed explicitly as 


j ARG Saag seo —n)!}* 
Se reser a 


% (cos B/2)\o a eeimangin B/2\" 9" 


where the sum is taken over all values a ¢ which lead to non- 
negative factorials. (Formulae for 7 = 4, 1, $, 2; Table 1.) 
In particular 
diuw(7) = (—1)'*™6(M, —N), 


de yy(2m) = (—1)%5(M, N). Ce 


Symmetry relations for the matrices dj,y are listed in 
Appendix V. The rotation matrices reduce to spherical har- 
monics when M or N = 0 


(x By) = (Cry (6 «))*. (2.19) 

The second of the os (2.18) gives 
dere (27) = (—1)" digye (0), 

thus if J is half-integral the rotation matrix for 8 = 27 has the 
opposite sign to the matrix for 6 = 0 and the rotation matrix 
is pericdic in 8 with period 47. Alternatively in the range 
(0, 277) the rotation matrices are double valued, the two values 
differing in phase by a factor of (—1). In terms of wave func- 
tions, a wave function corresponding to half-integral angular 
momentum changes sign on rotation through an angle 27 
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about any axis. This arbitrariness in the sign of the wave 
functions leads to no arbitrariness of observable quantities 
provided the initial choice of phase is used consistently 
throughout the calculation. Care should be taken to use the 
same phase for identical rotations. Wigner [78] gives a detailed 
discussion of the double valuedness of the half-integral 
representations. 

Besides unitarity the rotation matrices obey another 
orthogonality condition arising from a theorem of products of 
representations under group integration [Weyl [75]]. The 
theorem states that the products of matrix elements belonging 
to inequivalent representations of a group, and products of 
different elements of the same representation vanish on sum- 
mation over all group elements (integration for a continuous 
group). Applied to the rotation group the theorem gives 


Qr2Qr 7 
| ff Ohacl« By))* Diya By) sin B dB da dy 
000 
82? 
= Fag Dee AME, Bes ay 


The normalization factor 87?/27 +1 arises from the unitarity 
equation (2.16) of the matrices of the representation. 

The particular case of J = 4 deserves special mention 
because of its importance in discussion of spin. The com- 
ponents of J in the J = } representation are represented by 
2 x 2 matrices conveniently expressed in terms of the set of 
Pauli spin matrices o by J = }e. 

If J, is chosen to be diagonal and the choice of phases is 
made as in equation (2.4), then 


you 0 —-i i 
een | 0) em a) ee Oe 1)” 


The matrices o have the anticommutation properties 
9;0;-+0j0; = 264, 

and in particular of = a, = o? = 1. Together with the 2 x 2 

unit matrix 1 they are sufficient for a complete description of a 

J = 3 system. 


3 
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ie, Be 


TABLE 1 
Formulae for d,,,(B) for j =%, 1, 2 and 2 


$9 


pee MN 
ayy = dry = cos( 5) 


Ao ies Weer 7 
CES adie, ot Sie sin(§ 
Ci cos*(f) 


B 


ad, =da44= sin’(5] 
itil 11 9 


i sl a __ gil 
dy, = diy = —do4 


= —di, = sin B//2 


Lo diy, — cos%( 5) 
4 + i 
Ch = oe = 4; 
= —d2 4-4 
= —3 cost(5) sin( 5) 
ee 
yy = ay = hg 


eee 
a, = a24-5 
= cos(5)(3 cost( 5) 2] 
9 2 
en oy 
df= —dliy 


= sin( 5) (3 sine(®) —2) 


ON Ss Gea cost(5) 
dit ae —di, = Ss 
= d?,_, = —} sin B(1+cos B) 
din = do = de 29 
= di, = Vf sin’f 
doy aa di = —@i 


= —d?,, = } sin B(cos B—1) 


ale 
dss = dune = sin‘(5) 


di, = d*,4 
= 4(2 cos B—1)(cos B+1) 
di = dey 
= 4(2 cos B+1)(1—cos f) 
dio = doy ae —diy 


I 


—d*,, = ~v/2 sin B cos B 


4(3 cos?B—1) 


2 
doo 


I 


References for additional tables for d/,,,. (8) are as follows: 


we. 
tou 


j 


= 2, 4, 6: Buckmaster, H. A. (1964) Can. J. Phys. 42, 386 
1, 3, 5: —— (1966) Can. J. Phys. 44, 2525. 
3: Yine-Nan Cuiu (1966) J. chem. Phys. 45, 2969. 
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The rotation matrices take a particularly simple form 
because J, = 40, and of = 1. Thus 


exp(—ifJ,) = exp(—ifo,/2) = 1 cos B/2—io, sin B/2. 
Substitution of the explicit form of the matrices o, and 1 
gives dt,.(B) = ie 6/2 —sin a 
eM sin B/2 cos B/2/]° 


2.5. Rotation Matrices as Symmetric Top 
Eigenfunctions 

The rotation matrices D},(«, B, y) are eigenfunctions of the 
total angular momentum of a rigid body whose orientation is 
specified by the Euler angles «, f, y. (These angles measure the 
orientation of the principal axes (x’, y’, z’) fixed in the body 
relative to a set of axes (x, y, z) fixed in space.) The rotation 
matrices are also eigenfunctions of L, and L,, with eigen- 
values M and N respectively. If the rigid body has an axis of 
symmetry and the z’-body fixed axis is oriented in the direction 
of this axis then L,, as well as L, are constants of the motion, 
and the rotation matrices are eigenfunctions of the Hamil- 
tonian of the rigid rotator. These facts follow simply from the 
rotational properties of the wave function. Suppose ¢(f) 
where R = (« B y) is a wave function of the rigid rotator. 
A rotation of ¢(R) by R, = (x f, y,;) produces a new wave 
function ¢'(R) = D(R,)¢(£) and the value of the rotated wave 
function ¢'(R) at the point F# is the same as that of the old 
wave function at the point R’ which is carried into R by the 
rotation Ff, 1.e. 

¢' (Rk) = D(R,)¢(R) = o(F’). (2.20) 

If ¢(R) = ¢7y(R) is an eigenfunction of L*, L, and H the 
Hamiltonian with eigenvalues J([+1), N and EH, then ¢'(R) 
must also be an eigenfunction of L? and H with the same eigen- 
values. If the eigenvalue H has only rotational degeneracy 
the state ¢’(R) can be expanded in the set ¢,,,(f). Thus 


$(R) = brnl(R’) = & $r(R)(IM|D(R,)|IN), 
= » oru(R) Depy( Pr). 
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In this expression the ket |JN) in the matrix element stands 
for the state ¢;,(R). The relation takes an interesting form if 
R, = R = («f y), because then R’ = (0 0 0) = (0). Using the 
unitary property of the rotation matrices (2.16) we obtain 


bral R) = & (Phy B))*b1y(0). (2.21) 


Thus the 27 +1 wave functions ¢,,,;(#) are determined for all 
values of # by their value at R = (0 0 0) and the rotational 
invariance of the Hamiltonian. Equation (2.21) gives the general 
form of the wave function of an asymmetric rigid rotator, 
[42], [80]. When the rotator has an axis of symmetry (chosen 
to be the 2’-axis) there is a further specialization. An arbitrary 
rotation yy about the symmetry axis z’ should leave the wave 
function invariant up to a phase and L,, is conserved. This is 
possible if only one of the $,,(0) is non-zero and the wave 
function is (24,,(R))* apart from a normalization factor. 
The quantum numbers M and N are eigenvalues of L, and L,, 
respectively. 

Relation (2.20) may also be used to obtain explicit expressions 
for the components of L as differential operators. For example 
if #, is an infinitesimal rotation through an angle e, about the 
z-axis then R’ = (a—e,, B, y) if R = (a, B, y) and D(R#,) = 
(1—te,L,). Thus (2.20) becomes 


(1 —te,L,)b(a, B, y) fad o(a—e,, Bs y), 


od 
wpa, B, 2) Se By 
d i = =) g 
an eee 


Similar calculations using infinitesimal rotations about the z, 
y, and 2’-axes give§ 


F oe {- , Ae ag a1 a | . a 
= —1e —co = ,= -t 
* a ae sin B dyl’ dy’ 
t Bohr and Mottelson [13] use wave functions for the rigid rotator which are 
the complex conjugate of ours. 
§ To obtain LD, and L, we need expressions for R’ = 
(x-+da, B+dB, y +dy) 


in terms of the infinitesimal angles of rotation about the z, y, and z-axes. 
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while the expressions for L,, L, give ° 


WI) a f a 1 (= =) 2cos fp 0 
=| ap 88 ap aap 3e8" aya) Tes Be 


The above discussion also holds for half integral angular 
momenta; but then the rotation matrices lose their meaning 
as eigenfunctions of a classical rigid rotator. They are, how- 
ever, the approximate eigenfunctions of a particle of half 
integral spin coupled to a rigid rotator and in this form occur 
as the collective eigenfunctions of a deformed nucleus with odd 
atomic weight [13] [80] and of a molecule where the component 
of the electron angular momentum along the molecular 
symmetry axis is not zero. 


2.6. The Vector Model and Classical Limits 

As discussed in section 2.2 a representation of the rotation 
group is unique only up to a choice of basis for the manifold 
determining the representation. The basis vectors are chosen 
as eigenfunctions of the square of the total angular momentum 
J? and its z-component J,. There are, however, an infinite 
number of equivalent ways of choosing J, corresponding 
to all possible directions of the z-axis. If for example, we 
are applying a perturbation which has axial symmetry 
about some direction it is most convenient to take J, 
referred to this direction. Then the perturbing operator 
commutes with J, and the perturbed states remain diagonal 
Wat la 

If a set of axes (z’, y’, 2’) is obtained by a rotation F from a 
set (x, y, 2) then the eigenstates |JN)’ of J, are given by 
rotating the corresponding eigenstates |JN) of J, with the 


These may be obtained from expressions for the components of angular velocity 
of a rigid body and are 


de = «, cot B cosa + e, cot sina — e, 
dB = €, Sin a—e, COS & 


cos & sin % 
~ € sin B “*sin B* 


dy = 
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axes, i.e. by transforming the state |JN’) with the operator 


aay \JNy! = D(R)|JN) 
= S| MyIa| D(R|IN) 


= [J M) Dyy(R). (2.22) 


Equation (2.22) mee the eigenstates of Jz’ in terms of the 
eigenstates of Jz. The states (JN| conjugate to those in (2.22), 


rotate as “IN| = (IN|D+ = > (Diy) *(JM|, (2.23) 
M 


and we say they transform contragrediently (cogredience is 
defined by (2.22)). A symmetry property of rotation matrices 
(Appendix V) 
Dyy(B)* = (14 D2 57_y(R), 

shows that the transformation (2.23) for (J NV| is the same as that 
for (—1)*|J —N); that is these two quantities behave in the 
same way under coordinate rotations. 

Spherical harmonics afford a particular example of equation 
(2.22) 

Cr(8; o') = » Zl As B, v)Cim(9, 9). (2.24) 


The angles (0, ¢) and (6’, y’) are the angular coordinates of a 
point in the old and new coordinate systems. If n = 0 we 
obtain an addition theorem for spherical harmonics (cf. also 
section 4.6) 

P, (cos 6’) a CAB oe) GG, (2.25) 


using equations (2.19) and (2.10). 

Equation (2.22) leads to a geometrical interpretation of the 
rotation matrices. If we have a state with J? = J(J+1) and 
J, = M the indeterminacy of J, and J, is represented on the 
vector model by a vector J (of length +/{J(J +1)}) precessing 
about O,. If we make a measurement of the projection of J on 
another axis O,, inclined at an angle f to O, (Fig. 3) the prob- 
ability for finding a value M’ is just | Dip-y(« B y)|?. On the 
vector model we should expect to find the spread of values 
My, < M' < Mz shown in Fig. 3 due to the precession of J 
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about O,, with a probability P(M") of finding a value M’ 
assuming that J precesses uniformly about the O,-axis, 


PM ~(5°(1 —cos?f) —(M?+M'—2MM' cos f)]~*. (2.26) 


The effect of quantum indeterminacy is to allow values of M’ 
outside the limits (M,, M,) given by the vector model, but 


\ M’ 
Classical 
Ce limits, M,, M, 


Fie. 3. Rotation of the quantization axis from O, to O,. The classical limit 
indicates the spread of values for M’ to be expected from a classical precession 
about O,. 


| Fip-y)|2 falls off exponentially with M’ in this region. Also 
Dey oscillates within the limits (/,, M,); but for large J, ie. 
in the classical limit, | D},-,,|? averaged over several values of 
M’ to remove the oscillation is approximated by P(M’) 
obtained from the vector model. The asymptotic behaviour of 
the rotation matrices for large J has been discussed in detail in 
the W.K.B approximation by Brussaard and Tolhoek [15]. 

The general rotation includes the two Euler angles « and y 
(section 2.4) for azimuthal rotation about the old and new 
z-axis respectively. These play the role of phase angles only, 
occurring as a factor exp —7(Ma+WNy). Since 


| DE, ole: B y)|? = a nie 
they do not affect. the probability interpretation just described. 
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2.7. Coupling of Two Angular Momenta 

Often we work with systems made up of two or more parts, 
each with angular momenta. These may be different particles, 
or perhaps the spin and orbital properties of one particle. For 
the present we consider a system in which the total angular 
momentum J is the sum of components J, and J,. If an 
interaction between the two parts is such as to leave the indivi- 
dual angular momenta and their z-components constants of 
the motion, a complete set of commuting operators would 
include H, Jj, J,,, J and J,,. The corresponding eigen- 
functions |«j,j,m\m,) may always be written in the simple 
product form 

|aj,jgmyme) = x [2i1771 )| v7 2M2), 


where «, 8, and y represent any other quantum numbers needed 
to specify the states. In the following they will not be written 
explicitly. We have the eigenvalue equations 


TiljJemme) = jijit+ 1)|1J my), 

Fiji mye) = M]j,Jmymp), 
for Jj and J,, and similarly for J? and J,,. We could, however, 
choose a set including H, J?, Ji, J? = (J,+J,)? and J, = 
J,,+J2,, which contains as many physical observables as 
before. The eigenfunctions |7,j,J 17) now satisfy 


Jjjo7M) =ST+ jj), 
J lj 27M) = Mj.joJM), 


while Jj and Jj have the same eigenvalues as before. In 
physical applications this is often a more useful set. For 
instance when an interaction between the two parts of the 
system is introduced as a perturbation J? and J, may be 
conserved, but not the individual z-components J,, and J,,. 

The unitary transformation connecting these two represen- 
tations 


\jid27M) => | rJ2™MyMe) Hrjemymo|jyjoJM), 


MyM 


a ie. oe ae 2.27 
[jij gmyme) = 2 lisJxFM) aja FM|jjgmyme), 
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defines the vector-addition coefficient (j,j,m,m.|j,j.JM) = 
ijetM | jijemim2) (sometimes called a Wigner or Clebsch- 
Gordan coefficient). Often for brevity we shall write the vector 
addition coefficient as (j,j2m,m,|JM) or even as (mym,|JM) 
when confusion will not result. For givenj, andj, the values of 
J are restricted by the ‘triangular condition’ [Dirac [20]] 


jitje > J > |j:—Jal, 
and J ranges from j, +j. down to |7, —j,| in integer steps. Classi- 
cally J is the sum of j, and j., so the magnitude of the vectors 
must be such that they can form three sides of a triangle. 
The triangle condition is symmetric in j,j, and J as suggested 
by the classical vector picture. Since J, = J,,+J,, the vector 
addition coefficient vanishes unless M = m,+mz,. The ortho- 
normality of the eigenfunctions |JJ/) and |j,j7,m,m.) leads to 
the orthogonality relations for the coefficients 
> SMU jijomyme) jijemyma|J'M') = oJ, J’) 6M, M’), 

and 

PD grsateiat aid MET BA 9,920 73) = 6(m,, M;) 6(ms, my), (2.28) 


which express the unitary nature of the transformation (2.27). 
Since each coefficient vanishes unless M = m,+m, the sum 
over M is purely formal in the second orthogonality relation 
and in fact the sum is only over J. 

From a dynamical point of view the transformations (2.27) 
describe the addition of angular momentum. There is, however 
a geometric or group theoretic interpretation. The wave 
functions |j,j,m,m,) for a two-component system in angular 
momentum states j,, 7. span a (2j,+1)(27.+1) manifold. On 
rotation of the coordinate system the wave functions for the 
two components transform separately according to represen- 
tations 9, and 2, of the rotation group and the composite 
states bese forni as “equation (2.22) 

idem)’ = F DA, (RK) DE, (R)jjemm,). (2.29) 

Utes eery 
The basis states of the composite system transform according 
to a (2j,+1)(2j,+1) dimensional representation of the 
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rotation group denoted by 9; xB,.. This representation is 
reducible and the unitary transformation of equation (2.27) 
reduces it to its irreducible components and the states |J1/) 
are basis vectors of the reduced representation. It can be 
shown that J runs from j,+j, down to |7,—j,|. We may write 
the reduction symbolically as 
jitse 

Ix Bj, = PZ (2.30) 
Writing the reduction explicitly yields relations between 
rotation matrices. If D is an arbitrary rotation we have 


OR = IMAM) 
> UM|mym,)(mym,| D|nyn.)(mn2|JN ) 


MyMohyNe 
= 2 (IM |mym) Bin, Dring MUNIN ) (2.31) 
MyMoN Ng 


and similarly for the inverse relation 


Gr Di = DY RE De TN ts (2.32) 


MN; ~" Mens 


2.7.1. The Vector Model 

In terms of the vector model the state |j,j,J 1) is represen- 
ted by the two vectors j, and j, precessing in phase about their 
resultant J (which in turn precesses about the z-axis) (Fig. 4(a)). 
The precession of j, and j, about J and its projection on the 
z-axis represents the indeterminacy in their individual z- 
components m, and m, although their sum M remains con- 
stant. The square of the vector addition coefficient 


l(jrjomm,|J M |? 


is the probability that a measurement in the state |JJ.) 
gives the particular values m, and m, for J,, and Jg,. 
Conversely in the state |j,j.mym,) the two vectors precess 
independently about the z-axis (Fig. 4(b)) and | (J1Z|jijemime )|? 
is the probability that at any instant their resultant will be 
J. As in the case of rotation matrices the vector model 
suggests an expression for the squares of the vector coupling 
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coefficients in the limit of large quantum numbers. A compari- 
son of Figs. 3 and 4(a) leads to expressions for the vector- 
z 


aS 


(a) (b) 

Fic. 4. Vector picture of two angular momenta coupled to give a resultant J. 
(a) In the (JM) representation j, and j, are coupled to give a resultant J which 
is precessing about O,. m, and m, are undetermined. 

(b) In the (j,7,J,7,) representation j, and j, precess independently about O, 
leading to uncertainty in J. 


coupling coefficients which should hold in the limit of large 
quantum numbers. 
If j, > J), and hence J > j,, then 


GJ 2m,m,| JM) * d’ss-j{B) 
andl 705 Oandy) Si) 
(Jj,.MO|JM) = P, (cos B) 


where cos 8 = M/J. These and other relations are derived by 
Brussard and Tolhoek [15]; see also [96] and Appendix VII. 


2.7.2. Explicit Formula for Vector Addition Coefficients 
Recurrence relations for the coefficients can be obtained 
from the operator identities J, = J,,+J,,. In matrix form 
these become 
JM+1[J,|JM) 
= Y JM 41|nynzq)(nyng| J, +e 4.|mym2)(mm,|JM). 
MMs 


Matrix multiplication from the Icft by (mjm3|JM"), use of the 
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orthogonality relation equation (2.27) and the matrix elements 
of J, (2.4), give us 
{(J4£M+1)(JFM)} (mymJJM +1) 
= (Fm t+ ii +m,)}* (m,# 1m,|JM) + 
+{(JoF me +1) (jo +m ,)}* (mym > JM). (2.33) 
These relations are sufficient to determine the vector addition 
coefficients. 

The left-hand side of equation (2.33) vanishes if we take the 
upper sign and put M = J. With the normalization condition 
(2.28) this enables us to determine the various (mimo|JJ), 
apart from an overall phase; this we fix (following Racah [48]) f 
by the convention that (j1J—ji|JJ) is always real and posi- 
tive. The lower sign in equation (2.33) gives us (mime2|JM—1) 
in terms of (mimj3|JM), so by a ‘ladder calculation’ 
starting with 7 = J we get all the coefficients, which, we see, 
must all be real. After some heavy algebra along these lines 
Racah obtained the general formula 


(abaB|cy) = 6(a+8, y)A(a b c)x 

x [(2e-+1)(a-+-a)! (a—a)! (b+)! (bf)! (e+y)! (c—y) 1x 
x ¥ (—1)’[(a—a—»)! ieee (6+6—v)!x 

Cat Sheep nha ea. 

(ath —c\"(aee 26)! (bere ats 

se) = eae 8M 


and y runs over all values which do not led to negative facto- 
rials. 

The formula can be simplified when «=f=y =0 
(Racah [4]); if 2g =a+b-+c, (ab00|c0) = 0 if 2g is odd, and 
(ab00|c0) = (—1)?**(2c +-1)* A(abe)g! [(g —a)! (g —b)! (g—c)!]* 

(2.35) 


where 


if 2g is even (cf. section 4.7.2 and Table 2) 
t With identical results, Condon arfd Shortley [17] adopt the convention that 
iJaljt+e ji: +32) = +1 and that the matrix elements 
AIIM, lijed —1 M) 
are real and positive for all J. 
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Obtaining numerical values from ‘the general formula is 
tedious, but by assigning a definite value to one argument, 
say b, equation (2.34) reduces to simple closed forms. Table 3 
contains formulae for the more symmetric Wigner 3-7 
symbol defined in equation (3.3) for b = 0, 4, 1, 3, 2. Symmetry 
relations and other formulae are given in Appendix I. 
Extensive numerical tables have been published [58] 
[66], and for (ab}—}|c0) by de Shalit [18]. Values for 
ee) are fiver i laple 2 


TABLE 2 

BS z : 

Values of veal An asterisk means the symbol 1s 
negative. 

abe | abe | ab ¢ 
@ sk yee | 1 56 6 6/143 | 3 3 6  100/3003 
0 2 2 ibe mo 22a olan los one 20/1001 
O73 35 lite me 2 24 38185 3°55 6 7/420" 
Me gy io 253 3 4/105 4 4 4 18/1001 
OSoe 1) tiles 2s meal * 4 4 6 20/1287" 
OGG lle 2 2 4° 20/693" ao pools 
is es Te meray ets 4 6 6 28/2431 
} 2 8 3/458 2°59 5 10/229 5 5 6 80/7293 
Ge eo 6 leyfcla® 6 6 6 400/46189* 
1 4 5 5/99* aot 2a 7* 


2.7.5. Exchange Symmetry of Tivo-particle States 

When eigenstates of total angular momentum of two 
identical particles are constructed according to (2.27), and the 
individual angular momenta are the same, j,; = jo, the sym- 
metry of the states under exchange of the two particles is 
determined by the symmetry of the vector addition coeffi- . 
cients = eee 

Gjmm'|J My) = (—1)7 7 jjm'm| J aD). 

If we denote the state occupied by the 7th particle as |),, the 
coupled state may be written 
gM) = > |Jm)[ jm’) (gjmm'|J A). 


Pon 


= $5 {ljm),|jm’).+( —1)7 jm’), | jm) oh jjmm'| IM). 


mn 
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TABLE 3 
: a bee : 
Algebraic formulae for ce with c = t, 
1, 2 and 2. 

u (=P y ‘) = (Sse [ ut+tatl b 
(3 —a—$ }/ (2a +2)(2a+1) 
a a i _)a-a Eereateer + 
(; —a—l | = 2a(a+1)(2a+1) 
a a 1\ . a 
(7 —o Jee [a(a+1)(2a-+1)}+ 
(: al 7 ea anree [ (atat1)(a+a+2) 
a -a-—-l 1 (2a+1)(2a+2)(24+43) 
a a+l ) ye [Se 
(; —a 1) (a+1)(2a +1)(2a +3) 
@ att :) aon [Me rete et Dea) i 
(g Sa) > 2a(2a + 1)(2a-+2)(2a-+ 3) 
a a+s $\ _ ae ata+l ae 
(2 = 4) me ae TaT Tee | (Cae 
(: a+3 ;) Smears (at+a+1)(a+ta+2)(ata+3)7]¢ 
a ages (2a + 1)(2a-+2)(2a +3)(2a-+4) 
a a+3 : = (—)a-« 3(a—a+1)(a+a+1)(a+a+2)]t 
ie eo (2a +1)(2a + 2)(2a +3)(2a +4) 
ft a 5} See [ee ee Mea) t 
a —a—-2 2) a(2a +3)(2a+2)(2a+1)(2a—1) 
a a 2 eee 3(a—a)(at+a+]) 
( oo) ‘| =) | soa oe ee | a 
(: a 0) See 3a?—a(a+1) 
a —o Of} [a(a +1)(2a +3)(2a +1)(2a —1)]t 
he a+l1 3] jane (a+a-+-1)(a+a+2)(a+a+3)(a—a) 
a -—a—2 2 a(a+1)(2a +4)(2a+3)(2a4+1) 
a atl 2\ _ anes (ata+2)(a+a+1) = 
Ge ni Fearea ieee Coie 
(: a+l 3) eat [ 3(a+a+1)(a—a+1) + 
a —o Oo} e a(a+1)(a4+2)(2a+3)(2a 41) 
c a+2 >) Beas (a+a+1)(a+a+2)(a+a+3)(a+a+4)]¢ 
a —a—2 2 (2a + 5)(2a + 4)(2a + 3)(2a4-2)(2a+1) 
(; a+2 1) = (=)eret [teceack Me tee Metatee—at) t 
“a —a—l 1 (a+1)(a@+2)(2a+1)(2a +3)(2a +5) 
(E a+2 Al aya SS? + 
a ao oO; (a+1)(2a + 5)(2a +4)(2a +3)(2a +1) 
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Thus when particles 1 and 2 are interchanged it merely multi- 
plies the state vector |jjJm) by (—1)’~”. That is, the state is 
symmetric or antisymmetric under exchange as (J —2)) is 
even or odd. If 7 represents an orbital angular momentum, 2) 
is even, and the condition is whether the resultant J is even or 
odd. On the other hand if j is the total (spin plus orbit) for 
each particle (j7—7 coupling), 7 will be integral or half-integral 
according to whether the particles are bosons or fermions. 
Thus both the symmetric boson states and the antisymmetric 
fermion states will have J even only, odd J states having the 
wrong symmetry in both cases. 

A very simple example is given by the total spin states of 
two spin —4 fermions, 7 = 4. The singlet J = 0 state is anti- 
symmetric, the triplet J = 1 state is symmetric. 

Of course, this simple property no longer holds when), 4 j., 
and the exchange symmetry is no longer determined by the 
vector coupling. 


CHAPTER III 


COUPLING ANGULAR MOMENTUM VECTORS 
AND TRANSFORMATION THEORY 


3.1. Transformation Theory 

OrrEeNn there are several degenerate but independent states 
|x), |x), ... which are eigenstates of some operator a with the 
same eigenvalue «. Thus we require other labels to distinguish 
them. These may be provided by finding another operator B 
which commutes with a. Hence 8 has eigenstates which are 
linear superpositions of the |x); belonging to the same eigen- 
value «. It is then straightforward to diagonalize in this 
subspace to find these superpositions |«,6) which are now 
labelled by « and the various f. If two or more combinations 
are still degenerate, i.e. share the same eigenvalues « and f, 
we need yet another operator y which commutes with both 
a and @. So we proceed until we have resolved the original set 
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of states |x),, degenerate in «, into a set completely non- 
degenerate in the eigenvalues «, f, y, ... of a set of commuting 
operators a,B,y, ... Such a set is called complete and represents 
the maximum number of simultaneously measured data 
allowed by the uncertainty principle. 

The most important of these operators is usually the Hamil- 
tonian H of the system or at least the principal part of it, 
the rest perhaps to be treated as a perturbation later. The 
other operators then have to be chosen to commute with H. 

The need for the present chapter arises because there are 
often two or more ways of choosing our complete set. These 
sets of course are not independent but have eigenstates which 
are related by unitary linear transformation. If two such sets 
are denoted by A and B we may write 


|A) = Ye|B)(BA). (3.1) 
The expansion coefficients (B|A) are the transformation 
amplitudes and form a unitary matrix. We have already met 
with two examples in sections 2.6 and 2.7. One is the trans- 
formation (2.22) linking states referred to differently oriented 
quantization axes, the other is the change of representation 
(2.27) for states comprising two angular momenta. Below we 
extend the latter example to systems made up of more than 
two angular momenta. 


3.2. Scalar Contraction of Angular Momentum States 

We have so far met the vector addition coefficient as giving a 
rotationally cogredient linear superposition of the products of 
two functions which separately behave cogrediently under 
rotations. Alternatively, it projects out the various irreducible 
parts of such products. A different point of view put forward by 
Wigner [78] considers the coupling of three angular momentum 
vectors to a zero resultant; and because the wave function of 
the coupled state has zero angular momentum it is independent 
of the choice of axes, i.e. it is a scalar or invariant quantity. 
Equation (2.27) and the knowledge that 


(coy —y|00) = (—1)?/4/(2e +1) 
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make it easy to perform, the scalar contraction product of 
three state vectors and we find that 


Dasy 12x) BB ley )(abaBle—y)(—1)°~” (3.2) 
is an invariant. 

We shall meet the contraction (3.2) in another guise in section 
4.7 where we discuss the Wigner—Eckart theorem. 

Among the various notations (listed in Appendix I) used for 
vector addition coefficients the Wigner 3 —j symbol emphasises 
this contraction property. It is related to our present trans- 
formation coefficient by 


b 
(° 5 )veety = (I? "abaple—y), (8.3) 
provided, of course, that «+$8+y = 0. Now the result 3.2 
becomes 


Dee laxS|O ley) (2 ; ‘) = scalar invariant. (3.4) 


In addition, use of the 3—j7 symbol often facilitates algebraic 
manipulation because of its high degree of symmetry. It is 
invariant under an even—cyclic—permutation in the order of 
its arguments, while an odd—non-cyclic—permutation merely 
multiplies it by (—1)*t?**. It thus avoids the unsymmetrical 
surds and phases appearing in the corresponding Clebsch- 
Gordan symmetry relations. The origin of these symmetries 
also appears more clearly through the contraction (3.4). All three 
angular momenta a, b, c, enter on an equal footing, thus the 
scalar invariant should be unchanged within a factor +1 by 
re-ordering them. A special case with a = 6 =c = 1 is the 
triple scalar product of vectorst 


UAVv°>W=VvVAW:U=WAU:V 


. 111 
= —11/6 31, wr .0r( : ih 


The 3—j symbol also has the visual advantage of displaying 


{ The quantities u,, v,, w, are the spherical components of the vectors u, V, 
and w defined in equation (4.10). 
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z-components on a different level from the vectors in accord 
with their inferior role. 

Specializing equation (3.4) to the case c = 0 (and therefore 
a = b, «+f = 0) allows us to define also a 1—7 symbol, the 
analogue of the metric tensor. We have 


Dap laaoo)(,",) = a scalar, 
where (,",) (ee ley (3.5) 


In many physical problems measured quantities are inde- 
pendent of the choice of coordinate axes, hence these quanti- 
ties may be expressed in terms of scalars. In the course of a 
calculation of such a quantity we may be confronted with a 
complicated expression containing products of vector addition 
coefficients, and the solution of the problems often lies in 
extracting the various scalars. In the following sections we 
give some of the scalar invariants which are common to many 
problems. 


3.3. Coupling of Three Angular Momenta 
When we have three angular momentum vectors we may use 
an uncoupled representation 


|J171, JgMz, JgMz), 
or one in which the vectors couple to a resultant J and M, that 
is, an eigenstate of J? = (j,+J,+J,)? and J, =j,,+32,+39:- 
However, the latter is no longer unique, and we require a 
further quantum number. There are three possibilities: We 
may couple j, and j, to form J,,, then add j, vectorially to 
give J. First 


iJ eFi2My2) = re 1912 )|J 2M) (frJ aM M9| Jy 2M»), 
then 


H(jJ2)Fi12J33 JM) = Dane i13 2F12-My2)|j3™2)(Sr2j3 My gms| JM). 
(3.6) 


This state is also an eigenfunction of Jj, = (j,+j,)?, and J;, 
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provides the additional quantum number which specifies the 
state of the coupled 3-vector system completely. 

Alternatively, we may first combine Jj, and j, to give J, and 
then add j, to give J, that is 

[J2js723Mo3) = rn, [Jeol Ja73)(JoJaM2Ms| S03 Mos) 
then 
[Jus (Jeds\Je3; IM) = 2 are, 117% 19 2534 23 Mes) 
Jiticg M2 Mold M). (3.7) 

So we now have an eigenfunction of J5, = (j,+j3)%. Similarly 
we could have used J,3, the resultant of j, and js, to define a set 
of states |(j,j3)J13, J2; JM). Clearly, these three representations 
are not independent, and since they span the same subspace 
they must be connected by a linear transformation (equation 
(3.1)). For example 


(jJadi2 Js; IM) = ie [Jus (J2Ja)4033 JM) x 
X (ji (Joja)F 233 J |(jJa)Fr2Ja7) (3.8) 
The transformation coefficient is a scalar and independent of 
M. We use the coefficient to define the Racah W-function [48], 
hr Gojs)F 235 J |(hJe)4i2J33 J) = 
= [(20,g-+1)(2Jag +1) VW (jj IsiSr24ea), (3-9) 
whose normalization is chosen to simplify its symmetry 
properties. We may easily express it in terms of vector 
addition coefficients. We use (3.6) and (3.7) to expand both 
sides of (3.8) in the uncoupled representation |j1m1j2m2)3m:3 ). 
Equating coefficients we obtain 
ij ota s| Jae Mie) (Sie jg Miz mslJM) = 
= ae (JoJgmgms| Jo M3) (j1I 93 m,M»,|J M) x 
X (rs (J2Ja)F 033 J|(jid2) Sie, Jsi J). (3.10) 
For clarity, we rewrite this using a, 6, c, d, e and f for the 
vectors and a, f, y, 6, e, and for the corresponding z-com- 
ponents 


(abaBlex +f) (eda+p, y—«—Blcy) 
= >, dB, y—«—B| fy —«)(afa, y—aley) x 
x[(2e +1)(2f+1)]* W(abcd; ef). (Jan) 
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Multiplication of both sides by (6d 8, y—a—Alf’y—a) and 
summing over f leads to the further relation 


(afa, y —alcy)[2e-+1)(2f-+1)}*W(abed; ef) = 


= ¥, (abaflea+B)(eda+f,y —«—Bley)(bdB, y —« —B] fy —a). 
(3.12) 
In a similar manner we finally obtain 


[(2e +1)(2f-+1)]! W (abcd; ef) 6 oly | = 
= a (abap|ex +f) (eda +8, ee x (323) 


x (bdB, y—a —B| fy —«) (afay —ale’y’) 
in the summation of which y is held constant. In this last form 
the W-function appears not as transformation coefficient but 
as the scalar invariant obtained by contraction of four vector 
addition coefficients. In this, perhaps, lies its greatest value. 
The symmetry properties of the W listed in Appendix II 
can be obtained easily from equation (3.13) by considering the 
corresponding symmetries of the vector addition coefficients 
involved. Equation (3.13) also makes it clear that the following 
triads of vectors have to satisfy ‘triangular conditions’ (section 
(acf), (abe), (bdf), and (cde). In diagrams introduced by 
Levinson [85] lines represent angular momenta and triangular 
conditions must be obeyed at each vertex (cf. Chapter VII). 


: (3.14) 


Equation (3.13) may be put into a form which shows the 
symmetry of the Racah W-function in a more obvious way by. 
using the 3—j symbol defined in equation (3.3) (ef. Appendix 
II). This form is probably the most useful for calculations. 

A general formula for the Racah functions may be found 
from (3.13) by using expression (2.34) for the vector addition 
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coefficients. Racah [48] succeeded in reducing this to 


pel = A(a be) A(acf) A(bd f) Ate de) x 
x >, (—1/(a+b+e+d+1—2)!z!(e+f—a—d +z)! x 
i oem, soe (a+b —e—z)!(c+-d—e—z)! x 
x (a+e—f—z)\(b+d —f—z)!]>, (3:15) 


where A(a b ¢) is given by equation (2.34). 

In the same way as the general expression for vector addi- 
tion coefficients (3.15) may be reduced to simple closed forms by 
assigning a definite value to one argument, say e. Formulae 
fore = 0, 3, 1, 3/2, 2 have been given by Biedenharn e? al. [7], 


TABLE 4 
Formulae for W(abed; ef) with e = 4, and 1 
W(aa+4bb +43 4c) 


= (=yrrre | (a+b+e+42)(a+6—c+1) i 
= (2a + 1)(2a + 2)(26 +1)(26 + 2) 


W(aa-+ }bb—-}; te) 
= (—)atb-e Reesvessonenle 
(2a+ 1)(2a+ 2)26(2b +1) 


W(aa+1bb +1; 1c) 
= (—)at-e [a ee ESS Sy 
5 4(2a+3)(a+ 1)(2a+1)(2b-+3)(28 FYO+Y 


W(aa+1bb; le) 
a ae Pee ea 
4(2a-+3)(a+1)(2a-+1)0(b-+1)(20-++1) 


W(aa+ 1bb—1; le) 
= Aare RE 
aioe [ 4(2a+3)(a + 1)(2a4-1)(26 —1)(26 +1) 


W (aabb; 1c) 


= (—)yatb-ent a(a+1)+b(b+1)—e(e+1) 


[4a(a+1)(2a+ 1)b(b + 1)(2b-+1)]$ 


those fore = 0, 4, 1 can be found in Table 4. Extensive numeri- 
cal tables have been prepared by Biedenharn [4] Simon et 
al. [66], and Rotenberg e¢ al. [58]. 

The unitary nature of the transformation (3.8) which defines 
the Racah function leads to the orthogonality relation 


Dd, (2e+1)(2f+1) W(abed; ef) W(abcd; eg) = 6,,. (3.16) 
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We can use the expansions (3.6) and (3.7) of the three-vector 
states to get relationships between the transformation coeffi- 
cients (3.9). For example, the symmetry properties of the 
vector addition coefficient (Appendix I) show that 


la, (db) fe) = (—1¥** a, (bd) fi), (3.17) 
so that 
(a, (db) f; cl(ab)e, d; c) = (—1)-?- a, (bd) f; e[(ab)e, d; ¢) 
SS es adh I 
x W (abcd; ef) (3.18) 


from equation (3.9). Similarly we can show that 
((bd)f, a; e|(ab)e, d; c) = (—1)*%-* (a, (bd)f; ¢|(ab)e, d; c) 
(3.19) 
A fruitful source of sum rules is the closure relation for 
transformation coefficients 
(A|B) = Yo (A|C)(O|B), (3.20) 


where the sum runs over a complete set of eigenstates |C). 
Applying this to the change-of-coupling coefficients we can 
write, for an 


(a, (bd) f; c|(ab) e, dsc) = 
= (a, (bd)f; c|(ad)g, b; c)((ad)g, b; e|(ab)e, d; ¢), 

which using equation (3.19) gives the Racah sum rule 

De (29 +1)(—1)?+9 W(aded; gf) W(abde; eg) = W(abcd; ef), 

(3.21) 

where p = a+b+c+d+e+f. Other sum rules derived in a 
similar way and some other relations are displayed in Appen- 
dix II [5], [7], [23). 

Parallel to his treatment of the vector addition coefficients 
as contraction symbols Wigner [77] has defined a 6 —j symbol. 
This differs from the Racah function in phase only, 

ie be 


Bin | = (—1)2t?t+e+4 W (abcd; ef). (ee) 
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The 6—j symbol has somewhat higher‘symmetry being invari- 
ant under interchange of any two columns and also under the 
interchange of the upper and lower arguments in each of any 
two columns (Appendix II). Some extensive tables of the 6—j 
symbol are now available [38], [58]. 


3.4. Four Angular Momentum Vectors 

A state in which four angular momenta a, b, d, and e have a 
resultant 2 is specified by giving the resultants of two pairs of 
the vectors. For example, we may couple a and 6 to a resultant 
c, then d and e to f before finally adding c and f to give 7. This 
state with total z-component m can be written as 


\(ab)c, (de) f, im) 
being an eigenfunction of the angular momentum operators 

Ae De eC? — (A Bk — (DE) 

I? = (A+B+D+E) = (C+F)? and J,. 
These form a complete set for describing its angular mo- 
mentum properties. 

Clearly there are three ways of choosing pairs from a, 3, d, 
and e, but just as in the three vector case discussed above the 
corresponding eigenfunctions are not independent. They are 
connected by a linear transformation, e.g. 

\(ad)g, (be)h; im) = Yu, [(ab)e, (de)f; im) x 
x ((ab)c, (de)f; i|(ad)g, (be)h; 4). (3.23) 
The transformation coefficient in equation (3.23) that changes 


the coupling defines the 9—j symbol of Wigner [77]. 
((ab)e, (de)f; a|(ad)g, (be)h3 1) = 
abe 
= (een tiner nent fier (3.24) 
hi 


which is identical to the X-function of Fano [28]. For typo- 
graphical convenience it is often written X(abc, def, ght). We 
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may obtain an expression for the 9—7 symbol by expanding 
both sides of equation (3.23) with vector addition coefficients. 
The calculation is exactly similar to that for the W-function in 
equations (3.10)-(3.13), thus we shall not give it here. The result 
expressed as a contraction of six Wigner 3—j7 symbols is given 
in Appendix III together with some related formulae. 

The Racah W-function may be used to contract the sums 
over vector addition coefficients, leading to the form most suit- 
able for numerical evaluation 


X (abc, def, ght) = >, (2k+1) W(aidh; kg) x 
x W(bfhd; ke) W(aibf; ke). (3.25) 
There is often a small number of terms only in the sum over k, 
k being limited by triangular inequalities in the triads (kat), 
(kbf), and (kdh). In particular, if one of the arguments of X is 
zero the sum over & reduces to one term (Appendix IT). 
The orthogonality properties of the 9—j symbol follow from 
the unitary nature of the change of coupling transformation 
(3.23) 


@ 0 C)\e(d.Oac: 
Dor (26+-1)(2f+1)(29+1)(2h+1) «de feyde fi? = 
Galt = 
=> Oce’ Off’ - (3.26) 
The 9—j symbol is highly symmetrical. Interchange of any two 
adjacent rows or columns multiplies it by (—1)? where 
p=a+b+e+d+et+ftg+hti, 
i.e. the sum of all its arguments. It is also invariant under 
reflection about either diagonal. 
The closure relation (3.20) can be used to generate sum rules 


just as in the case of Racah functions in equation (3.21). In 
particular, we use 


((ab)e, (de)f; j|(ad)g, (be)h; j) = 
= 2 a((ab)c, (de)f; j|(ae)k, (bd)l; J) x 
x ((ae)k, (bd)l; j|(ad)g, (be)h; J), 
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to give the analogue of the Racah sum rule tf (3.21) 


abe 
de f> = dei (—1)-- 242k +-1)(2141) x 
ghj 
abc) (adg 
xsjedfrsebh (3.27) 
Uj) lel 


Other such relations are left to Appendix III. 

The large number of arguments makes tabulation of the 9—j 
difficult, but despite this fairly extensive numerical tables are 
now available [67]. The 9—j symbols with two arguments 
equal to 4 are of particular importance. They represent the 
coefficients in the change from L—S to j—j coupling of two 
spin —4 particles and have simple explicit forms. 


3.5. More Complex Coupling Schemes 

The number of possible coupling modes rapidly increases 
with increase in the number 7 of angular momentum vectors. 
Corresponding to transformations between these modes we 
may define more complex coefficients or invariants. Owing to 
the difficulty of numerical tabulation however, their usefulness 
is confined to arguments based on symmetry and orthogonality 
properties. Two 12—j symbols have been defined and their 
properties discussed [24], [39], [47], [61]; they correspond to 
changes of coupling of 5 vectors. 

For n > 4 it can be shown [Sharp 62] that there exist two 
3n —j symbols whose symmetries can be displayed by writing 
their arguments on a cylindrical band, in one case twisted (a 
Mobius strip), in the other not. When 7 is greater than 4 there 
exist other less symmetric symbols as well. Wigner [77] has 
shown that these coefficients are not specific to the group of 
real rotations in space with which we are concerned here, but 
have their analogues for any arbitrary simply reducible group. 


t We have again used the change of phase that occurs when the order of the 


vectors is changed, e.g. 
|(ab)e, (ed)f; %) = (—1)-*-4|(ab)e, (de)f; 4). 


CHAPTER IV 


TENSORS AND TENSOR OPERATORS 


4.1. Scalars and Vectors 

WE are accustomed to deal with certain physical quantities 
such as mass and energy which are in no way connected with 
the orientation of our coordinate system and which have no 
directional properties. Such quantities are scalars, or tensors of 
rank zero.{ 

Other quantities, such as the position of a point in space or 
the velocity of a particle, have associated with them a direc- 
tion as well as a magnitude. These quantities are vectors or 
tensors of rank one and in order to specify them we must give 
the direction as well as the magnitude. Alternatively, a vector 
can be specified by giving its components along the directions 
Ox, Oy, Oz of a set of orthogonal axes. The components repre- 
senting the vector depend upon the particular choice of axes, 
and changing the axes changes the components in a specific 
way. The components (,, 29, 73), or (%,), of the position vector 
r with respect to a new set of axes are obtained from the 
components (%{, x}, 23), or (%;{), referred to the original axes by 
a linear transformation 


oH 2 54%. (4.1) 
J 


The coefficients a,, are the direction cosines of the new axes 
with respect to the old and are definite functions of the Euler 
angles specifying the rotation which takes the old axes into 
the new. If A is an arbitrary vector with components (A,), 
then these components transform under rotation of axes 
according to the same law as the components of r, 


A; = 2 4,4, (4.2) 
A 


t Pseudoscalars, which depend upon the handedness of the coordinate 
system but not its orientation, are counted as scalar for the purpose of this 
discussion. 
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where the coefficients a,; are the same as in (4.1). This transfor- 
mation leads to a new definition of a vector, as a quantity 
represented by three components which transform according 
to (4.2) when the coordinate axes are changed. 

Let us consider a property of a physical system represented 
by a vector A with components (A,;) with respect to some 
chosen set of axes, and enquire how the vector A transforms 
when the system is rotated. The discussion of section 1.4 shows 
that the change in A produced by rotating the system can be 
described by making an opposite rotation of the coordinate 
axes. Thus A transforms to a vector A’ with components 


A; as > a;;A;. 
j 


The transformation matrix d,, 1s the inverse of the matrix a,, 
of (4.2) for the same rotation. 

On either view of a rotation, considering it as a rotation of a 
physical system or as a coordinate transformation, the new 
components of a vector are given in terms of the old by a 
matrix transformation which is specified by the Euler angles 
of the rotation. These matrices form a representation of the 
rotation group of dimension 3, and the components of the 
vector are the basis for the representation. 


4.2. Tensors of Higher Rank 
The mass moments of a system of particles are given by 


M;,=% MN g% aiX a3: 
a 


where (x,;) is the position vector of the particle « and m, is its 
mass. The mass moments are a set of 6 linearly independent 
quantities (Z,, = M,,) and represent a symmetric tensor of 


rank 2. On rotation of axes the particle coordinates transform 
as in equation (4.1) hence the mass moments transform as 


Mi; = > C55, Mj). (4.3) 
kt 


A general (non-symmetric) tensor of rank 2 is represented by 
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9 components 7,,; with the same transformation law on 
rotation of axes as equation (4.3) 


T= 2 tad P er (4.4) 


Similarly the general tensor of rank is represented by 3” 
components which transform according to the obvious 
generalization of equation (4.4). 

If the nine components 7’,; of a general second-rank-tensor 
are written as a column vector, the transformation coefficients 
form a 9X9 matrix; the set of these matrices for all rotations 
forms a representation of the rotation group of dimension 9. 
It is well-known that one can form from the general second- 
rank tensor 7’,, a scalar > T,,,, an antisymmetric tensor 
T= HP Gls 


cP) 


and a symmetric tensor with zero trace, 

de cat UP s+ Py) — P52 Tix: 
Conversely the 5 independent components 7',,, the 3 inde- 
pendent components 7,,, and p T., together specify the 


nine independent components 7,,;, Under a rotation the 


components iT transform amongst themselves as do the com- 
ponents 7,,, while > 7, is invariant. Thus the above reduc- 
k 


tion of the general second-rank tensor 7',, corresponds to the 
reduction of the 9 dimensional representation of the rotation 
group of equation (4.4) to its irreducible components of dimen- 


sion 5, 3, and 1. The tensors is i and > 7’,, are the irreduc- 
k 


ible components of the general second rank tensor belonging to 
the representations 2,, Z,, and Y, of the rotation group. 


4.3. Irreducible Spherical Tensors 

Cartesian tensors are defined in section 4.2 as quantities 
represented by a set of components which have a definite 
transformation law under rotations of the coordinate system. 
It is natural to define a general spherical tensor T, of rank kas 
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a quantity represented by 24 +1 components 7’,, which trans- 
form according to the irreducible representation Y, of the 
rotation group.{ 


ius aa > a DE (apy) (4.5) 


(afy) are the Euler angles (section 2.4) of the rotation taking 
the old, unprimed, axes into the new, primed, axes. Thus (4.5) 
expresses a component 7';, with respect to the new axes in 
terms of the components 7’, defined with respect to the old 
axes. 

It is an immediate consequence of the irreducibility of Z, 
that the tensor T, is irreducible. 


4.4, Products of Tensors 

In the theory of Cartesian tensors the simple (uncontracted) 
product of two tensors of rank m and n respectively yields a 
tensor of tank m-+n. For example, the nine products a,b, of 
the components of two vector a and b are the components of a 
second rank tensor, which transforms according to the repre- 
sentation J, x Z, of the rotation group. The tensor a,b, gives 
rise to irreducible tensors }4(a,b;+a,b,)—36,a.b, aAb and 
a.b, and this reduction corresponds to the reduction 


DXB, = D.4+ PD, +B, 


of the rotation group (equation (2.30)). 

If c is a third different vector, the set of 27 products a,b,c, 
represent a third rank tensor transforming according to the 
representation 2Y,x2,xQD, of the rotation group. This 
representation is reducible as follows, 


De i — Pex De), 
that is, the 27 components of the tensor a,b,c, give rise to one 
irreducible tensor of rank 3, two of rank 2, three vectors and 


¢ The definition (4.5) is chosen to agree with the transformation law (2.24) 
for spherical harmonics on rotation of axes. This means that spherical har- 
monies provide a special example of spherical tensors. 
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one scalar. Explicitly the three vectors are (a. b)c, (b.- c)a, 
(c.a)b, and the scalar is (aA b).c. In this way we can count 
the number of different irreducible parts of a general cartesian 
tensor of rank n. 

If R,, and S,, are irreducible tensors of rank k and k’ re- 
spectively in the spherical representation, the (2k-+1)(2k’ +1) 
products R,,S,7 form a tensor transforming under the repre- 
sentation Y, xQ,, of the rotation group. This representation 
is reducible (equation (2.30)) and its reduction gives the irre- 


ducible tensors{ 


To (kk’) =5 ReiSyq ek‘ 9q'|KQ), (4.6) 


with K running from k+k’ to |k—k’| and Q = q+¢q’. 

Unless the operators R, and S, commute there is no simple 
relation between 7',9(F,, S,) and T9(S,, R,). An interesting 
special case occurs when R = S and k = k’ which is analogous 
to the exchange symmetry property of two-particle states 
(section 2.7.3). Provided the components R,, and Ryo_, 
commute interchanging them merely multiplies T-g(kk) by 
(—1)*-* (= (—1)* if & is integral). Thus only product 
tensors with K even do not vanish. This property is a generali- 
zation of the vector relation vA v = 0.8 

When k = k’ and K = 0 the product (4.6) is the generali- 
zation of the scalar product of two vectors. When k is integral 
another phase and normalization is used in the definition of 
this mae 


R,.S, => (—1)*R,8,_, = (—-)) ¥/ (2k +1) Ty(R, S)- (4.7) 


4.5. Tensor Operators 
The notion of tensors can be extended directly to quantum 
mechanical tensor operators. An irreducible tensor operator 
{ Sometimes we use as arguments of a composite tensor the actual tensors 
out of which it is constructed writing it as 
TrQ(Fi, Sy) instead of Tro(k, k’). 


§ The angular momentum J does not obey this relation but rather 
Ja J = iJ, because its components do not commute. 
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of rank k is an operator with 2k +1: components Tq Which 
transform under a rotation («fy) of axes ast 


Lae = DA lire = > Ne ASHE (4.8) 
Pp 


Relation (4.8) with D representing an infinitesimal rotation 
leads to commutation laws of the components of the angular 
momentum J with the components of T,. 

Let D be the infinitesimal rotation (1—iaJ,). Thus from 
equation (2.15) 

DE = (kp|1 —iaJd,|kq) = 6,,—ia(kp|J,|kq). 
With this rotation equation (4.8) becomes 
(L—tod,)T (1 +tad,) = 2 IN ee 
Pp 
or 


Beate. shy = be T,, (kp|J,|kq). 
P 


Putting J, = J, and J, in turn and using the matrix elements 
of J, from equations (2.4) one obtains 


(J, Pq] = oat. 
[Fa Tool = [LG +IMEF OTT py (4.9) 
The commutation rules (4.9) of J with spherical tensor 
operators can be used in finding the spherical equivalents of 
Cartesian tensors. For example, if A is a vector, 
[J., A,] — 0 
from the commutation relations in equations (1.9). Thus 
A, = A,. From (4.9) we have 


A ean 


1 
+1 = val 
1 
a/2 


I 


+ A rotation R of a quantum system transforms an operator A as 
A’ = Dt(R) AD (RB), cf. section 1.5. 


However, a rotation R of axes is equivalent to an inverse rotation R- of the 


[J., A,]+aJ,, 4,1, (4.10) 


(A,+14,). 
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Again, if A and B are vectors the spherical tensor of rank 2 
formed by their product is 


GAB), = > A,B lmij2g), 


mn 


In particular, using the value of (22|1111) = 1 
T(AB)o2 = A,B,. 
Then from the commutation relations (4.9) 


T(AB), = el A,B,] 
= #((J_, AB, +A,{J_, By) 


1 
= (A, Beane) (4.11) 
1 
and T(AB)oy = 75 (34,B,—A . B). 


4.6. Spherical Harmonics as Tensors 

Spherical harmonics transform on rotation of axes according 
to equation (4.5), hence are examples of spherical tensors. 
There are some relations involving products of spherical 
harmonics which follow simply from section 4.4. These take 
their simplest form when expressed in terms of the modified 


spherical harmonic 
47 | 
Cn = (a5) Vig 


If C,,, (8, p) and C,,(9’, $’) are modified spherical harmonics then 
the scalar product (equation (4.7)) 
C, is C,, = (1) 0 g)C,_ (6, Y’); 
q 
is invariant with respect to rotation of axes. It follows that the 
scalar product must be a function of the angle © between the 
directions (6, ) and (6’, ¢’), this angle being the only quantity 


system so that the transformation operator in equation (4.8) is D(R-1) or 
D+(R) rather than D(R). If |0) is a spherically symmetric wave function 
T',q |0) is a wave function with angular momentum quantum numbers (f, q). 
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independent of choice of axes. Choosing axes so that the 
direction (6’, ¢’) becomes the new z-axis 


Creal’, $’) > C,4(0, 0) = 6(¢, 0), (cf. equation (2.10)) 
C.o(9, ) > P,(cos 6) = P,(cos ©), (ef. equation (2.11)) 

and the scalar product (cf. equation (2.25)) 
CG, 2—oP (cos ©). (4.12) 


Again, if C,,(6,p) and C,,,.(8,p) are spherical harmonics of 
the same angles (6, ») then 


¥ (KQ|KR' a0 Oy Oey 
aq 


is a tensor of rank K. Since it is a function of (6, m) only it must 
be proportional to C,9(6, 4). That is 


2 (KQ|KK'qq') Cyq(9, $)Oy(9, 6) = AgCxg(, @); 
= C4) 


kk’00) Cico(9, 4). (4.13) 


The value of A, in equation (4.13) is found by putting 9= d= 0 
and making use of equation (2.10). Equation (4.13) follows also 
as a special case of the rule (2.31) for the combination of 
rotation matrices. 

Examples of the spherical harmonics polarized by the vec- 
tor operators V and L are discussed in section 4.10.2 in 
relation to vector fields. 

Spherical harmonics with different arguments may also be 
used to construct product tensors 


Bxo(bh’) = ¥ (KQ|kk'qq') Cp(W)Cya(w). (4.14) 


These tensors, which we may call bipolar harmonics, con- 
tinually appear in problems involving two directions, 
(u and u’ being unit vectors along those directions). For 
example, they describe the angle dependence of two particles 
moving in a central field in an eigenstate of total angular 
momentum JL, 

(u,u’|ll'LM) = [(20-+1)(20’ +1)]}} Bz y(ll’)/47. 
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Often they are required in problems of the angular correlation 
of nuclear radiations [5], [16]. 
The bipolar harmonics are orthogonal for integration over 
the angles of u and u’ 
6(1,1%) 6(Ip1) 6(LL’) 6( MM’) 
'IR* yey 2 SS 
[auf au BE yal) Bry (Gh) = 162 (1, +1)(2l,41) 


(they would also be normalized to unity if they had been 
defined with the normalized Y,, rather than with the C,,) and 
have the closure property 


> |Bru(hl.)|? = 1. 
LM 
When K =Q = 0, (4.14) reduces to the spherical harmonic 
addition theorem (4.12) 
Boo( kk’) = (—)F 8(kk")P, (a. u')/4/(2k-+1). 


Again, when k = k’ = 1, we retrieve the vector relations 
Bo (11) = —u. w’//3 
and By (a Aw Fi 2 


4.7. Matrix Elements of Tensor Operators 


4.7.1. Wigner-Eckart Theorem 

We wish to evaluate matrix elements of a tensor operator 
with respect to the state vectors of a dynamical system. When 
angular momentum is conserved, so that the state vectors are 
eigenfunctions of J? and J, the matrix elements of a tensor 
operator have a simple geometrical dependence on the 
magnetic quantum numbers. Let 7',, be a tensor operator of 
rank k and consider the matrix element («J M|7,,|«'J'M’). 
The state vector 7',,|«’J’M’) transforms according to the 
representation J, x Y,. of the rotation group. We reduce this 
representation to its irreducible components by forming the 
vectors with angular momentum (K, Q) 


|BKQ) = > (KQ|I'EM'q)T yal’! M’). 
q. 
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Inverting this relation and taking matrix elements with 
(2J M| gives 


(oJ MT yqle’J'M’) = > (uJ M|BKQ)(KQ|J'kM'q) 
KQ 


= (at M|pJM)(JM|J'k.M'q) 

= (—1)* (ad ||T,||o'F’) (J M|J'kM'q). (4.15) 
The inner product of state vectors («JI|BJ.M) is independent 
of M. Thus the matrix element of the tensor operator factorizes 
into two parts. The directional properties are contained in the 
Clebsch-Gordan cocfficients and the dynamics of the system 
appear only in the scalar matrix element («J M|8JM) usually 
written as («J||T,||«’J’) and called the reduced matrix 
element.{ 

Inverting (4.15) gives the analogue to the scalar contraction 

of vectors (section 3.2) 


(aJ||T,|Ja’F’) = (—D*® YS IM|I'EM'ghoI M| To I'M’). 
M’q@ 


The prototype of such matrix elements is the integration 
over three spherical harmonics: 


(lm| Vy yl’) = § Yn (8, $)* Yer(9: P) Yoml9, p)dQ. (4.16) 
This is evaluated by combining the last two harmonics in the 
way described in the last section (equation (4.13)) 

Via Yeu = > Yael, p)kq| LV Msn") (ke0| L100) x 

Pera y(2b+1y} 
| 47(2k-+-1) | 

From the orthogonality properties of the harmonics we see 

that only the k = 1, m = q term has a non-zero integral. Thus 

(ln| Yz4;|l'm’) = (lm|l’ Lm‘ M\X]]| LA 

’ $ 

C]| Y;||2) [ae] (lO|L1'00). (4.17) 


t The definition used here is that of Wigner [78] and Rose [54]. Wigner’s 
tensors T¥,, however, are equivalent to the adjoints, T*, of ours (cf. section 
4.8) because of the rotation conventions he adopts (Appendix IV). Racah [48], 
[31], and Edmonds [22] define reduced elements which are (2J +1)? times those 
in (4.15). (The factor (—)?* makes the phases of the two definitions identical.) 


where 
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The physical significance of the Wigner-Eckart theorem is 
now clear. The ‘integrand’ of the matrix element on the left 
of (4.16) must be a scalar if the element is not to vanish. So the 
only part of the product of 7’',, and |J’M’) which can contri- 
bute is that which rotates contravariantly to (JJ/|; that is, 
like |JM). By definition its amplitude is just the Clebsch— 
Gordan coefficient, and this is the geometrical factor contain- 
ing all dependence on magnetic quantum numbers. A further 
important example is the matrix of the angular momentum 
operator J whose spherical components J, form a tensor of 
rank 1 as discussed in section 4.5, The reduced matrix is easily 
obtained by considering J, = J,: 


(IMJ ,|J’M’) = Mo(JJ') 5( MM’) 
= JM|J'1M'0)(||J| IJ’). 
The explicit form for the Clebsch-Gordan coefficient (Table 3) 
gives immediately 


(J'||J|| 7) = 6(JI’) / (J +1)}. (4.18) 


All processes which do not involve a definite spatial direc- 
tion (such as radiation transition probabilities) are indepen- 
dent of the magnetic quantum numbers and are described by 
reduced matrix elements. This can be seen explicitly by 
supposing that matrix element (JU|T.|F'M’) describes a 
radiative transition from state |JJ/) to |J’M’). If there is no 
preferred direction then all orientations JM’ of the final state 
are equivalent and the total transition probability is 


IM [Ped M)P 
a 
= [APOE orale aay 
Sy 
= |All? 


since the Clebsch-Gordan coefficients are normalized. 

The above theorem due to Wigner and Eckart [78], [21] 
represents an important step in the program outlined in the 
introduction to make a division between the geometrical 
and physical properties of a system, or to remove from the 
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problem those aspects which relate to the symmetry of the 
system. The Wigner-Eckart theorem makes this division 
explicit, the matrix elements factorizing into a Clebsch-Gordan 
coefficient containing purely geometrical information concern- 
ing the orientation of the system and a reduced matrix 
element depending on the detailed physical structure of the 
system. 


4.7.2. Selection Rules 

Selection rules connected with the conservation of angular 
momentum are mostly derived by using the Wigner—Eckart 
theorem. It follows from the properties of the Clebsch-Gordan 
coefficient in section 2.7 that the matrix element 


(IM|T,,|J’M’) 


vanishes unless 


M=q+M' 
and the triangular conditions are satisfied, 
\J—J'|] <k < JJ". (4.19) 


Clearly for 4 — 0 we must have / = J’, M = 7 and’ the 
corresponding Clebsch—Gordan coefficient is unity. This shows 
that the matrix element of a scalar is independent of the 
magnetic quantum number (i.e. of the choice of quantization 
axis), as it must be, and is identical to the reduced matrix 
element with our definition 


(JM|P ol’ M') = o(FJ') 6M’) (|| To[U’).- 


The spherical harmonic integra] (4.16) also contains the parity 
selection rule for the spatial part of a single particle matrix 
element; (/+/'+L) must be even, otherwise (11'00|L0) 
vanishes (cf. equation (2.35)). 


4.8. The Adjoint of a Tensor Operator 
The hermitian conjugate or adjoint 7’+ of an operator T is 
defined by expressing its matrix elements in terms of matrix 


elements of 7' as q|r+|2) = (2)P]1)*. (4.20) 
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The operator 7 is called hermitian if 7 = T+. A necessary 
and sufficient condition for an operator to be hermitian is that 
all its eigenvalues should be real. For example, the compo- 
nents of J,, J,, J, of the angular momentum operator are 
hermitian and from section 2.2 they have real eigenvalues. 
If A and B are operators and a is a complex number it fol- 
lows from the definition (4.20) of the adjoint operator that 


(AB)+t = BtAt+ and (aA) =a*At (4.21) 
hence 
LAR BB eA) eae 


From (4.21) and from the hermitian property of the compo- 
nents of J it follows that 

Cl Sai 
and from the commutation relations (4.9) of the shift operators 


Js and Jo with the components of a spherical tensor operator 
T; that 


[For ie = leh ss as —qT;, 
(Ja, Tig) = Je, Pl” = bFQt KAO Tae 


Thus T; transforms contragradiently to T, (section 2.6). 
Hence if we define an operator 7’, by 


a (4.22) 


then the 2k+1 components of ie transform as a tensor 
operator of rank k, (p is an arbitrary integer or half integer 
depending on whether k is integral or half integral). 

As a result of the property (4.22) at most only the q = 0 
component of a tensor operator can be hermitian. It is possible, 
however, to extend the notion of hermitian operators to cover 
tensor operators by defining a hermitian tensor operator as 
one with the property 

Tin ae Te (4.23) 


The choice of phase p is arbitrary, although some authors 
demand that p =k (Edmonds [22]) and others use 
p = 0 (Schwinger [60]). Then operators which have p = k+1 
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or 1 respectively are called anti-hermitian. On the other hand 
the phase (—1)* is essential to preserve the correct rotational 
properties (4.8). The spherical harmonics Y, and the spherical 


components Jg of J (Jo = Jz, Jar = FS 
(4.10)) show this hermitian property with p = 0. 
Yiu a Yiu a (oot) Ye 


Je = (1 


(Jz+tJy), ef. equation 


There is a simple conjugation property for reduced matrix 
elements of a hermitian tensor operator which corresponds to 
equation (4.20). From equations (4.22) and (4.20) 


(JM|T,,|J'M’) = (-1P*UM|Tt_|J' mM’) 
= (—1)-*'M'|T,_|FM)*, 
or in terms of reduced matrix elements 
VJ (2F +1) J||7,||7) 
= (—1 7? Vf (2J' $1) "||P, ||7>*. (4.24) 


4.9. Time Reversal 
The results of section 1.8 do not apply directly to a system 

with angular momentum because J is not invariant but 

changes sign on time reversal. If @ is the time reversal 

operator then 

P 07,071 = —J 


z 


OO ee) 
(The anti-linear character of @ requires that the complex con- 
jugate of any number which appears in the transformed 
expressions should be taken, hence J, — —J_.) It follows 
from the above transformation equations that if |jm) is a set of 


states with angular momentum (j, m) then 
J, Olajm) = —OJ,|ajm) = —mBlajm) 
J, 8lajm) = —OJ,|ajm) 
= —[GFm+l)(jtm)} OlajmF 1). (4.25) 
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Thus @|ojm) transforms contragradiently to |ajm). If we put 
|Bjm) = (—1)?-"6|aj—m) then |fjm) has the correct trans- 
formation laws for a wave function of angular momentum 
(j, m). 

We consider a system with Hamiltonian H which is time 
reversal invariant. The requirement that wave-functions are 
eigen-states of angular momentum makes it impossible to 
choose them invariant under time reversal, but the trans- 
formation properties of a time-reversed wave function 
discussed in the last paragraph make it possible to choose a set 
of eigenstates of H, J? and J, for which 


Blajm) = (—1)?-™|aj —m). (4.26) 


The phase 7 is so far not specified but it is usual to require that 
p =J, for then the choice of phases of the vector addition 
coefficients (section 2.3) determines that 
|7 M1) = > jij ot») (mm,|J M1) 

LOG 
satisfies equation (4.26) if it is satisfied by |j,m,) and j.m,). If 
we say that a tensor operator T, has a definite transformation 
under time reversal when 


e7',.0-) = (—1)7,_, (4.27) 


then by an argument similar to that given in equations (1.13) 
and (1.14) it can be shown that the reduced matrix elements of 
T,, are real if (p —k) is even and pure imaginary if (p —k) is odd. 
In each case only one real parameter is necessary to specify 
the value of the matrix element. 

The requirement that (4.26) and (4.27) should be satisfied 
with p = 1 is the origin of the factor i‘ sometimes introduced 
into the definition of spherical harmonics Y,,, (section 2.3). 

Equation (4.26) implies that 


8? |oxjan) = (—1)”” |axjm) = (—1)* |ajm) 
= |ajm) if j is an integer 


= —|ojm) if j is a half-integer. 
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4.10. Multipole Expansions 


4.10.1. Scalar Fields 
Let V(r, 8, y) be an arbitrary scalar field. Rotating the field 


through a small angle « about the z-axis produces a new 
field V’ given by 


V(r, 6, ~) = V(r, 9, p—«) 
a 
= (1 -155) V 
= (l—iaL,)V. (4.28) 
mC oa 2a 
The operator L, = as = -i{xs ay — tee is the infinitesi- 


mal rotation operator for the field about the z-axis. The 
remaining components of L, L, and L, are defined similarly. 

A multipole expansion aims at expressing the field V as a 
sum of components with rotational properties of spherical 
tensors, and for a scalar field this is achieved by expanding V 
as a series of spherical harmonics: 


v( (7, 6, @) =2 Vint) tO, ¢)* (4.29) 


where Vive [Tint 6, 6) V(r, 6, 6) dQ. 


A rotation of the field V(r, 6,4) specified by Euler angles 
(«, 8, vy) is given by transforming Y,,, according to equation 
(2.22) while a transformation of the field produced by a 
rotation of axes through angles («, 6, y) is given by trans- 
forming the coefficients V,,, according to equation (4.5). A 
simultaneous transformation of Y,,, and V,,, corresponds to 
a rotation of axes with the field and produces no change in V. 
The components of V,,, transform as components of a tensor 
of rank / and are the multipole components of the field V. 

As examples we consider the expansion of a plane wave and 
of a potential function. A plane wave travelling in the z-direc- 
tion is symmetrical about the z-axis and can be expanded as a 
series of Legendre polynomials referred to this axis, 


ett? — Y 7'(2141)j,(kr)P,(cos 8). (4.30) 
I 


64 TENSORS AND TENSOR OPERATORS IV, § 4.10 


The radial functions j,(kr) can be expressed as Bessel functions 
of half-odd integer order 


jkr) = (57. ).F att 


We may obtain the equation for a plane wave travelling in an 
arbitrary direction (8, «) by rotating the wave travelling in 
the z-direction. Equations (4.30) and (4.5) (or alternatively the 
addition theorem (4.12) for spherical harmonics) give 
exptk.r = > 714(2141)j,(kr)C,,,(8, ACE (B, «). (4.31) 


im 
The electrostatic field produced by a system of charges 
provides an example of a scalar potential field. If r andr’ are 
two vectors with direction (6, y) and (6’, ¢’) and if © is the 
angle between the two vectors then for r > r’ 


1 Eee 
jr—r'| = yaa (cos Q) 


l 
re : 
7 » AFC im(9, P)Cim(8, $'). (4.32) 
im 


Thus the potential at r due to a charge distribution of density 
p(r’) is given by 


a "<eee $), (4.33) 


ms 


-(< 
where 


ee 7 Cin (8's 4") d 


The tensors Q,,, are the multipole moments of the charge 
distribution. For a system of point charges 


Qin an >» er Crm(9;, $;)s (4.34) 


where (7,, 6;,¢,) are the spherical coordinates of the :th particle 
and é, is its charge. For a quantum system the /th multipole 
moment of a state is the expectation value of Q,,, in that state. 


| 


| 
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4.10.2. Vector Fields 
Let A(r, 0,¢) = >» Ae; be a vector field. The e, are unit 


a 
vectors along the coordinate axes and the A, are the com- 
ponents of A along those axes. A rotation of the field through 
a small angle « about the z-axis produces a new field 


A'(r, 8, @) = Ba (r, 0,6— 
where e; are unit vectors along a set of axes rotated with the 
field. These unit vectors expressed in terms of e, are 
e, =e,+ae,,e€, =e,—ae, and ee, 
or e, = €; +a, A e,. 


Expanding the A, as a power series in « and substituting for 
e, in terms of e, we obtain an expression for the rotated field to 
the first order in « 


oA 
A’(r, 8, ¢) =A+a(e, AA-— si) 


og 
or Ge ine) = (ketal Wake (4.35) 
0 
where 1 an vA; 
7] 7] 
= =i a -v-) +1e, A, (4.36) 
= 1s 


Thus one part of the change produced by rotating a vector 
field is due to the variation of the field components at different 
field points. The differential operator L, generates this part of 
the transformation. The other part of the change generated by 
S, is due to a re-resolution of the vector field components 
when the field is rotated. The operator S, can be written as a 
3x3 matrix and is one of the infinitesimal transformation 
matrices of the representation Y, of the rotation group, thus 
the ‘intrinsic spin’ S = 1 properties associated with a vector 
field. 
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Let us define a set of three vector fields from the unit 
vectors e,, e,, and e, by 


@) = en (4.37) 


S, = 1e, A. (4.38) 


The operators (4.38) obey the commutation relations of angular 
momentum operators, and, acting on the set (4.37) of vector 
fields they transform them as components of a tensor of rank 
one, as may be proved by testing the relations (4.9). For 
example 


S,e) = (e, +ie,) Ae, = —(e,+te,) = 24a. 


The field r provides an interesting example of the above 
analysis. 


J,v cam ay FS (ee yey 2@7)) 
Fam atl Jee +xe,) +4e, A (xe, +ye, +2e,), 


from the definition (4.36) of L, and Sz. 
Hence it follows that 
oot 0, 


a result which is not surprising because the field r is spherically 
symmetrical and a rotation transforms it into itself. 

Vector fields which are generalizations of spherical harmonics 
may be formed by taking products of the vectors e, defined in 
(4.37) with spherical harmonics and by using equation (4.6) 
to give an irreducible tensor. 


NG == OD, UES) On (4.39) 
The set of 2L+1 fields Y}}, transform amongst themselves 


as components of a tensor rank J, i.e. according to equation 
(4.5). They are products of the tensors e, and Y,,, of rank 1 and 
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l respectively. As an example the vector field r can be written 


47s 
as rere > (0 iti 


= —14/4r Y51,(09). 
Thus the field r is the scalar product of the spherical tensors 
rY,,, and e,, (both of rank 1) and is invariant for rotations as 
concluded above. In addition the vector spherical harmonics 
Y¥ (9) satisfy the orthogonality conditions, 
[V2 (06))*. V3 (09) dQ = dU’) (LL') (MM), (4.40) 


and they form a complete set for expanding the angular 
dependence of an arbitrary vector field. 

It has been shown that the vector field r is invariant for 
rotations of axes. If ¢,,, are a set of 2L-+1 scalar fields 
forming a spherical tensor of rank L (e.g. Y,,,) then the set of 
vector fields r¢,,, also transforms as a spherical tensor of rank 
L. The vector field operators 


= 0 0 
~ eae TPH ay Toray’ 
L= -ir,V, 


and VAL, 


share with the field r the property of being invariant under 
rotation. Hence each of the three sets of vector fields, 


Vor Ldpy and Valdzy, (4.41) 
form spherical tensors of rank L. In the special case when 
=e 
Pie = Su L Ying = {L(L +1] Y%q 


is a vector spherical harmonic. The fields (4.41) share with the 
vector spherical harmonics the property of orthogonality with 
respect to angular integrations and sometimes are a more 
convenient set than the vector spherical harmonics for making 
a multipole expansion of a vector field. 

Such a case arises in considering solutions of the wave 
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equation for a vector field A with wave number k, 
VA+A = 0. (4.42) 

The functions ¢,,, = i2(2L4+1)5,(kr)C,4(9¢) where j,(kr) are 
the spherical Bessel functions introduced in section 4.10.1 
satisfy the scalar wave equation corresponding to equation 
(4.42), and with this choice of 4,,, the vector fields (4.41) form 
a complete set of expanding solutions of equation (4.42). In 
addition the set A;y,, Az» Atm With 

Ary = (th) Very 

Aba = [AV{L(L +) 17 V ALdry 

Atu = (V{L(L4+1)})7* Udi (4.43) 
have the same normalization as ¢,,. The fields A,,, are 
irrotational and have parity (—1)”, while the fields A%,, and 
A”,, are solenoidal and have parity (—1)” and (—1)7*? 
respectively. The notation looks forward to the application of 
these fields in expanding the vector potential of the electro- 
magnetic field. The fields A,,, are longitudinal fields whereas 
the fields A‘, and AZ, will represent the electric and magnetic 


components of the transverse field. It follows from equations 
(4.43) and (4.42) that 
Aba = EUV AAT) (4.44) 
and INE tes POY Io he 
As an example of the expansion properties of the vector 
spherical harmonics and of the vector multipole fields (4.41) 
we obtain expansions of a plane wave in a vector field moving 
in the z-direction and with polarization vector e,. If g = 0 
the wave has longitudinal polarization, if g = +1 the wave is 
transverse and has right or left circular polarization. Remem- 
bering the expansion of a scalar plane wave given in equation 
(4.30) and the definition (4.39) of vector spherical harmonics 
we get an equation 
ee = 2 Uy/{4r(20+ 1) hi (kr) Yip(9.b)e,, 
= SEV E4a(A+M Chr) Mg |LOG, (4.45) 


as a series of vector spherical harmonics. 
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Alternatively we may expand the plane wave in terms of the 
fields (4.43). For the longitudinal component 


= Aj. (4.46) 
L 


The transverse wave is a little more complicated. The sym- 
metry about the z-axis determines that e,e** is a sum of 
components with M = q = +land since the fields (4.43) form 
a complete set 


Taking the divergence of both sides of this equation deter- 
mines that a; = 0, taking the curl of both sides (V A) and 
using equations (4.37) and (4.44) gives b; = qc, and, finally, 
taking the scalar product of the operator L with both sides and 
using the expansion (4.30) for exp (tkz) gives ¢, = —1/,/2. 
Hence the transverse circularly polarized wave has expansion 


ee = —, > (GAL +A): (4.47) 
L 


Equations (4.45) and (4.47) give the expansion of a plane polar- 
ized vector wave travelling in the direction of the z-axis. The 
expansions for a wave travelling in an arbitrary direction 
with wave number k can be obtained by rotating these 
expansions. For example the transverse wave has the expan- 
sion 
ike 1 m é L 
eet! = ==> (Gary tAS yD (Rk) (4.48) 
V2 tt 

where fF is the rotation taking the z-axis to the direction k 
(equation (2.22)). These calculations show that the Azm are the 
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natural fields for expanding a longitudinal vector wave. In the 
expansion of the vector potential in a transverse electro- 
magnetic wave the components Aj, and AT, correspond 
directly to its division into electric and magnetic multipoles. 


4.10.3. Spinor Fields 
The general spin-half field has two components and may be 


written Srv) 


where y, is the two-component spinor, 


nal’) ot c=() 


The two components S, and S_, ofS are defined with respect 
to a particular z-axis. A rotation changes the spatial depen- 
dence of the S, and also mixes the two components. Just as 
in the case of a vector field the infinitesimal rotation operator 
is a sum of t art 

is a su wo parts Fares 


where L, acts on the spacial dependence of the field and S;, 
mixes its components. For the spin half-field the spin operators 
can be represented by the Pauli matrices S; = 4e,. 

The analogues of the vector spherical harmonics are 


gin = > (jn|tlom)7,i'C im( 8). 


Note that j is always half integral and can never be zero. 
Hence there is no spherically symmetric spinor field analogous 
to the vector field r. 

There is also an expansion of a spinor plane wave similar to 
(4.45) 


Loee® = FY (2+1)(4100| jo) 75,(kr) by, D2,( B). 


ijn 


The suffix o refers to the polarization of the wave along k, its 
direction of propagation. 
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4.10.4. Electromagnetic Multipoles , 
The electric and magnetic fields are given in terms of a 
vector potential A obeying the gauge condition div A = 0 by 


1 0A 
ial Sentgl 73) 10, = a 
c Ot 


The component of the vector potential which oscillates with 
wave number & = ~ can be expanded in terms of the multi- 
pole fields (4.43) as 
Ay Ae Am 
= X Qi Aia + 9a AT) (4.49) 


The fields V¢,,, do not obey the gauge condition div A = 0, 
therefore do not appear in the expansion. The first terms 
of expression (4.49) denoted by A? are the electric multipoles 
and the terms denoted by A” are the magnetic multipoles. 

On quantizing the electromagnetic field the Fourier coeffi- 
cients 95, and 74, in (4.49) can be written in terms of creation 
and annihilation operators for photons (Heitler [34]). 

If the photon states are normalized in a spherical box of 
radius #, then the electric multipole part of the vector poten- 
tial takes the form 


ho\} 
Att) => (FZ) Able) Heeb AG(D)) — (4.50) 
LM 
where a¢3, and a%,, are the creation and annihilation opera- 
tors of electric 2”-pole photons.{ There is a similar expansion 
for the magnetic multipole part of the field. 


+ An electromagnetic transition probability is given by 
2 


where p(Z) is the density of final photon states. For waves of a definite angular 
momentum and for a spherical box of radius R, p(#) = R/whc. The radius of 
R of the region of quantization therefore cancels with the factor R in A (equa- 
tion (4.50). 


6 
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4.10.5. Multipole Sources and Hyperfine Interactions 

In potential problems where the extension of the source of a 
field may be neglected it is convenient to introduce the idea 
of a point source, for example a point charge or a point 
dipole in an electric field. Mathematically the source of the 
field produced by a unit point charge is represented by a 
6-functiont and the potential equation for the field 1/r is 

ve = —47 o(r). 
r 

If we differentiate both sides of this equation with respect to 
z we obtain the result.that the source for the dipole field 
Zo 06(r) 


ols 4a De 


This idea of a point source can be generalized by introducing 


1 
a source function for the multipole field = iC rm (6¢) by the 
equation 


Il 
We (<aCiu(09)) = —47n 6,,,(r). (4.51) 


The source functions 6,,,(r)$ are tensors of rank L, and 
are generalizations of the three-dimensional 6-function which 
vanish if r ~ 0 and have a singularity at r = 0 such that if 
f(r) is an arbitrary scalar function regular at the origin 


[f #0) duadt)dr = Fhy (4.52) 


where F',,, is the coefficient - r”C,4 in a power series expan- 
sion for f(r) about the origin. 
In particular 


| r’C*, 4 Ony(F) dr = d(L'L) (MM). 


t A 6-function is defined by the property that for any function f(x) regular 
atz = 0 


[re S(x) da = f(0) 
and the derivative of a 6-function 6’(z) by 
dé a 
ce \dz = — | 6(2) 2 a = —f’(0) 


a result obtained by a formal partial integration. 


§ 6,4 is an Lth derivative of the 3-dimensional 6-function. 


IV,§ 4.10 TENSORS AND TENSOR OPERATORS 73 


The set of singular functions may be further augmented by 
introducing functions 


; (m+ DI (2E+1)) 
Ona = ae Ons th b(t) (4.53) 


and the resultant set includes all derivatives of the 6-function. 
The functions 67, have the property that for a scalar function 
f(r) 

[f*cr) T) Oy (0) dr = PThy (4.54) 


where F,, is the coefficient of r”t°"C,,, in a power series 
expansion of f(r) about the origin. This result can be demon- 
strated by a series of formal partial integrations, reducing 
equation (4.54) to an integral of the form (4.52). 

If we have an arbitrary scalar density distribution p(r) 
concentrated near the origin we may introduce the formal 
series 


(r) = 2 Qin Oiu(t) (4.55) 
with coefficients Q7,,, given by 


bn = |e Onn: (4.56) 


Equation (4.55) can be considered as an expansion of the 
function p(r) as a series of derivatives of the 6-function in the 
following sense. If ¢(r) is a function of r which is expandable in 
a power series about the origin with a sphere of convergence 
which includes p(r) then 


[d(r)p(r) dr ={9(r)p'(r) dr (4.57) 


where p’(r) is the series given in equation (4.55). The relation 
4.57 may be proved by substituting the series for p’(r) and 
using the property of the 6-functions given in equation (4.54). 
It may seem strange at first that the function p(r) of finite 
extent may be represented as a series of 6-functions which take 
values only at the origin; but if it is remembered that equation 
(4.57) defines the meaning of the series p’(r) and that ¢(r) must 
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be of finite extent then the difficulty is removed. The expres- 
sion (4.55) may be used to find the potential V(r) produced 
by a charge distribution p(r). We have 


V2V = —4rp 
and combining equations (4.51), (4.53), (4.55) we get 
Ci = 6? V(r) 
as fa IM = pk IM ; 4. 
V(r) » Keres » dar LM In (In +2L 41) (4.58) 


IM p=l 
The first terms in (4.58) give the potential outside the charge 
distribution (equation (4.33)), while the further terms represent 
the modification of the potential inside the charge distribution. 
The interaction energy between an extended charge 
distribution p,(r) and a concentrated distribution p,(r) (e.g. 
the potential energy of an electron in the field of a nucleus) 
can be found by expanding p,(r) as in equation (4.55) and 
obtaining V(r) by equation (4.58). The interaction energy is 


W = | or) V,(r) dr (4.59) 
=) OF u(eOza(), 
IM 


retaining only the first term in the expansion (4.58). The 
Qu(n) are given by equation (4.56) and 


Oat 
Qzule) = | se dr. (4.60) 


Inclusion of higher terms in equation (4.58) gives an expansion 
of the interaction energy W as a power series in radius of the 
charge distribution p,(r). 

A similar expansion can be made for the magnetic field 
produced by a static magnetization distribution M and a 
stationary current distribution of density j. The magnetic 
induction B and the magnetic field H are given by 


B = H+47M, 
curl H = 47j, 
B = curl A, 
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and the vector potential A is chosen,to satisfy the gauge con- 
dition div A = 0. These equations yield an equation for the 
vector potential 


WA = —4n(j+curl M), 
= —47J, 
introducing an effective current density J which is the sum of 
the charge current and the ‘magnetization current’. Remem- 


bering the stationary current condition div J = 0, J may be 
written in terms of two scalar fields as 


J = L¢,+(V A L)Q,. 
If we retain only the lowest terms of a power series in the 
radius of the charge distribution it can be shown that the 
field ¢, gives no contribution. Expanding ¢, as a dé-function 
series (equation (4.55)) and keeping only the lowest terms for 
each multipole we can write 


1 * 
J=> ap ld Ona (4.61) 
aT a 
hence 


1 
L.J => 7 ial* Oyy(0) 


LM 
= td > Mee (Lal) Ope lt). 
iM 
Therefore the coefficients M/,,, are given with the help of 


equation (4.56) as 


iy = an { O14 L.J dr. (4.62) 
Since J = j5+V ,M we see that the multipole moments M,,, 
have contributions from the current j and the magnetization 
M. The current contribution has exactly the form of equation 
(4.62) with j replacing J, while the magnetization part can be 
transformed by partial integration. Using the result of equa- 
tion (6.6) (ef. Appendix VI) we get 


Di gerre| Mec 5) Madr. (4.63) 
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In a quantum system the current j is due to motion of charges 
while the magnetization M is associated with particle spins. 
For a single particle the current is represented by the operator 


° é oe UV 
Ve ee a a IL 


and the magnetization by 


eh 
M = oe Is, 

where m is the particle mass and g, and gg are its orbital 
and spin g-factors, (eg, is the particle charge). In a many 
particle system j and M are sums of the corresponding single 
particle operators. With this identification the moments (4.62) 
and (4.63) are exactly equivalent to the dynamic magnetic 
moments (6.12b) and (6.12d). 

Given the magnetic moments M,,, of the source the vector 
potential A can be calculated by using relations (4.51) and 
(4.61) giving 


a 
A(r) => eM Lr OPCy(8, 9). (4.64) 
IM 
The induction B is found by evaluating curl A and making use 


of a result analogous to equation (6.6) 
(V ALG OC, LV eG ea rircoeee 
giving 


; 4a 
B => re huge Cine) sae Sua} (4.65) 


IM 
We shall show that the 6-function singularity in B at the 
origin is responsible for the spin contact interaction in hyper- 
fine structure. 

The Hamiltonian for the magnetic interaction between an 
electron and the nucleus of an atom is 


W = |(j,-A,-+M,-B,) dr (4.66) 


where j, and M, are operators representing the electron cur- 
rent and intrinsic magnetization and A, and B,, are the vector 
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potential and magnetic induction of the nuclear magnetic field. 
By replacing the fields by their multipole expansions we obtain 


W = > Mile) My a(). (4.67) 
LM 


The nuclear moments M,,,(n) are given by equations (4.62) 
and (4.63) while M,,,(e) is a sum of two parts, one being due 
to the electron orbital motion, the other to the spin interaction 
coming from the terms of equation (4.66) containing j, and M, 
respectively. The orbital part is (from equation (4.64)) 


4 
M,y(e) = sa Le 20, 4, dr. (4.68) 
and the spin part is (from equation (4.65)) 
“b 
MEE A) AiO, .M, ar |, (Ilene, 
(4.69) 


For L = 1 the last term in equation (4.69) is the familiar con- 
tact interaction between the electron spin and the nuclear spin 
when the electron is in an S-orbit. These are the same as the 
formulae for multipole interactions given by Schwartz [59]. 

Equation (4.66) represents the magnetic interaction between 
a classical electron and the nuclear magnetic field. The 
corresponding relativistic interaction is 


W = fea. A,dr, (4.70) 
where a is the Dirac [20] matrix giving the electron velocity. 
The multipcle expansions still hold if j, is replaced by ea 


and M, -> 0 because the spin interaction is already included 
in (4.70). 


CHAPTER V 
MATRIX ELEMENTS OF TENSOR OPERATORS 


WE are now in a position to calculate those parts of any 
physical problem which are concerned with angular integra- 
tions, the coupling of angular momentum vectors, or our 
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choice of coordinate system—in fact general rotational 
symmetries—which we might call the ‘geometrical’ parts. 
The basic technique is to exploit the simplest matrix reduction 
which we have already met in Chapter IV, namely the Wigner— 
Eckart theorem for the matrix elements of spherical tensors 


(JU\T golJ’M’) = (— 1K TMK M'Q)I||Txl|J"). (5.1) 


The aim is to derive some general formulae for matrix elements 
which depend only on the geometrical structure of the states 
and tensors involved, so then they can be applied to a variety 
of physical situations. Some of these applications are described 
in the next chapter. The results of such calculations are 
expressed in terms of the algebraic functions discussed in 
Chapter III. Owing to the extensive tabulation of these 
functions for a wide range of values of their arguments, this 
is a form very suitable for numerical work. 


5.1. Projection Theorem for Vector Operators 

First, we draw attention to the so-called Projection Theorem 
for vector operators V, which is the special case of the Wigner- 
Eckart theorem for tensors of rank one. It concerns matrix 
elements between states with the same J, 


JM'\V|JM) = JM'|I(S.V)|JM)/J(J +1) 
= 9,(V)JM'|J|J I). (5.2) 
This may be proved by expanding the right hand side, with 


the help of the Wigner—-Eckart theorem. The physical signi- 
ficance of (5.2) becomes obvious when we realize that 


J(J .V)/J(J +1) 
is the component of V along the unit vector 
I/V I (J +0}. 


According to the vector model, V is precessing about J, and 
its component perpendicular to J averages to zero (cf. Fig. 
5). The theorem states that the expectation value of.a vector 
operator in a state of sharp J is always proportional to 


V,§5.2 MATRIX ELEMENTS OF TENSOR OPERATORS 79 


the expectation value of J, as the second form of (5.2) emphasis. 
The coefficient g(V) is a generalization of the g-factor familiar 
in magnetic moment problems. In reduced form, (5.2) becomes 


SIV|[7> = 9V) VI (J +1}. (5.3) 


5.2. Matrix Elements of Tensor Products 
Next we consider a tensor 7’';9(k,k,.) which is itself a tensor 
product of tensors of rank k, and k,, as described in Chapter IV. 


z 


J-V 
Vid (J +1)} 


“S 
Fic. 5. The vector V precesses rapidly about J so that its component at right 
angles to J averages to zero. Then the component of V along O, is determined 
by the direction of J. 
We may write such a tensor 
T rglk,k.) = > By Sq—7hik2d@ —q|KQ). (5.4) 
qd 


The reduced matrix of the composit tensor T may be evaluated 
in terms of the matrices of the R and S by using (5.4) and 
introducing intermediate states between the two tensors, 


GT xll7’) = S (-—D*¥ IMT cold’M —Q) J MJ’ K M —QQ), 
Q 
= > (-1)?* {JM |J'KM —QQ)(KQ|k,k.qQ —9) x 
QqJ” 


x JMR, |J”M —q J" M —4|S;,,9-qJ'M —@)}, 
= ¥(—1)¥ hs (27 +.1)(2K +1)} WSS’ ky ky; KJ") x 
es 


x (IIR 97S, (5.5) 
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The geometrical significance of the Racah function is readily 
understood when we notice we have carried out a recoupling 
of three angular momentum vectors, expressing a scheme like 
[(k,k,)K,J’; J] in terms of one like [k,,(k,J’)J"; J]. The Racah 
W is just the transformation coefficient between these two 
schemes. 

When K = 0, (5.4) defines the scalar product of two 
tensors, although a different normalization is convenient 
(cf. equation (4.7)) 


R, Sp = (—DV(2k + Pool). (5.6) 


Equation (5.5) then reduces to 


Sy ee ay 
(PRs Sul) =D, (<0 (S55) VIRUS 
(5.7) 


IfS, = R,, and is Hermitian in the sense of (4.23), we may use 
(4.24) to reduce (5.7) further 


(J|[R, . R,||7) = 2 (—1)?|(J|| BIJ’) |. (5.8) 


We have said nothing about the nature of the tensors R, and 
S,; in the next section we shall consider the special case of a 
system composed of two parts, with R acting on one,S on the 
other. In such a case (5.5) can be expanded further. 


5.3. Reduction of Matrices for Composite Systems 


We are often faced with the problem of computing the value 
of an operator which refers to only one part of a composite 
system. Such a composite system may consist of two or more 
separate systems, or merely different aspects of the same 
system, such as the spin and orbital properties. Suppose we 
have a two-component system with angular momenta j, and 
Jo, and resultant J, while the tensor 7,(1) acts only on part 1 of 
the system. The coupled state vectors are given by (2.27), and 
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the reduced matrix of 7,(1) found by expanding them 
de [TQ )|\i527") = Y (—1)?* JM J"k.M —qq) x 
q 
XGJoI M|Tyq(1)| 91527’ M —q) 
= > (JM |J'kM —qq)(j,j.mM —m|J M) x 
qm 
X<(ijgm —qM —m|J'M —q) x 
x (jm jikm —qqK All TQ) ID 6(j252), 
= (1 P27! +1)(25, +1) WG GITs bed) X 
x AllT |i Oied2). (5.9) 


The Racah function appears for precisely the same reason as 
it does in (5.5). The mirror formula to (5.9), when T,(2) acts only 
on part 2 of the system, is easily obtained by remembering 
that changing the order of coupling of parts 1 and 2 only intro- 
Suces & PNAS, gM) = (14 jaf 

from the symmetry of the Clebsch-Gordan.coefficient. 

A simple example is that of a particle with spin, and an 
operator which acts only on the spatial coordinates. Such an 
operator could be the spherical harmonic C,,,. If the particle’s 
total angular momentum j is made up of orbital] J and spin s, 


(Isj||Cz||U'sf’) = (—1P P(t +1027 +P x 
x W(Il'jj’s Ls) U|Cz|U) (5.10) 


where (J||C,||J’) has been given in (4.17.) When/+1’+Liseven 
and s = } this simplifies considerably when the relation for W 
in Appendix II is used, 

aagliealras) =< (yey ah 4 a). oan 
which is independent of / and I’. 

More generally we may have a product of tensors with 
R(1) acting on part 1 and S(2) on part 2 of our composite 
system. The matrix of such an operator may be expressed in 
terms of the matrices of the component systems in the same 
way as (5.9), or by noticing that in this case each of the matrix 
elements on the right hand of (5.5) are of the form (5.9). Either 
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way we have to change the coupling of the four vectors j;. 
from a scheme [(j, J), (J1$}/’s K] to one like [(j,);),(dods)he3 
K] which indicates that an X-coefficient is involved. We soon. 
find 


idol [IT (erka)|[5t957") = (27 +1)(2K +.1)(2j, +1) x 
Sh ae 
x (2j2+1)} x fer | (jilfBu lt) GellSell92 >. (5-12) 
Jods he 
Putting k, = 0 of course, (5.12) reduces to (5.9). When A = 0, 
we have the scalar product (5.6), and.(5.12) becomes 


Goel Ria S20 ie eye a 
x 8S, 7’) —1 WHI Jod03 FAV ||Ral li) GellS,1|32)- 
(5.13) 


(This may also be obtained from (5.7) by using the Racah sum 
rule (3.21).) 

There is a simple geometric interpretation of the diagonal 
elements of (5.13) with 3, = j,,j, = 9,. The reduced elements 
on the right hand are the average or expectation values of R 
along j, andS along j,, respectively (see 5.2), while in the limit 
of large j,, 3. and J [Biedenharn 5}, 


(17-72 (25, +1)(27.+1)]' Wj jojo, KJ) & P,(cos 0). 
Here @ is the angle between the vectors j, and jg. 
cos 8 = {J(J +1) 9,9, +1) dele + 1) he. 


For k = 1, this is the familiar form for the scalar product of 
two vectors. 

There are many such scalar product operators, one example 
being the scalar product of the renormalized spherical har- 
monies of (2.9). 


Gl, L\|C,(1) . C,(2)[UUL') 
= 6(L; L)(—1)***[(20, +1)(20, +1)(20; +1)(21; +1)} x 


; eugl 91, [ RR aL: 
 Weiblss kE)(* 7 i) (' ae. (5.14) 
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Another important example is the spin-orbit coupling 
operator 1. s, with k = 1. This has matrix elements 


(isj||1. s||l’s’j’) = 6(s, 8’) 6(2, 1’) 6(j, 9’)(—1)?'-* W (Iss, 1j) x 
x [U0 +-1)(21-+1) s(s +1)(28 +1)]#. 
The explicit form for the Racah function form Table 4 leads to 


(Isj||l. s|[tsj) = 4G +1) —-U+1)—s(s+1)]. (5.15) 


In this particular case, of course, this result could have been 
obtained far more easily by using the vector relation 
21.s = j?—I?—s?. (5.16) 

It should perhaps be emphasized that in the explicit 
examples (5.11), (5.14), and (5.15) above we are concerned only 
with angular and spin integration; any radial integration has 
not been included explicitly. 

Other more complex systems with more than two com- 
ponents may be dealt with by repeated application of the 
techniques described here, or by the methods described in the 
next section. 


5.4. Systems of Many Particles 
When dealing with systems of similar particles, we are 
interested mainly in the matrix elements of two simple types 
n 


of operators. The first is the one-body operator F = > f(t), 
i=1 


where f(+) acts only on the coordinates of the ith particle. The 
other is the two-body operator G = > g(ij), where 


i>j 
g(t9) = g(9t) 
acts only between particles 1 and j, and 1 > 7 means a sum 
over all pairs. 

It is not appropriate here to go into the details of the 
construction of multi-particle wave functions; the reader is 
referred to other standard works [14], [17], [26], [49]. We 
shall merely indicate briefly the technique for evaluating 
matrix elements of operators like F and G. 

The wave function for a number 7 of identical fermions 
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(bosons) must be written so that it is antisymmetric (sym- 
metric) under the exchange of the coordinates of any two 
particles. Thus it is not possible in general to write it as the 
simple product of an antisymmetric (symmetric) function for 
nm—1 of the- particles times a wave function for the nth. 
However, it must be possible to write it as a suitable sum of 
such products. Suppose the system had total angular momen- 
tum J and z-component M, we may write (summed over «,, 
J,, a and j) 


la(n), JM) = > la,(n—1)J,, aj; JM) (a,(n—1)J,, aj|}a(n)J ), 
(5.17) 
where |a,(%—1)J,, aj; JM) = > |a,(n—1)J, M —m)lajm), x 


x (WJ, jM—mm|JM), (5.18) 


and «, «,, @ are any additional quantum numbers needed to 
specify the states fully. The parent states |a,J, M,) are 
orthonormal and fully antisymmetric (symmetric) in the 
first n—1 particles, while the mth occupies the state |ajm),,. 
The individual terms of (5.17) are not fully antisymmetric 
(symmetric) in all particles, but their sum, weighted by the 
numerical coefficients of fractional parentage (cfp) 


(Xp Fp» aj] }oT ) 


pp? 

must be. The cfp describe how the state |«J J) may be built 
up from its possible parent states obtained by the removal 
of one particle. The second relation (5.18) merely describes 
the vector coupling of the two parts to the correct total 
angular momentum. The orthonormality of the various states 
in the expansion (5.17) ensures that the cfp obey 

2 (oF {otg Ig aja yJ, aj|}a'T") = d(a, a’) (J, J’). (5.19) 

Kp pf) 

The expansion 5.17 is of principal interest when the many- 
body wave function |«JI/) describes a single ‘independent 
particle’ configuration, that is, one in which each particle may 
be assigned a set of quantum numbers (aj) (an ‘orbit’) inde- 
pendently of the others, the whole suitably vector-coupled to 
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give the correct total angular momentum. Configurations of a 
number of ‘equivalent’ fermions, that is, all occupying the 
same orbit (aj), have been widely studied in the literature and 
many of the cfp tabulated [26]. 

Since the n similar particles are indistinguishable, the 
matrix element of the one-body operator F’ must be just 2 
times that of f(n). Suppose F were a tensor of rank k, F,. When 
we use the wave function (5.17), f,(m) only acts on the nth 
particle, so we have a sum of matrix elements of the form 
(5.9) discussed above. We find 


n)J||F,,||a’(m)J ) 
= nD (—1)fotk-J-F'{(2I' + 1)(27 +L} WGI; ben) x 
x (ajl|fxlla’y’ ) (a(n) J {|ap(n—1) Jp, a7) (ap(m —1)J p a'7" |S’ (nT 
(5.20) 


summed over Oy oy a’, j,j'. Because of the orthogonality of 
the parent wave ‘auvnearhlenne, only those parent states common 
to |a(n)JM) and |a'(n)J’M) can contribute to (5.20). A very 
interesting result appears when we have n equivalent particles, 
all with angular momentum Jj, so the sum over j, j’ reduces to 
one term and 


aj) [Fy |’ (9")J’) = n(ajl|f,(2)||a7) x 

x DY (lye tt 77 (2d +1) (27 +1) W (Gj; eT) x 

X (ad {]a,Fp, 29) (pI y, aj|}a’J’)). (5.21) 
When F is scalar, k = 0, this becomes 

a(j")J||Folla’(g2)J’) = n da, a’) 6, J’) (aj||folla7). (5-22) 


The form of (5.21) shows that for these equivalent particle 
configurations the ratio of the matrix elements of two such 
one-body operators, F, and R,, is equal to the ratio of the 
corresponding single particle matrix elements. 


PI Fille’) — @ill Alley) 
ae 7 ~ (ajliry|lagy’ 


(5.23) 


86 MATRIX ELEMENTS OF TENSOR OPERATORS V, § 5.4 


independently of «, «’, J and J’. In particular, with the vector 
operator V we may take R, to be J, obtaining the simple result 


J(J +1) 


p)FL[V | a’ a) ean leet | an (G41) 2) calves 


(5.24) 


When discussing the projection theorem for vector operators 
in section 5.1 we defined a generalized g-factor by (5.2), (5.3). 
The relation (5.24) shows that for one-body operators the g- 
factor for n equivalent particles is the same as for one part- 
icle, independent of J and «. Putting F, =J in equation 
(5.21) yields a sum rule for cfp. 

Matrix elements of the two-body operator G may be evalu- 
ated in the same way, except that now we need to apply the 
fractional parentage expansion twice in order to ‘peel off’ the 
nth and (n—1)th particles. The value of G is then the value of 
g(n, n—1), times the number of pairs, 4n(n —1). In practice the 
two steps are often combined, and cfp defined for the reduction 
of the n-particle wave function to products of a function for 
the first n —2 times one for the last pair, appropriately sym- 
metrized [26]. Again, only parent states of the n—2 particles 
common to both n particle states can contribute to matrix 
elements of G between the latter. With independent particle 
wave functions this means the configurations for the two 
states can differ at most in two sets of single particle quantum 
numbers if the matrix element is not to vanish. 


5.5. Isotopic Spin Formalism 

In nuclear problems, neutrons and protons are often treated 
on the same footing by introducing isotopic (or ‘isobaric’) 
spin. 

The 1-spin operator + is identical with the Pauli spin 
operator o, except that it is said to operate in an abstract 
charge, or t-spin, space. The eigenvalues +1 of +, are 
then used to denote neutron or proton respectively. Charge 
independence of nuclear forces then leads to conservation of 
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total i-spin, T = } t;, and nuclear state vectors are simple 
i 


products of ordinary space-spin functions with an 7-spin 
function, |«JM, «;7M,) = |«JM)\«;1M;). Since the opera- 
tors t and T behave in charge space just as angular momenta 
in ordinary space, all the vector coupling and other techniques 
described in this book apply equally to them and their 
eigenfunctions. In particular, we may construct operators in 
this space (which, just as in Pauli spin space, can always be 
expressed in terms of t and the unit operator 1), whose matrix 
elements are to be evaluated as described in this chapter. If 
the interactions are charge-independent the operators repre- 
senting them must be scalar in the product charge space for 
the whole system, but if not (such as for electromagnetic 
interactions) higher rank tensors will appear. For example, 
the two-body charge exchange operator P’ can be written in 
the scalar form ey ee ey 

while the electrostatic interaction of two nucleons take the 
form 


(1 —a,)() Tada 


= [l—1,—7e,t F(t, T2) + V/ $7 29(71, T2)] e*/4145, 


in which scalar, vector, and second rank tensors appear. 


e” 


CHAPTER VI 
APPLICATIONS TO PHYSICAL SYSTEMS 


THE purpose of the present chapter is to apply the results of 
the previous two chapters to some physical systems of interest 
in atomic and nuclear physics, as examples of the usefulness of 
these techniques. 


6.1. Electromagnetic Radiative Transitions 
In this section we shall use the multipole expansion 
of a vector field, described in section 4.10.2, to define the 


7 
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electromagnetic interaction tensors and to calculate their mat- 
rix elements for the emission of a photon. We confine ourselves 
to non-relativistic particles, which is usually sufficient for 
nuclear problems. The treatment of relativistic motion follows 
very similar lines [34], [53], [54]. 

A particle of spin s and momentum p moving in an electro- 
magnetic field whose vector potential is A experiences (in the 
Lorentz gauge) an interaction [59], 


Hr) = = to fA(e) -p-+p.A(r)]-+9,8s. V A A(r)}. 
(6.1) 


g, is the spin g-factor so that the intrinsic magnetic moment of 
the particle is uw = gsS(eh/2mc), while gz is its orbital g-factor 
and eg, its charge. Terms quadratic in A are neglected in (6.1) 
since they lead to two-photon emission. Also, the first two 
terms of (6.1) are identical since A commutes with p in this 
gauge (div A = 0). 


6.1.1. The Multipole Operators 


The multipole expansion of the vector potential A, for a 
circularly polarized plane wave has already been described in 
section 4.10.2. It remains to insert these multipoles into (6.1) in 
order to find the corresponding interaction tensors. From 
(4.48) we have 


: ] 
A,(k, r) = e,e*" = aa (qAT u(r) ae eS) ne) 
IM 
(6.2) 
where V{L(L+1)AR = (2L41)97(kr) LC pa, 
and 
Aba = (HAV AAT: (6.3) 


corresponding to a flux of (k/47h) photons per second along k. 
F is the rotation taking the quantization axes into those with 
z-axis along k, and g = +1 for left/right circular polarization 
about k. 
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The charge dependent terms of (6.1) are proportional to 
A.V. The magnetic potential of (6.3) can be used as it stands, 
but it is more convenient to use the vector relation which 
follows from L.V = 0, 


Lf(r). V = —Vf(r).L, (6.4) 


because the operator L commutes with the radial part of the 
particle wave function while V does not. This is particularly 
advantageous when the long wavelength approximation is 
used, 

jrtkr) = (kr)*(2L41)', if kr<1 (6.5) 


Then f(r) is proportional to the solid harmonic r*C,,,, and the 
gradient formula of Appendix VI can be used. 

The long wavelength approximation also simplifies the 
electric radiation field in (6.3), for we have the relation 


VAL(rC yy) = (L+1)V(r7C zy) (6.6) 


and the electric part of the vector potential is the gradient of a 
scalar potential. We then use the anti-hermitian property of V 
and the Schrédinger equation for the states |x) of the radiating 
ticl 

oe ae [ —(H?/2m)V2+V—E,]]a) = 0, (6.7) 
to simplify the electric radiation matrix element. With (6.6) 
the electric component of the vector potential (6.3) is propor- 
tional to Vf(r) where f(r) is the solid harmonic r’C,,,. The 
matrix element of this part of (6.1) is then proportional to 


(B|Vf. V-+V. Vila) 

(B|Vf . Via) —(a|Vf. V|B)* 

— (B|fV?|a) + (al fW?|B)* 

(H, —B,)(2m/h2) (6 fx). (6.8) 
We have assumed the potential V in (6.7) is hermitian. Thus 
the effective electric multipole transition operator for long 
wavelengths is the solid harmonic r”C,,,, which is the same 


as the operator for the electrostatic moments, @,,,,, derived in 
section 4.10.1. This result is easily shown to hold for systems of 


I 


i 


ll 
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more than one particle (the energies L,, H, are then the total 
energies of the system). A derivation parallel to (6.8) uses the 
continuity equation [9]. 

For the spin terms of the interaction (6.1) we require the 
magnetic field vector H = V AA. From (6.3) we see immedi- 


ately that H",, = kAS y. (6.9) 


An analogous result follows for H* from (6.3) if we remember 
the wave equation for A, VA(V AA) = FA, 


Hoa (6.10) 


Collecting these results together we can write down the 
multipole tensor expansion of the interaction (6.1) when A, 
represents a plane wave circularly polarized about its direction 
of motion. 

clips 
avk L+1 
———_ _/____9"" ( R) x 
» (22 A ee 


x {(Qra + Oru) —1(M py +Myy)}- (6.11) 


Q, and Q;, are used to denote, respectively, the charge and 
spin contributions of parity (—)”, leading to electric radiation, 
while similarly M, and M‘, have parity (—)”*’ and lead to 
magnetic radiation. Explicitly, in the long wavelength 
approximation (kr < 1), and with the magneton 6 = eh/2mce, 


Ay(k,r) = — 


Qim(t) = egz (7 Cz), (6.12a) 
Myy(t) = 28g,V(r°C,y,). L/(L+1), (6.12b) 
Qta(t) = —kBgL(r7C,,,).S/(L+1), (6.12c) 
Mia lt) = B9.V(r°C yy). S. (6.124) 


These spherical tensors are of quite general application;{ the 
coefficients outside the brackets in (6.11) merely select the 


t The multipole tensors (6.12) agree with those used by some authors [11], 
[50], [59]. They are [47/(2L-4-1)]! times the tensors defined by others [9], 
[13], [63], because of the use in (6.12) of the more convenient Cz, instead of 
Yz. In addition to this, the Blatt and Weisskopf [9] multipole tensors are 
the hermitian adjoints of ours. 
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combination appropriate to the emission or absorption of a 
photon with definite linear momentum k and circular polari- 
zation q about k. The interaction for other polarization states is 
given by a simple combination of terms like (6.11). For ex- 
ample, from (4.37) we see the operator for photons with linear 
polarization in the x-direction is just 


1 
Hk, r) = vgltt alk, r) —H,(k, r)]. 
Of course, if the polarization state is of no interest, we use an 
incoherent sum of the g = +1 terms. 
The tensors (6.12) are hermitian in the sense of section (4.8). 
and have simple time reversal properties (section 4.9). If T+ 
stands for any tensor of (6.12),t 


Dg = —67T;,,97 = (a ae (6.13) 


where p = 0 for the magnetic and p =1 for the electric 
terms. 

For systems of more than one particle the interaction 
operator will be a sum of terms (6.11), one for each particle. 


6.1.2. The Matrix Elements and Transition Probabilities 

The selection rules operating for matrix elements of the 
tensors (6.12) are given immediately by section 4.7.2. In 
particular the different parity of electric and magnetic multi- 
poles of the same rank JL ensures that only one can contribute 
for each allowed ZL, to transitions between states of definite 
parity. For the same reason, when two multipoles of rank LZ 
and L+1 contribute, one will be electric and the other mag- 
netic in order to conserve parity. The time reversal properties 
(6.13) ensure that reduced matrix elements of H(k, r) will be 
real (section 4.9.). 

The static electric and magnetic multipole moments of a 
system with total angular momentum J are given by the 


{ Apparently (6.13) is not obeyed by the Qzy of (6.12a), but it should be re- 
membered that this is only an effective operator in the sense of equation (6.8). 
When the hermitian conjugate of a matrix element (A|Z'|x) is taken, the 
factor (HZ, — Eg) in (6.8) also changes sign. 
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diagonal, M = J, elements of the tensors (6.12) (summed over 
all particles in the system). The magnetic dipole moment is 


w= BUS 9g, Lb +9 8I) = VI WIFY} S|[4+ My) 
(6. 


14) 
while the conventional electric quadrupole moment is 
Q — 2(Qo)/é; 
J(2J —1) ) 
= 2__y2 = 2{|—___——__} J Hb 
eQ = e(JJ |g 1 (32? —1?)| JJ) le aee (J||Qel|7) 
(6.15) 


The other static moments of interest are the magnetic octu- 

pole moment Q [41] [59], and electric 2*-pole, or hexadecapole, 

moment Q [74] (this reference uses the symbol M,, = 8Q‘) 
— —Bp JT |(4giL +9sS).V(r3C30)|JJ ) 

( i ieee) ) J||Mg+Ms||7) (6.16 

 ) eastgaeaersey) Ie ft 

Qa) = i 

a 


= ae GIO. (urea 
= rare nessa ar45y) VICI» (6.17) 


(JF |g (3524 —302?r? +3r4)| JJ) 


To first order, the probability amplitude for a transition 
from a state |o,J,J1,) to |a.J,M,) with emission of a circular 
polarized photon along the direction & is proportional to the 
matrix element of the interaction (6.11), 

q k \i pl 
Abt a(K) = — (<3) SS GaP x14(R) SiMi| Tia 


LMn 


(6.18) 


where z is e for electric (q” = 1), 2 is m for magnetic (qg” = q), 


(iy 1 La) ee ee 
(2L—1)!! (St ee se (6.19) 


a= 


and 7',,, stands for the tensors (6.12). If we do not observe the 
orientation of the spin J, of the final system we must sum over 
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M,, and the probability P4, (k) of emission from the initial 
substate with z—component M 1 is equal Ho 2 > lASear (k)]*. If 


linear, not circular, polarization is observed, we Retake a coherent 
superposition of g = +1 terms A? as described in the previous 
section, while a polarization insensitive measurement requires 
an incoherent sum over q. If the radiating system is in a 
cylindrically symmetric environment, so that M, is a constant 
of the motion if the z-axis is chosen along the symmetry 
axis, the total radiative probability for photons along k is 
obtained by weighting each P{, with the population or 
relative probability w(M/,) of the substate , 
Pk) = ¥ w( I) P5,(k) = ¥ ~My) |44, ae (K)? 
M My M, 

(6.20) 
(If the system does not possess cylindrical symmetry, J,, 
cannot be made diagonal, and w(//,) has to be replaced by a 
more general density matrix p(M,, M,): (see section (6.4).) 
Using the Wigner—Eckart theorem of section 4.7.1 and com- 
bining the rotation matrixes according to (2.32), we soon reduce 
the aan distribution (6.20) to 


Pah) = = iD, Be Bx(J1)Px(cos B) Wii LL’; KJ2) x 
x (204+ 1)} (—1) etd eti-L'-Kang’ (LL'q—q|K0) x 
x Slee Talo) Silat Tele )*. (6.21) 


B is the angle between k and the symmetry axis, and B, 
describes the orientation of the initial system (for example, an 
assembly of oriented nuclei [12], 


B,(J) = > w(M)\(-1) -™4/ (20 +1) (KOJI M —M), 
M 
BJ) =1. (6.22) 


If the system was initially randomly oriented, so w(M,) is just 
the statistical weight (2/, +-1)~1, only the K zero term survives, 
and (6.21) becomes 


k & 
Puy) = 55 >| VilleiTie)P/@L+1) (6.23) 
Ln 
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which is independent of the polarization g and the direction of 
emission k, as expected. The total decay rate is obtained by in- 
tegrating (6.21) over the direction of emission k and summing 
over the polarization q. We note that, while different multipoles 
interfere coherently in the angular distribution (6.21), there are 
no interference contributions to the total decay probability 
(6.23). The decay probability per second is 


al EAU E41) ||| TE| oI? 
7 AS T(2L—1)iN}? (2L41) 


while the ‘reduced’ transition probabilities introduced by 
Bohr and Mottelson are just 


L+1 
B,_..(7L) =| ||T7||7 2)? (= ) (6.24) 


Comparison with (6.21) also shows that the matrix element 
which appears in formulae for the angular correlation of 
gamma rays is, in the conventional notation of Biedenharn 
and Rose [8], [19], 


(GLI) = N Get Tal9, /2L+1) 


where the normalization factor N is independent of L and 7. 


6.1.3. Single-particle Matrix Elements 

We can now calculate the matrix elements of the multipole 
tensors (6.12) between the states of a single particle in a central 
field. For a spin—4 particle with spin-orbit coupling these 
states are 


4jm) = ae 2 |4o)u,(r) Ges (Op) dm — —oo|jm) 


where 2’ is included to give the required time reversal prop- 
erties (section 4.9.). The matrix elements now have the forms 
discussed in Chapter V. The elements of Q,,, follow immedi- 
ately from (5.11), while M,,, has the form of (5.5). The spin 
termsQ,,, and M{,, both have the form (5.12). Using these and 
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some of the standard reduced matrix elements listed in 
Appendix VI we obtain 
(451||Qr||lo470) = eg, J (L) dy, 
L451||Miz||lo390) re 28g p08 ( —1)3 (LI Wh ole; 3L) x 


Le FI + DANAE AYL/E-+1 (2 ‘ | 


I 
ht jl|Qr||lo4J2) = 49PkI(L)(a,—a2)byo/(L +1) 
(1,494||M{||l2392) = 39,8 J(L—1)(L—a, —ag)by, 

yi be jn } j i L 
where by, = 1-2 (—1) + ait 7 i): 
a = (I-j)(2j+1), 
and J(L) is the radial integral 
J(L) = fey? un(r) dr. 
To conserve parity (L+/,+/,) must be even for electric 
operators and odd for magnetic. 
For M1, El, and £2 transitions the elements are very 
simple. For M1 we have l, # I, transitions are forbidden, and 
(14 j|[My + Myl[39) = (9 +) GF 0d, 
(44M +My +1) = BG.—91)V (I+ D/(21-+1). 
(jl) is the Schmidt value for the magnetic moment of a 
particle in an orbit (4, 7) [11]. 
(Gl) = jBlgc(ge—gz)/(21+1)} as jf = 143. 


The dominant charge contribution to #1 transitions has 


poms (45|| Qulllod7) = degrF (2 -/GG +1), 
G4 QlZ+1, 4941) = —diegrJ(1)[(27 +3)/G+1)}F. 
Similarly the charge contributions to H2 elements are 
ee) 
ae | 
-1, (_ 3(2+3)_\ 
= ter 4 (| 
C35] Qallied5 +1) = dee T(2)~% | Te 
; oe el 
(4g O42, 4.542) = dea) (oes 


(4jl] QuI|133) = —teg,J(2)( 
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6.1.4. Systems of More Than One Particle 

Systems of more than one particle may be treated as 
described in section 5.4, since the electromagnetic interaction 
is a one-body operator. To emphasize the possible importance 
of the coupling to other particles we consider here just two 
particles, with a transition from a state | jjI,M,) to | WB). 
The matrix element is given directly by (5.9). In particular let 
us compare the 2”-pole transitions (A) j,j > j2and (B) j,j > 9}, 
when J, = J,+L. Expanding the Racah functions in (5.9) the 
ratio of the square of the reduced matrix for transition A to 


that for B is (2) —Jy)! (Jot 2f, +D)! 
(2, — Je)! (Fg +2541)! 

The interesting cases arise when j, and J, are large, and 9, J, 
are small. For example take 7, = 33, 9 =3, J, =7, and 
J, = 2, so that L = 5. Then, although the same single 
particle matrix element (|| 7 a8 is involved in both transi- 
tions, the ratio of their probabilities is 1/728, while transition 
A is reduced by the angular momentum coupling to 1/2070 
the intensity of the single particle L = 5 transition ? > 43. 

More detailed considerations of transitions betitean many - 
particle states have been given elsewhere for both 7—j and 
L—S coupling [14], [43], [46], [49]. Some interesting results 
arise from the general properties of matrix elements discussed 
in Chapter V. For example, the relative probabilities (or ‘line 
strengths’ [17]) for electric transitions between the com- 
ponents of two Russell-Saunders multiplets follow immedi- 
ately from the decoupling relation (5.9). The reduced transition 
probability B(EL) of (6.24) between the L—S states |«,L,8,J,) 
and |a,L,S,J,) becomes 
47 By (EL) = dg 5 (2L,+1)(2J,+1)(2L+1) x 

x WJ, JL, L,; LS,)| (a, L,||T4\o,L,)|. (6.25) 

So between states of the same two multiplets, but with 
different J, and J,, we have the purely geometrical ratio 


By 5 (EL) = (23, +1) WJ J, 0,L,; LS;) 
B, 5 (£L) (27, +1)W*(J,J,L,L,; LS,) 
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Again, if we sum (6.25) over final angular momenta J,, we get 
ae By .7 (EL) = (2D+1)| (oq, Ly||T4]|.L,)|?. 


So, between two L—S multiplets the total intensity from an 
initial Zeeman level with J,, J, is independent of both J, and 
M,. Hence the sum of intensities from the Zeeman compo- 
nents of a given J, is proportional to its statistical weight 
(2J,-+1). Conversely, the sum of intensities feeding the Zeeman 
levels of a given final J, is proportional to (2J,+1). 


6.2. Interaction between two Systems 

In this and the following section we shall discuss the inter- 
action of two systems, each of which may itself possess some 
structure. If the combined system is isolated, the Hamiltonian 
is invariant under rotations and any interaction terms coup- 
ling parts 1 and 2 must be expressible as scalar products of 
tensors as in (4.4), 


VQ, 2) = p3 Rx(1).Sx(2). (6.26) 


The matrix elements of such products are given by (5.13) in 
terms of the matrix elements for the component parts 1 and 2; 
it merely remains to discuss some specific forms for the 
tensors Ry and Sz.t We shall also mention briefly a case, 
anisotropic hyperfine structure, where the effective interaction 
has to be represented by a tensor product rather than the 
scalar form (6.26). 


6.2.1. Interaction of Nuclei with Atomic Fields 
The interaction between two charge distributions was 
discussed in section 4.10.5, and the form (6.26) found for both 


t With a second quantization treatment of radiation fields, as mentioned in 
section 4.10.4, radiative transitions are also induced by scalar product inter- 
action terms like (6.26). For example, the electromagnetic interaction will 
have the same form (6.11) as in the previous section, except that, following 
(4.50), each multipole tensor will be associated with a photon creation or 
annihilation operator in the proper scalar combination. Matrix elements are 
then taken between product states for field and radiating system. 
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the electric interaction (4.59) and the magnetic interaction 
(4.66). Important dynamic applications of these interactions 
occur in the theory of internal conversion [53], [55], whereby a 
nucleus transfers excitation energy to an atomic electron, and 
Coulomb excitation of a nucleus by bombardment with 
charged particles [1]. Here we shall consider briefly the station- 
ery effect of the interaction between a nucleus and its elec- 
tronic environment which leads to the hyperfine splitting of 
atomic spectral lines [41] [59]. Ignoring penetration into the 
nucleus by the electrons (sor, > r,), the interaction for iso- 
lated atoms is given by (4.59) and (4.67) 


Ve, n) = ZC ral) zac) + MG acle)M rarlm)]. (6.27) 


The nuclear electric moments Q(n) are just those defined by 
(6.12a) in the previous section, while the electric moments 
M(n) are the sum of the charge and spin contributions (6.12b) 
and (6.12d). (Summed, of course, over all nucleons.) The 
corresponding moments Q(e) and M(e) for the electron 
cloud are given by (4.60) and (4.68), and are similar to the 
nuclear moments except in radial dependence. The quad- 
rupole operator Q,,(e), for example, is just 40?V/dz? at the 
origin, where V is the electrostatic potential due to the 
electrons. 

The small contributions from r, < 7, may be deduced from 
section 4.10.5, for example the higher terms in the expansion 
(4.58). 

The first-order hyperfine energy shift W,,, for an atomic 
level with total angular momentum F (vector sum of nuclear J 
and electronic J) is just the expectation value of V(e, ~)+. It 

{ Second order effects of low multipoles may be important in the interpre- 
tation of higher multipoles; for example, the second order term for the dipole 
I = 1 has a part which looks like a first order quadrupole, L = 2, term [41], 
[59]. Also the much closer approach of the meson in y-mesic atoms allows 
higher order effects to be important which involve excitation of low-lying 


nuclear levels (nuclear ‘polarization’) through the off-diagonal elements of 
V(e, n) [79]. 


| 
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follows immediately from (5.13) that , 
Wisp = UIF|V(e, n)|LIP) 
= 2 tee wee 27) WIL. FL) x 


x U|TH(m)|| [THe (6.28) 
where T* = Q, T™ = M, and parity conservation ensures 
that only Q, with Z even, and M, with L odd, can contribute. 
We can invert (6.28) to express the reduced matrix elements in 
terms of the observed W,,;,, using the orthogonality (3.17) of 
the Racah coefficients 


[22 +127 +I) FAI PFIE) ATA) 
= ¥(-1) YY -F((2F+1)(2D41)] WLS, FL) Wyzp. 
F 


(6.29) 


The evaluation of the reduced matrix elements has been 
discussed in the previous section and elsewhere [11], [41], [59]. 

Allied problems, such as in molecular spectra, or with the 
application of external magnetic fields, can be treated by 
similar techniques [41], [50]. We shall confine ourselves to a 
few remarks on an ion which is not isolated, but situated in a 
crystalline electric field. Such an electrostatic field can be 
expanded in solid harmonics about the position of the ion 
(following section 4.10) 


a A, 1*C,(9, @). (6.29) 
kg 


The coefficients a,, are determined by the nature and sym- 
metry of the crystal lattice [10]. The matrix elements of (6.29) 
are readily evaluated using the Wigner—Eckart theorem and 
the techniques of Chapter V [25]. In the ‘single-electron’ 
approximation for the ionic states this is almost trivial. An 
alternative method for diagonal matrix elements between 
many-electron states is to find an ‘operator equivalent’ for 
each term in (6.29) [68]. Between states of the same J we use 
J,, J,, and J, to construct operators which have the same 
matrix elements as (6.29); for example 


PO oy > % (r*)[3I2 —J(J +1). 
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Then «,(r?) is proportional to the reduced matrix element of 
rC,, and techniques have been devised for their evaluation [68]. 
The other part of the operator equivalent gives matrix 
elements proportional to the Clebsch-Gordan coefficient of the 
Wigner-Eckart theorem; for the example above, 


(JM|3J7 —J(J +1)|J 1’) 
= (MM )J(J +1)(2J -1)(2J +3) JM J2M0). 


This approach, however, will not give off-diagonal matrix 
elements. 

This overall crystal field polarizes individual ions and 
induces anisotropic hyperfine structure. The environment is 
no longer spherically symmetric and the effective nucleus- 
electron interaction can include tensor products of higher 
rank than the scalars in (6.27). A particular example is the 
spin-Hamiltonian of Abragam and Pryce [10] which includes 
a term 


Ha, = AS. Ie = BSS IeeCs 
= AS ol see —C)(S, ee! 4) oe +C)(S,1_y asi gh 


S is an effective spin of the atomic electrons. The symmetry 
characteristics of Hg; are made evident if we express it in 
terms of product tensors (4.6) built from S and I. 


Hy, = oS .1+67,,(S, 1) +y7(7,,(S, 1) +7,_,(8, 1] 
where (6.30) 


«= HAFB4C), P=—Z2A—B-0), y= H(B-0). 

The effect of the anisotropy is to introduce second rank 
tensors into the interaction which will connect states of 
different F’; that is, F is no longer a good quantum number in 
the presence of this interaction. However, if axial symmetry 
is restored by putting B =C, the T,,, terms vanish and 
F, =I1,+8, does remain a constant of the motion. Hg, is, of 
course, diagonal in S and J even though it mixes different S,, I, 
values. 
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This interaction is the basis of: the Bleaney method of 
orienting atomic nuclei [12]. The crystal field orients the 
ionic electron clouds which in turn, through the anisotropic 
interaction Hg,, orient the nuclei. Of course, to be effective 
the sample must be cooled to a temperature 7 such that kT ~ 
the hyperfine splittings induced by Hg,;. Because Hg, is 
symmetricin +J,, only nuclear ‘alignment’ can be produced; 
that is, only even-order polarization moments are induced, 
and (I”) vanishes for n odd. 


6.3. Interactions between two Particles in a 
Central Field 
A particular case of the interaction of two systems which is 
sufficiently important in the theory of atomic and nuclear 
structure to warrant separate treatment, is the mutual- 
interaction of two particles moving in orbits in a central field. 
The particles may possess spin. Then the various interaction 
terms, although scalar in the overall product space of spins 
and orbits, may be classified according to their properties 
under spatial rotations in the spin and orbit spaces separately. 
If the position coordinates are denoted r, and r, and the spin 
operators S, and s,, the interaction can be written in the scalar 
product form 
[Ay = p3 R,(T,, Pe) - Sz(S,) So). (6.31) 


Rx(rir2) operates only on the position coordinates, Sx(si, se) 
only on the spins; each is built up from its arguments in the 
way discussed in Chapter IV. 

A nuclear system of neutrons and protons may be described 
by the isotopic spin formalism; the extension of (6.31) and the 
discussion below to this case is straightforward (see section 
5.4.2 and [26]) and will not be considered here. 


6.3.1. Spinless Particles and Central Forces 

Simplest is the scalar K =0 (or central) force .acting 
between particles without spins. This will be just a function 
of the distance r,. between the particles, J(r,,) say, where 
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Ty. = T,—T,. The matrix element of this may be evaluated by 
the Slater method [26], [69] by expressing J(r,.) itself in the 
form (6.26). We write 


De — > toe Cost) (GRY) 


where 2 
Tilt %) = }(2k +1) [J(ry)P,(oos w)sinw dw. (6.33) 

0 
Expressions for J, for various interactions have been given 
[69], and the particular case J = e*/r,, has been discussed in 
Chapter IV. w is the angle between r, and r,. The Legendre 
polynomials may be expanded by the addition theorem (2.25), 
P,(cos w) = C,(1). C,(2). If the particles are spinless and their 
orbital angular momenta are coupled to a resultant L, we may 
immediately carry out the angular integrations of the matrix 
element of (6.32) using (5.14). In the usual notation, and using 


single particle wave functions 
|tm) = UnplT) Fo: g) 
(L1,LM|J(r,.)\l1,L'M’) 
= (LL (MM) Sf llylelily; L)R (Ulf I) (6.34) 
k 


where f, is just (1,1,L|C,(1). C,,(2)|0,15.L) given in equation (5.14). 


The f, have the following symmetrics: 


ALLL; L) = f(GLL 1; L) = f,(Io 1,1, ; L) 
and (6.35) 
fo = 6(1,0;) d(Ia1). 


The condition that both (k+1,+1,) and (k+1,+15) be even 
ensures that J(r,,.) does not connect states of different total 
parity, and k is limited to <1,+1, and <l,+1}. The R* are 
Slater radial integrals, parameters independent of the total 
angular momentum value L, 

BL =F dr,{r3 Ory TM yy Ty) My, (Ty )Uy, (71) Uy, (72) My, (72). 


(6.36) 
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u,,(r) is the single particle radial wave function; the indices n 
denoting the principal quantum numbers of the orbits have 
been suppressed for clarity. In the limit of very long range 
(compared to the extension of the wave functions w/(r)) 
J(r,.) becomes constant, independent of r, and r,, and (6.36) 
ensures that J, vanishes unless k = 0. In this limit all matrix 
elements vanish which are off-diagonal in the single particle 
orbits, and the others are independent of L: 


Qylo LMF (1,9) (Ul, LM) 2 WUE) CMU EOI es Gain) 
The other limit of zero-range is of interest for nuclear structure. 


The expansion (6.32) for a delta function is 


d(r,—Yr,) = d(7, —7,) d(cos w —1) 


Qari; 
: 2k Cc Cc 
== arr) D, (2k+1)C,(1).C,(2), (6.38) 
then immediately in this limit 
+1 
ee AT 6(r,—7,) = (2k+1)J, (71, 72), (6.39) 
mrt 


if J(r,,) = A 6(r,—r.). Thus we have R® = (2k+1)R© and 
this allows us to carry out the sum over k in (6.34) using (3.11). 
The zero range matrix element becomes 


LEE OB FIR, — TRAP OWE (6.40) 
where 


P(UlgtI,; L) 
7 ; ; a ie ; 4] , ’ + 
i i: 0 ol ho o of lea + IA + 1H (2h +1)(205 +1] 


which, of course, vanishes unless (J, +1,+Z) and (1, +1,+L) 
are both even. The radial overlap RF) has become simply 


RO = = fr dr w, (7), (r) ty (7 uy, (7). (6.41) 
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Occasionally it is more convenient to work with the uncoup- 
led matrix elements, (1,m,1,m,|J|l;m,1,m,). This will be so if the 
environment is not spherically symmetric, such as for the 
motion of individual nucleons in a strongly deformed nucleus 
[13], so the total L is not a constant of the motion any more 
than the individual J. The same multipole expansion (6.34) is 
used, and the angle integrations carried out by use of (4.16) 


(l,m,l,me|J(r12)|L, mlm.) = 0(m, +45, m, +My) X 
x > ck(l,m,l,m,) c(l_mlim,) RO; (6.42) 
k ea 


c(i’) = [20-412 +1) —1)"(F : Alle a 


m—m' m'—m 
= (—1)"—™ c*(l'm'Im). 


Values of c¥ have been tabulated [17], [64]. 

Finally, it should be remarked that a technique different 
from the application of (6.32) has been developed by Talmi 
[32], [70], [73] for use when the functions w,,(r) are eigen- 
functions of an harmonic oscillator potential well. 


6.3.2. Particles with Spin; Central Forces 

When the particles have spin, K ~ 0 terms may appear in 
(6.31), and the scalar K = 0 matrix element is itself modified. 
Because of the form of (6.31) it is most convenient to use an 
L-S or Russel-Saunders coupling scheme, when we get an 
example of the general matrix element (5.13){. 


(LSJIM|V(12)|L'S'J'M') = 6(JJ') 6(MM’)(—-1)"-*-* x 
Se W(LL'SS'; KJ)4/(2L+1)(L||Rx(1, r2)|[L’) x 
K / (28 +1)(SI]Sx(s,, 8,)||S’). (6.43) 
If the scalar part (K = 0) of V(1, 2) is spin independent 


{ Matrix elements in the j7—j coupling representation are easily obtained 
from (6.43) by using the transformation (3.23), (3.24) 


GrjeI|V|jv5e'I) = YS CLSI|V|L‘S’I) x 
LL’ SS’ 


X C(8i) Ir» (2a8a)J a3 T\(Als)L, (318 a)S3 J)<(1’83)91’, (ba’8a)5a's T\(by0q')L, (318,)S’3 J). 
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(So(s,8.) = 1, (SI|1||S) = 1), the decoupling factor represented 
by the Racah function gives unity. In addition, however, 
there may be a scalar spin interaction, 

So = 8, Se; 

(Sl|si. Sel[S) = 3[S(S+1)—si(si +1) —so(so+1)]. (6.44) 
When s, = s, = s we may introduce also the spin exchange 
operator Ps, Because of the symmetry induced by the Clebsch- 
Gordan coefficient in the coupled state |ssSJ2) when the 
angular momenta are equal (see section 2.7.3.), P*® has the 
eigenvalues (—)°~**, 

PlssS Bf) = (=1)° “|sa8 11) (6.45) 
hence 
(S||P']S) = (—1)8~**. 
In particular, when s = 4 we can use (6.45) to give a represen- 
tation for P*: 
Pt =442s,.s, = }(1+0,.4,). (6.46) 
The second form in terms of the Pauli spin matrices follows 
since S = }o here. This allows us to evaluate the expectation 
value of the spin exchange operator for  spin-} particles, 


Ah EDP iy ( (known as Hund’s operator) without recourse 


to the snetnodls of Section 5.3. Because of the meaning of P* 
n (6.45), the value of M is the difference between the numbers 
o symmetric and antisymmetric pairs of spins. Using 6.46 
and remembering that 
Se (sa se 
a t#j 
we get 


M = [}n(n—1) +28*—25 $i] 
and since s = 3, for a state of total spin 8 


(S|M|S) = 4[}n(n —4)+28(S+1)). (6.47) 


6.3.3. Vector and Tensor Forces 
The terms with K = 1 and 2 in (6.31) are referred to as 
vector and tensor forces, respectively. The only vector term 
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linear in the spins and momenta which satisfies the require- 
ments of symmetry and space- and time-inversion invariance 
is (Wigner [76]) 


J(r12)L(1, 2) .S, (6.48) 
where 
So EY Sas, 
and 
L(1, 2) =3(r,—P_) A (P, —Pe) (6.49) 


is the relative orbital angular momentum. The matrix elements 
of S,(s,s,) =S are simply dealt with, but R(r,r,) = I(r.) 
L(1, 2) requires more attention. The terms in the expansion 
(6.32) of J(r,;.) have to be coupled to those from the corre- 
sponding expansion of L(1, 2), 


L(1, 2) =4{L(1) +L(2)—1r, AP2+P, AT2}. (6.50) 


The first two terms, the one-particle orbital operators, are 
straightforward. To manipulate the other two terms we use 
the relations p = —iV and 


a 
Wee (6.51) 


v—¢ 
where (C, is the unit vector along r; L(1, 2) then becomes 


re] re] 
L(1, 2) =1{La) +12) +4G,(1) AG) (n= 


ar, 2 ar, 
+€,(1) A(Cx(2) AL(2)) 2 
2 
+0,2) A(G,(1) AL) 2}. (6.52) 
1 


When combined with the expansion (6.31) of J(r,,), this is a 
series of tensor products. These may be evaluated by straight- 
forward, although somewhat lengthy, applications of the 
techniques described here and in the previous chapter [23], 
[37]. 

The tensor force with K = 2 is simpler because of the 
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absence of the gradient operator. It is usually written 


J (1 15)S(12) 
where 
S(12) = (8,. e)(S.. @)—4(S, - SQ) (6.53) 
= R,(e). S,(s,8,), 


pe being the unit vector along rj,. The interaction of two 
magnetic dipoles p, and p,, for example, would have 
I(r) = 919o/Tto, Where g, is the g-factor for the ith particle, 
wt; = 9,8,;. From (5.12) the matrix elements of the spin tensor 
for two spin } particles aret 


(S||S_(818,)||S’) = 6(S8’)6(81)+/ &. (6.54) 


It is diagonal and vanishes for singlet, S = 0, states. The 
coordinate tensor is just a spherical harmonic and is readily 
re-written in terms of tensors acting on the coordinates r, and 
r, separately, 


rieC,(e) = r2C,(1) +13C,(2) —+/6R,(r,1,). (6.55) 


It is convenient to expand, not J(r,,), but 


] 
in = > Hem Cetllo Ge (6.56) 
12 k 


Again we have a series of tensor products to be evaluated by 
repeated application of the standard techniques. Because of 
the 1/r?, in (6.56), the separate radial integrals contain diver- 
gent parts, but these cancel in the final result [23], [36]. 

Again, when the radial functions for the single particle 
motion are eigen functions of an harmonic oscillator well, 
the technique of Talmi [71] greatly simplifies the matrix 
elements for vector and tensor forces. 


6.4. Multipole Expansion of the Density Matrix 
It is beyond the scope of this book to give a detailed 
discussion of either the basic theory or the widespread 
t S(12) is defined in (6.53) with the actual spin operators s. Thus for spin } 


it is } of the tensor operator often defined using instead the Pauli spin oper- 
ators, o = 28. 
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applications of the density or statistical matrix and its associ- 
ated statistical tensors. We can only outline the techniques, 
briefly referring the reader to the literature for more details 


[8], [19], [30], [81], [72]. 


6.4.1. The Density Matrix 
If all the (identical) component systems of an assembly are 
described by the same wave function 


pao, wvfaheal (6.57) 


the assembly is said to be in a pure state, and the density 
matrix p for the assembly is defined by 


12 * 
Pam zie A, ans 


tre = 1. (6.58) 


It often happens that we do not have such complete infor- 
mation about the assembly. An example is our knowledge of 
the individual spin orientations in a partially polarized beam of 
particles or assembly of non-interacting atoms or nuclei. This 
more general mixed state can always be described as a weighted 
mixture of the pure states (6.57) and (6.58), so that p becomes 
an average over all N component systems of the assembly 


1 N 
Pram = N » Ay(¥)Q,,(¥)* = Any, : (6.59) 


Clearly ep is Hermitian, and from (6.57) and (6.59) the expecta- 
tion value in the assembly of some operator O is given by 


(O) = > Pam mn = tr(pO). (6.60) 


If the representation we have chosen makes p diagonal (for 
example if p describes a paramagnetic gas in a uniform 
magnetic field, and m is the component of spin along the field 
direction) we can write 


Pam = WM) Onn: (6.61) 


w(m) is then the population function, or probability of 
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finding one of the component systems:in the state |m). (6.60) 
then takes a particularly transparent form 


(O) = ¥ w(m) Onm: (6.62) 


m 


6.4.2. Statistical Tensors or State Multipoles 
When the angular symmetries of the assembly are of 
interest (that is, its properties under spatial rotations), the 
natural choice for the basis ED) in (6.57) are the eigenstates 
\xJ M) of the angular momentum J and its z-component M 
(« denoting any other quantum numbers required). From (6.58) 
and (6.59) we see p transforms like an operator under a change 
of representation p’ = RtpR if |n)’ = Rin). In particular, 
upon rotating the coordinate axes through Euler angles (af), 
an element of p referred to the new axes is expressed in terms 
of those referred to the old by 
Pramas a = & (PinlOBy))* rosnaswDurelby)- (8.63) 
otal 
It is more convenient however to choose linear combinations 
of the elements (6.63) that form an irreducible representation; 
that is, carry out a multipole expansion of the elements of p. 
This is clearly just a matter of vector addition since from (6.59) 
the elements of e are bilinear combinations of tensor com- 
ponents @,,,,. So we can write 


PadM,os’M’ — 


> prolad, oJ’ )(KQ|JJ’, —MM')(—1)*-"-™. (6.64) 
KQ 


The pgg were called statistical tensors [28], or state multipoles 
[29], by Fano, and by inverting (6.64) may be written 


pxg (ad, a'J’) = | 
ee er KO) (6.65) 
M 


$ Unfortunately a variety of notations and definitions are in use for the 
statistical tensors. Ours agrees with some others [8], [6], [28], [29]. Several 
authors define them so that the tensor of rank K behaves under rotations like 
.Y%, contragrediently to ours [81], [19]. This is appropriate if the statistical 
tensors are regarded as the coefficients in an expansion of the density operator 
into a set of multipole operators. 


110 APPLICATIONS TO PHYSICAL SYSTEMS VI, § 6.4 


where 
‘yy? ¢ —4 
Poo Od aI") = b,,(2F +1) 2 PaJM,a'JM. 


When «, J have unique values, the latter becomes 


Poo = (27 +1)“ire. 


From its construction pxg is a component of a spherical 
tensor, ; 
PEo = 2 ProZaq(*By) (6.66) 


which may be confirmed directly from (6.63) using (2.31). 
The Hermitian property of p reappears as 
Pro el cdl) (mes ore ee (6.67) 

The definition (6.65) shows we can define tensors of rank K 
where |J—J'| < K < J+J’. It is then consistent to talk of 
the assembly possessing dipolarization (K = 1), quadripolari- 
zation (K = 2), ... 2X-polarization, if the corresponding 
tensors do not vanish. Only for spin-} systems does the term 
‘polarization’ have a unique meaning.t The monopole tensor 
(K = 0) is merely a normalization constant. For systems of 
sharp J the other low rank tensors have simple interpretations 
[6]. For example the dipole tensor is just the expectation value 
of the angular momentum operator 


PiJS) = J.) +2414, 
and the Q = 0 component of the quadrupole tensor is 
Po JS) = (8I7—I(J +1)) RI (J +1)(27 —1)(27 +1)(27 +3)]*. 


If it is possible by a rotation of axes (6.66) to make p diagonal 
in M, then with this choice of axes only tensors pyg with 
@ = 0 do not vanish. Physically this implies an axis of cylin- 
drical symmetry for the assembly (for example, the direction of 
an applied magnetic field in a gas of paramagnetic atoms). If 
there is no such axis there will always be tensors with Q < 0. 
If the assembly is isotropic, so that p does not depend on M or 
M’, (6.65) and the orthogonality of the Clebsch-Gordan 

{ Some authors restrict the use of ‘polarization’ to systems with odd-order 


polarization, using ‘alignment’ for even-order. Dipolarization is also called 
‘vector polarization.’ 
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coefficients shows that only the monopole (K = 0) tensor does 
not vanish, and the assembly is said to be unpolarised. 
If some operator O is expressed as asum of spherical tensor 


operators 


its expectation value in the assembly is found to be 


(O) = tr(pO) 
Bal ; , ea 
=>($eq) AxepKaled, aJ’)(a'I'||O xllat). (6.68) 


The sum runs over K,Q, «, «’, J, J’. We see that all dependence 
on the orientation of the coordinate axes (that is, on Q) is 
thrown into the statistical tensors pxg. It is this which makes 
these tensors such a convenient way of describing the assembly. 
Further, we see from (6.68) that the expectation value of a 
multipole operator O, of rank K depends only on the statis- 
tical tensor of the same rank. 

When the systems of the assembly are composite, so that 
the angular momentum J is the resultant of two or more com- 
ponent angular momenta (for example the spin and orbital 
momenta of a particle), the statistical tensors will show a 
corresponding structure. They will be expressible in terms of 
the tensors for the component parts, coupled in a way closely 
analogous to the matrix elements discussed in section 5.3 


[6], [19]. 


6.4.3. Development in Time 

The statistical tensors are a very convenient way of 
expressing the ‘angular information’ contained in the assem- 
bly. Now we look at the changes that can occur as the system 
evolves from time ¢, to ¢,, such as in a nuclear reaction or 
radioactive decay. These are described by the unitary trans- 
formation induced by the time-development operator of 


_ Dirac [20], {t) = U (tyto)|to) 


ie e(t,) = U (tty) (to) U (tito) *. (6.69) 
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As a special case we have the S-matrix [33] or collision matrix 
[44], S = U( o,—oo), connecting the initial and final states 
of some reaction or decay process. In an interaction represen- 
tation, for example, U can be written 


4 
Bey) = exp] ~6 i wat 
fo 


where H’ is the perturbing interaction. 

The important thing for our purpose is that for isolated 
systems H’, and thus U or S, is scalar under spatial rotations 
(conservation of total angular momentum). We establish 
from this that the statistical tensors pz, transform in the same 
way (6.69) as p, and that the tensor rank is conserved, 


pxalBJ, BI"; 4) = ~ Pral%J, a'J’; ty) U5 tito) Ugiqr(tato) *.(6.70) 


This means that the polarization or angular complexity of the 
final system (measured by the maximum rank of tensor 
required for its description) can never exceed that of the 
initial. For example a nuclear reaction initiated by s-waves 
will always display an isotropic angular distribution unless the 
colliding particles are polarized [19]. 

Perturbation theory is often used to calculate transition 
probabilities, in angular correlation problems for example [8]. 
This corresponds to an expansion of U or S in powers of the 
perturbation H’, so that the spherical symmetry remains and 
the transformation (6.70) still holds. 


CHAPTER VII 


GRAPHICAL METHODS 
IN ANGULAR MOMENTUM 


In any angular momentum coupling problem it is necessary to 
evaluate expressions containing sums of products of Clebsch- 
Gordan coefficients or Wigner 3j-coefficients. Such calculations 
may often be simplified by using the graphical methods of 
Levinson [85]. These methods have been extended by a number 
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of authors and the results are collected in two books by Yutsis, 
Levinson, and Vanagas [86] (this reference is denoted by YLV) 
and by Yutsis and Bandzaitis [87]. The present chapter gives a 
brief account of the graphical method which is complete enough 
to be used for solving simple angular momentum coupling 
problems. 

The rules presented in this chapter for constructing and 
manipulating graphs are not identical to those given by YLV, 
but are related to them in a simple and well defined way (cf. 
section 7.2.1). Our rules are somewhat more flexible than those 
of YLV; they may be used to evaluate algebraic expressions 
involving both 3j-symbols and Clebsch—Gordan coefficients, 
while the YLV rules can be used only for formulae expressed in 
terms of 3j-symbols. The graphical methods of YLV are similar 
to those introduced by Edmonds [22] and Judd [84], but the 
former give the correct signs as well as the magnitudes of 
expressions, while the latter give only the magnitudes. Exten- 
sions of the graphical method to include tensor operators and 
rotations matrices have been given by several authors [82, 83]. 


7.1. The Basic Components of the Graphical Represen- 
tation 

A graphical representation is a correspondence between 
diagrams and algebraic formulae. Each term in an algebraic 
formula is represented by a component of an appropriate 
graph. In a consistent graphical representation it must be 
possible to write down the algebraic formula corresponding to 
a given diagram in a unique, unambiguous way. 

The basic components of our graphical representation are as 
follows. 

(1) The Wigner 3j-symbol is represented by a node or vertex 
with three lines joined to it. These lines stand for the angular 
momenta which are coupled by a 37-symbol. 


eal (a) 
(Cay) 
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It is convenient to denote the orientation of the node by a sign; 
an anti-clockwise orientation is denoted by a + sign and a 
clockwise orientation by a — sign. Rotating a diagram does not 
change the cyclic order of the lines. The 3j-coefficient has simple 
symmetry properties (cf. Appendix I) and remains unchanged 
by a cyclic permutation of the columns in the symbol., There- 
fore, a rotated diagram represents the same 3j-symbol as the 
original diagram. The angles between the lines and the lengths 
of the lines have no significance. Consequently, any geometrical 
deformation of the diagram which preserves the orientation of 
the node does not change the 37-symbol represented by the 
diagram. A deformation which changes the cyclic order changes 
the orientation of the node and if the deformed diagram is to 
represent the same 37-symbol then the sign of the node must be 
changed. The symmetry relation 


ag) PEGs 


an ax 


implies that 


Ci (— )atd+e x SM 7) 
bB bf 
(2) The anti-symmetric symbol or ‘metric tensor’ 
(%,) = 8 - Ay 


which carries a phase present in many angular momentum sum- 
mations is denoted by a line with an arrow on it: 


aa bp 


='8(ab)/a 
| “) 
In particular (7.3) 
an a,—« =(—)te aa a,— “= ( ja~a 
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The 3j-coefficient reduces to the anti-symmetric symbol when 
one of the angular momenta is zero: 


(i ; i} — (2a+1)- 6(ab) (52) (7.4) 


This relation is represented graphically by 


bp 0 
(7.5) 


(3) An undirected line (a line with no arrow) represents the 
expression 0(a 6) d(« f), 


aa 6B (ab) (af). Woe 


(4) More complicated diagrams may be constructed by join- 
ing the three basic components together. Two lines representing 
the same total angular momentum can be joined. Joining two 
lines implies that the z-components of the two angular momenta 
should be set equal and summed over. It is not necessary to 
write these z-components explicitly in the diagrams and we will 
omit them. Often we will not even write the z-component of an 
angular momentum corresponding to a free end of a line 
explicitly, and will assume that the Roman letters a, b,c ... 
denoting total angular momenta have z-components denoted 
by the corresponding Greek letters a, 6B, y.... 

Lines which join nodes are called internal lines. External lines 
have one end connected to a node and one end free. Closed 
diagrams have no external lines. We will illustrate these defini- 
tions by constructing a few simple diagrams. 

(a) The first orthogonality relation for 37-symbols 


» é ; ‘) (: ; 3 _ sey d(ce’) d(yy") (7.7) 


ap 
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is represented by the graphical equation 
b 


a (7.8) 


The cyclic order of angular momenta (a } c) is the same in both 
3j-coefficients in equation (7.6). This order corresponds to an 
anti-clockwise orientation of the first node and a- clockwise 
orientation of the second. Hence the first node has a positive 
sign, while the second has a negative sign. 

If we put y = y’ in equation (7.6) and sum over y we get the 


result 
—- = (7.9) 


The signs of the nodes have been changed because their orienta- 

tions have been reversed. The factor (2c+1)-! has cancelled 

because the summation over y contains (2c+1) equal terms. 
The second orthogonality relation for 3j-coefficients 


abe\fabe 
y I —_ , f 
Hy 5 (cr pry) = Sax) 8108) (7.10) 
is represented by the graphical equation 


aa ad 


S (2c+1) > —= = 


c bB....-~SC« 
bp bp’ a ie (7.11) 


(b) The contraction of a 3j-coefficient with an anti-symmetric 
symbol is represented by a node with one arrow 


aa aa 


aa 


mess) 2086) < — 


(7.12) 
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This graph is very useful because it gives a way of representing 
a Clebsch—Gordan coefficient. This advantage is the main reason 
for modifying the original YLV scheme. A Clebsch-Gordan 
coefficient is related to the 3j-symbol by equation (3.3) 


(abapley) = (—yrre(aet ah (9°). (7.13) 


Comparing equation (7.12) with equation (7.11) we get two 
graphical representations of the Clebsch—Gordan coefficient: 


a 
<abaBley> = (—)8--¢(2e+1)* x p> 
b (7.14) 
a 
= (—)?4(2e+1)% x ya 
b (7.15) 


In the last equation the sign of the node has been changed. This 
corresponds to a change in the cyclic order of the angular 
momenta in the 3j-symbol and gives a factor (—)#+°+¢ because 
of the symmetry relation (7.2). 

(c) As a final example we give graphical representations of 
the Racah W-function and the Wigner 67-symbol. 


W (abcd; ef) = 
Se al earl Se | eae 


(7.16) 
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a ee lb eC ballad 
¢ a Cua) 


The sums in equations (7.16) and (7.17) are to be taken over all 
magnetic quantum numbers. Graphical representations of the 
9-7 symbol are given at the end of the chapter. 


7.2. Simple Rules for Transforming Graphs 


A calculation using the graphical technique often contains 
the following steps. 

(i) An investigation of a physical problem leads to a formula 
involving sums over Clebsch—Gordan or 37-symbols. 

(ii) The formula is represented by a graph using the rules 
defined in section 7.1. (Sometimes it may even be possible to go 
directly from the physical problem to the graph without writing 
down the algebraic formula.) 

(iii) The graph is transformed using rules derived in this 
section and in section 7.3. Each transformation of the graph 
corresponds uniquely to some algebraic manipulation of the 
formula. The aim of the transformation is to isolate various 
components of the graph which may be identified with standard 
invariant functions such as the Racah W-function, the 6)- 
symbol or the 9j-symbol. 

(iv) The transformed graph is reconverted to an algebraic 
formula. This can be done by comparison with standard graphs 
for 3nj-symbols. Often the last manipulations of a graph involve 
adding or removing arrows or deforming the diagram to get it 
into some standard form. These operations are discussed in the 
present section. 
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The rules for adding or removing arrows in a graph are: 

A. A line with two oppositely directed arrows is equivalent 
to a line with no arrows. This is a graphical expression of the 
result. 


aie )(3) = DA ema yae =e) fan") 
, (7.18) 


aa Aa aa aa 


B. A line corresponding to an angular momentum a with 
two arrows in the same direction is equivalent to a line with no 
arrows times a factor (—)?¢. 


Delors) (ge) = Dal = Irate? le —B) 88a") = (—)2% Ba!) 


f c 


aa aa ( ian y2 aa aa 


(7.19) 


C. If an arrow on a line with angular momentum a is reversed 
the graph must be multiplied by a factor (—)2¢. This result 
follows from the symmetry relation 


peal) 


aa aa (= )2a aa aa 


(7.20) 


D. Three arrows may be added at a node, one to each line 
joined to the node, without changing the value of the graph 
provided the arrows are directed either all away from or all 
towards the node. 


aD 


(7.21) 
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The proof of the first of these relations is 


Ble e Mewd lap) y) = rere) 


_ (’ b ’) 
apy 
because «+f8+y = 0 and the 3j-coefficient has the symmetry 
property 
a be ie b ‘| 
— (—)a+bd+e i (7.22) 
(a) = aa 


The second relation may be proved from the first and rule C. 

The number and orientation of arrows in a graph may be 
changed in many ways using rules A-D. A graph is in normal 
form if there is exactly one arrow on every internal line. YLV 
have shown that only those diagrams which can be put into 
normal form represent formulae arising from coupling of 
angular momenta. A diagram with no external lines represents 
an invariant 3nj-symbol only if it can be put into normal form. 

The normal form of a graph is not unique as the directions of 
arrows may be changed in many ways without altering the 
value of the graph. 

EK. The direction of all arrows and the signs of all nodes may 
be changed simultaneously in a closed diagram without altering 
the value of the diagram. Let J be the sum of the total angular 
momenta of the internal lines. Reversing the direction of all 
arrows in the diagram gives a factor (—)2/. This result follows 
from rule C if the graph is in normal form, because then there is 
exactly one arrow on every line. Adding or removing arrows by 
rules A, B, or D does not change the value of the graph, hence 
it holds for any graph which can be put into normal form. 
Reversing the sign at a vertex (abc) gives a factor (—)#+o+¢ 
(equation (7.15)). In a closed diagram each line is connected to 
two vertices. Hence changing the sign of every vertex produces 
another factor (—)?7. This cancels the factor (—)27 coming 
from reversing the direction of all arrows. 
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A graph may be deformed in any ‘way without altering its 
value provided 


(i) the direction of any arrow relative to the nodes it con- 
nects is not changed, 

(ii) the sign of a node is changed if the cyclic order of the 
angular momenta at the node is reversed. 


7.2.1. Relation with the Graphs of Yutsis, Levinson, and Vanagas 


The YLV graphs have an arrow on every line. A closed 
diagram in normal form of the type discussed in this chapter is 
completely equivalent to a YLV graph. Geometrically similar 
closed graphs in the two schemes represent the same algebraic 
formula. 

A YLV graph with external lines can be converted to a graph 
of the type discussed here if 


(i) any internal line of the YLV graph is left unchanged, 
(ii) an external line with an arrow directed out of the YLV 
graph is replaced by a line with no arrow, 
(ili) an external line (ax) with an arrow directed into the 
YLV diagram is retained with the arrow and the graph 
is multiplied by a factor (—)#-=. 


The graphical method may be used to simplify angular 
momentum formulae at two different levels: 


(a) if care is taken with the directions of arrows and the signs 
of nodes, the method gives the correct magnitude and 
sign of the result, 

(b) if reductions are made without worrying about signs of 
nodes or directions of arrows the graphical method gives 
the correct magnitudes. In some problems where the sign 
of the result is not important or where it can be deter- 
mined by physical arguments the magnitude of the result 
is all that is required. 


It should be emphasized that any calculation made using 
graphical methods can also be made using conventional alge- 
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braic techniques. To every graphical reduction there is a corre- 
sponding algebraic reduction because of the correspondence 
between graphs and algebraic formulae. The graphical 
method has two advantages over the algebraic method: 


(i) the notation is more compact because the redundant 
magnetic quantum numbers need not be written ex- 
plicitly, and 

(ii) reductions can be made by recognizing geometrical 
patterns. 


We conclude this section with some illustrations of applica- 
tions of the results presented so far. 


Example 1 

Verify that the graphical representation (7.16) of the Racah 
W-function is equivalent to the definition in terms of Clebsch— 
Gordan coefficients given in equation (3.13). 

If equation (3.13) is summed over the magnetic quantum 
number y then we get 


{(2e+1)(2f+1)(2c+1)2}! W(abcd,ef) 
= } (abpalee) (eded|cy ) (bd Bd] fd ) (afag|cy ), 


where the sum is taken over all magnetic quantum numbers. 
We write this equation as a graph using the representation 
(7.18) for each Clebsch-Gordan coefficient. 


{(2e+1)(2f+1)}% JN x [(2e + 1)(2f+ 1)(2¢ + 1)?]” 
(2c+1)W (abed;ef)= — = 


x ( = )2a+2e+2b+2a 


(7.23) 


The diagram in (7.23) can be reduced by omitting the oppositely 
directed arrows on the line c (rule A). Reversing the arrow on 
line e and inserting two arrows in the same direction on line } 
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cancels the factor (—)?¢+2>, Hence we get after cancelling the 
square roots 


Wabcd;ef)=<« 


In the first step of the reduction in equation (7.24) arrows 
directed out of the vertices (abe) and (afc) are inserted using 
rule D, and then pairs of oppositely directed arrows are removed 
from the lines a,b,e,f by rule A. These operations leave arrows 
directed in a clockwise sense on lines c and b and negative signs 
on the vertices. Finally the signs of all the vertices and the 
directions of all arrows are reversed using rule E to obtain the 
graph of the W-function in equation (7.16). 


Example 2 


Prove the symmetry relation 
W (abcd ;ef) = (—)+*-e-t W(aefd;bc). 
First draw the graph for the W-function given in equation (7.16) 


W (abcd; ef) = + 


(7.25) 


The diagram has been deformed so that e and f form two 
opposite edges. The deformation changes the cyclic order at the 
nodes (afc) and (abe) so the signs of those nodes must be reversed. 
The arrows may be transferred from the lines b and c to the 
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lines e and f by inserting arrows directed into the node (afc) and 
out of the node (aeb). Pairs of arrows in the same direction on 
lines a, b, and c can then be removed giving a factor (—)24*20+2¢, 
Hence 


. — { — \2a+2b+2 
W (abed; ef) = (-) ar 


In the final step the signs of the vertices (afc) and (abe) have 
been reversed giving an additional factor (—)@+f+¢)+(@+e+e) and 
the desired result obtained by comparison of the last diagram 
with equation (7.16). 
Example 3 

Prove 

abe 
W(abcd;ef) = (—)atorera ‘i ; A 

where the W-function and the 67-symbol are defined in equa- 
tions (7.16) and (7.17). 

The result can be proved by the following sequence of 
graphical transformations. 
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an é + 
= (— )atbteta We 7 


c a 


+ 
A phase (—)#+o+etc+at+e comes from reversing the signs of the 


nodes (abe) and (cde) and a factor (—)2¢ from changing the 
direction of the arrow on the line e. 


Example 4 
The following graphs all represent the same 9-7 symbol: 
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7.3. Theorems on Block Diagrams 


This section begins with a discussion of certain generalized 
angular momentum coupling coefficients constructed from 3j- 
coefficients and anti-symmetric symbols which have been called 
jm-coefficients by Yutsis, Levinson, and Vanagas (YLV). Let 
di(jimi), ... dn(jnMn) be the wave functions of the components 

Hil own JE 
WSL 6 6 6 | 
is called a jm-coefficient if it couples the angular momenta of the 
states $1 ... dn to a zero resultant, that is if the state 


of a quantum mechanical system. The function F'n ( 


oe DF n(,o °° 28) dalsums) «+ boldnm) (7.26) 


is a scalar invariant. In general F', will be a sum of products of 
3j-coefficients and anti-symmetric symbols which would be 
represented graphically by a diagram with 7 external lines. The 
detailed internal structure of the graph is of no importance for 
the questions considered in this section, and it is convenient to 
denote the graph of F, by a block with n external lines 


JnMn 


IV (7.97) 

Conditions which an expression F’, must satisfy in order that 
it should be a jm-coefficient have been investigated by YLV. 
They may be stated in terms of the graphical representation 
presented in this chapter in the following way. 

The expression /’, is a jm-coefficient if, by using rules A-D 
of section 7.2, its graph can be put into a form in which every 
internal line has exactly one arrow on it and every external line 
has no arrow on it. 


The anti-symmetric symbol (,) and the 37-symbol ( oa 


ap apy 
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are the simplest examples of jm-coefficients. General jm- 
coefficients retain some of the properties of these simple func- 
tions. For example 


(i) The jm-coefficient 


Fa J) la ‘a — 4) (7.28) 
M1... Mn 
unless }j; is integral and Yim; = 0. 
(ii) 
PF, dh oe: ) = (—yietia a DTS i). (7.29) 
—M1... —Mn M1... Mn 


The following theorems hold for jm-coefficients Fy, with 
ma— 12, and 3. 


THEOREM I. If r(? is a jm-coefficient associated with a 


graph with one external line then the function is zero unless 
— 7 — 0 


F(2 = F,(}) 8(j0) 6(m0). (7.30) 


THEOREM II. If Pa a ; "| is a ym-coefficient associated with 
12g 


a graph with two external lines then 


an je = ] Pos ji Bs. pe 
P(J) ~ (29, +1) 6(j1J2) (ona) t (7.31) 


a ji Seg 
where a 2 (ages x F223") 
Tunorem I. If Fa(? ua | 2s Siw ech able cme eee 
with a graph with three external lines then 
Fa(2” je rl ss (2 je i) Fs (7.32) 
mimems3 MimMeMs 
= ji je ds). p( J 32 i). 
where i ys (23) x a Saal? 


Mm, MMs 
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The proof of Theorem I follows directly from the definition 
(7.26) of a jm-coefficient. The only case in which one angular | 
momentum (jm) can be coupled to give a zero angular momen- 
tum occurs if 7 = 0. Theorem II can be reduced to a special 
case of Theorem I by coupling the angular momenta (j171) and 
(jamz) to a resultant (jm). Graphically 2 is represented by a 
block with two external lines. If these lines are coupled using 
the relation (7.11) we get 


(Arrows have been added at the vertices by rule D and removed 
from the line j by rule A in order to get the left-hand part of 
the graph in normal form.) By Theorem I this graph is zero 
unless 7 = 0. Hence, using relation (7.5) and cancelling re- 
dundant arrows, 


Fa( 2! 8 (jade) 1D: = 
m1 M2 (29, +1) 


which is the graphical form of equation (7.31). Theorem III can 
be proved by reducing it to a special case of Theorem II in a 
similar way. Details of this proof will be left as a problem for 
the reader. 

The properties of jm-coefficients discussed in the first part of 
this section may sometimes be used to decompose a complicated 
graph into simpler components. Suppose a graph can be 
separated into two blocks F and G where F has either one, two, 
or three external lines and is joined to G by these lines. It 
follows from Theorems I, II, and III of this section that such a 
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graph will split up into a product of disconnected components, 
one associated with the block F and the other with G. If the 
graph has one external line and is connected to G by this line, 
then Theorem I shows that the angular momentum associated 
with this line must be zero. The graphs F and G@ are dis- 
connected by simply omitting this line and using relation (7.4) 
at the associated nodes in F and G. The cases where F’ has two 
or three external lines will be considered in more detail. 

Theorems II and III can be used only if the graph F repre- 
sents a jm-coefficient. This condition is satisfied automatically 
if the graph F is in normal form with an arrow on every internal 
line and any arrows on the lines joining F' to G' are incorporated 
in G. If these conditions are satisfied, and if F’ and G are con- 
nected by either two or three lines, then the graphs decompose 
as follows. 


(a) Two connecting lines 


(7.34) 


These results follow directly from equations (7.31) and (7.32). 
In applications it is not necessary to put F explicitly into 
normal form. It is sufficient to be certain that any graphs 
formed by decomposing a larger graph can be put into normal 
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form by using the rules A-D. To satisfy this condition it is 
sometimes necessary to add pairs of arrows on the connecting 
lines before breaking them (cf. problem 6). We shall illustrate 
the decomposition (7.34) by an example. 


Example 5 
Prove the relation 


> (Pies) Sere iar-) a 
—feedea a: 0 


The result can be proved graphically by using (7.34) in the 
following way: 


Graphs which can be separated into blocks F and G where F 
has three or less external lines are called reducible graphs. We 
have seen that such a graph splits up into a product of factors. 
Trreducible graphs cannot be reduced to a product of factors, 
but often they can be written as a sum of products of simpler 
components. We will discuss the most important case where a 
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graph may be divided into blocks F'4 and G where F'4 has four 
external lines. The block F'4 will be assumed to have an arrow on 
every internal line and arrows on the four lines joining F'4 to G 
will be associated with G. If these conditions are satisfied, F'4 
represents a jm-coefficient. The angular momenta associated 
with the external lines of F'4 are denoted by a, b, c, and d. The 
graph may be decomposed by coupling the angular momenta 
a and 6 to a resultant x and c and d to a resultant y by using 
the graphical relation (7.10). 


(7.37) 


The result (7.37) follows from (7.36) by using Theorem II or the 
equivalent graphical decomposition (7.33). Arrows may be 
added to vertices using rule D of section 7.2 if desired and the 
signs of the nodes may be chosen in several different ways. We 
illustrate this result by an example: 
Example 6 

The jm-coefficients 


cary ep amare ere 
Hop = OUG e aMl4 7 a 


represent two ways of coupling four angular momenta to a zero 
resultant. Find the relation between the two coupling schemes. 


and 
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Drawing the graph for # and using the decomposition (7.37) 
we get 


a 
Q 
i 
“| 
= 
+ 
a 
>< 
Some 
x 
+t 
ore 


a 
«ESE err Dees 
(7.38) 
Equation (7.38) is another of the familiar relations between 
6)-coefficients and 37-coefficients (Appendix IT). 
An important special case of the result (7.37) arises when a 
sequence of blocks each with four external lines is connected as 
in the following graph. 


We assume that every internal line in F, G, ... has exactly one 
arrow on it, and that arrows on connecting lines are shown 
explicitly. Hence each of the blocks F,, G, ... represents a ym- 
coefficient and the diagram D is in normal form. 

This graph reduces to the following sum: 


D ae (2y+1)F(y)G(y) ... K(y), (7.39) 


where 


and similarly for G(y), etc. The result may be proved in the 
same way as (7.37). 

The results of this section are illustrated by a graphical proof 
of the Biedenharn-Elliott sum rule [5], [23]. 
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Example 7 
Consider the graph 


This graph may be decomposed by separating it on the lines 
(ghj) using the result (7.34) 


(7.40) 


See 


It may also be decomposed by using the result (7.39) 
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The 6j-coefficients can be extracted from the graphs in (7.42) by 
comparing them with the standard graph (7.17): 


D= 2 (Qe 1)(—jerorerdverzegeneten Me ; ane : iy 
(7.43) 


For example the first graph in equation (7.42) has a value 
(—yreeal y AR We use the facts that g+A-+7 is an integer 


and (—)-3* = (—)* in obtaining the final form of equation 
(7.43). 


Example 8 
Evaluate the expression 
P(kq,k'q’) =  (—)h mtd (' —mm'|kq) JI’ —NN'|k'q' ) 
x (lymm,| LM ) ilemyme|KQ) (Liz Mm2|JN ) DK 
x Um'm,|L'M’) Ulimyms|KQ ) (LL M'm3|J'N’ ), 
where the sum is taken over all magnetic quantum numbers 


except for g and q’. 
Using the relations (7.12) and (7.14) this expression reduces to 


F = Fy(—){(2L+1)(2L' +1)(2J +1)(2J’ +1)(2k+1)(2k' +1)}* 
x (2K+1) 

where 

= 2-2 oO ope oy op no 


(the phase (—)* = (—)?!+24 because (JL1,) and (l’L'l;) satisfy 
triangular conditions) and F; is represented by the graph 
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From Theorem ITI the graph for F may be separated on the 
lines J, J, and K, and on the lines /’, J’, and K giving a product 
of three factors G1, Go, and G3. These factors are 


= (—)2Y W(L, lol J; KL). 


_ 


—)h 
2k 
ihe hy 22 Dated ANT yaar en 


W(WII'skK) 5(kk’) 6(9q'). 


10 
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Collecting the results we get 
F(kq,k'q') = WhlelJ; KL) WLU; KL’) WWI RB) x 
x O(hk’) 6(gq’)(—ttK-#(2K -EI)[(2L+41)(2L' +1) x 
x (27 (27 =e). 


APPENDIX I 
3-7 AND CLEBSCH-GORDAN COEFFICIENTS 


THE CG coefficient is defined by the transformation (2.30) 
labey) = > labaB)(abap|cy) 
af 


and vanishes unless «+f = y. Other authors use the notations 
(abaplabey) [17], (axbplabey) [22], C.,(cy; «B) [9], Clabc; «B) 
[54], Ce [65], C2, [39], Si [78], and C%, [42] for the same 


quantity. The Wigner 3-7 is related to it by 


— _ ¢__\a-b-y 4 (@ b Cc 
(abuplo—y) = (—*%2e4 uy (FF ‘}. 


Note the appearance of y with a minus sign on the left, so that 
now a+f+y = 0. Related quantities are defined by Racah 
[31] 


ab ¢ Ce 
v( = V(abe; = (— | i: 
es (abe; aBy) = (—) - 
Orthogonality relations: 


> ec-+1)(% 


af 


> ec+n(4 } aE _ l= as bgy 


cy 


~ D om 
~~ 
ee 
a 
R & 
Do 
So ek 
— 
or) 
° 
Ge 
mo 
= 
Ss 


Symmetry: If the 3-7 is rewritten [51] 


(: b :) fae ct+ta—b cr | 
=| a—-« b—B oa ||. 


ae ata b+ 8 e+y 
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it is invariant under interchange of rows and columns (reflec- 
tion about diagonals), and is multiplied by (—)*+®+¢ upon inter- 
change of two adjacent rows or columns, giving 72 equivalent 
symbols. In particular this means 


& S {y= i: é | 


and that the 3-7 is invariant under cyclic permutation of its 
columns and multipied by (—)*+°+* by non-cyclic ones 


5 Ib Deere td 


Recurrence relations: 


me |  Oeeee 
UcFrety+UM yg yaa) t 


s(t eae” 
+[laF a)(ate+I1)] (om a ogee 


b 
(OFA o+b+0r(% B41 “| =o, 


[latb+e+iyb+e—ay (% % {)=(0-+8xe—rIx 


—4 —4 b—3 =) 
10 a a -10—Aye-ron (( p+} amy 
,(@ b ctl 
((at+b+c+2)(b+c—a+1)(a+e—b+1)(a+b—c)] oe 
,(@ b c 
= [(6—B)(b+B+1)(e+yVe+y FIV). 
b 
—2ptlety+ite—y+ut (% & 2) W+Ax 
, (4 b c 


Algebraic formulae for the general 3-7 are given by equation 
(2.34), and for « = 6 =0 by equation (2.35). Formulae for 


10* 
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c = t, l, 3, 2 are given in Table 3. Other special cases are 


abo ao ‘ =) 
(; B J (a 0 a Ly 


(” b aed 

Gee) Gy 

hae (2a)! (2b)! (a +b+y)! (atb—y)! : 

= ee sorscy (a+a)! (a—a)! (b+8)! a 

€ ee Il 
aw y 

(2a—1)! (26-1)! (a+b+y—1)! (a+b—y—1)!9! 

| (2a +26)! (a+«)! (a—a)! (6+ 8)! (b—8B)! | 

(: bab — 2 
a fp y 


= (—)*°-72(ba — af) x 


ore: 


(2+b—y—2)! (a+b+y—2)l (2a—2)! (262)? 
= (a+a)! (6+ 8)! (6—8)! (2a+26—1)! 


x [a-+a)(a +a—1)(b—f)(b—B—1) +(a—a)(a—a—1) x 
x (b+A)(b+B—1)—2(a —a)(a+a)(b +6)(b —B)), 


a bo ¢ eta 5 c=) 2A(abc) 
( (~) [(2a1)(2b-+1)]} * 
[}(k-+a+8))! 
Hatb—h]! [Ha+k—b—1)} (O+k—a—1)} 


ai 


with k =c if a+b+e even, k =c+1 if a+6-+c odd, and 
A(abc) as in (2.34). 


(2c)! (a+b—c)! (a—a)! (c+b+a)! a 
(a+b+c-+1)! (c—a+b)! (c+a—b)! (b—c—a)! aa 
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By specializing some of the relations in Appendix II: 
Gee Ge EN. 
(20-+1)( = (20-+1)( =} 
Pi pa 2 a? 


2 
a 


SS aya » — _e(¢+1)+6(6+1)—a(a+1) 
DA2+0(, g 5) = Bra. =7 : 


2c(e+1) 
Ce 0a a b c\c(c+1)—a(a+1)—b(b+1) 
( = 0)=(0 0 ) 2[a(a+1)b(b-+1)}* 


if a+b-+tc even, 
Geueye e ete (c—2)!q# 
(1 —3)=(1 1 JO —anaro-+0 sy) 
if a+b-+c odd. 
Gane a— la ie fee") 
( 0 0) =( Sell ) c(e-+1) 


“ar iia 5 eetee ty 


aa—le 3 (a+1)(a—1) i 
( ll 0} Ss 


aa Nl c(e+1) 


I 


+ 
000 | if c is even. 


= (c—1)(¢ +2) 
ie > lee es Vaan 
i; b i)= -a(} j ee 

a c(c 
eee oe | 
= 2 1 

ote =: dS 
fe = a+2 ) eo Bess eee 


1 
b 4 0 


e(c+1)(2a+1)(2a+c¢+43)(2a—c+ 2) 
if c is even. 


1 
| 
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Clebsch-Gordan coefficients, symmetries: 
(abaBley ) = (—)**?-* (ab—a—Ble—y), 
= (—)#+o-¢ baBaley ), 


2e+1)\?, \. _ 
= (Fey) (Hie ace —ylb—B), 


= (223) (-pnetcrtes 


Special cases: 
(abaB|00) = (—)*-2(2a+-1)-*ba0da,-2; 
(a0a0|cy ) = Oac Ons 


APPENDIX II 
6-7 SYMBOLS AND RACAH COEFFICIENTS 


THE Racah coefficient is defined by the transformation (3.8), 


\(ab)e, d;c) = > |a, (bd)f; 6)[(2e +1)(2f+1)}* Wabed; ef). 
f 


Orthogonahty: 
¥ (2e+1)(2f+1) W (abcd; ef)W (abcd; eg) = 6,,. 


Symmetry: (giving 144 equivalent coefficients [38], [52]) 
W (abcd; ef) = W(badc; ef) = W(cdab; ef) = W(achd; fe), etc. 
= (—)+-¢-f W(aefd; bc), ete. 


also = W(aBCd; EF), where 


B= }(b+ce+e—f),C = tb+c+f—e), # = tb+e+f—c), 
i ie --e-j 5) 
Sum rules [5], [7], [23], [48] 
W(abed; ef) = > (—)?*?-*(2g +1) W(gabf; de) W(gdbe; ac), 
W (abed; ef)W(abgh; et) = > (27 +1) W(jgfa; cc) W (gdb; hf) x 
j 
x W(jgde; ch) = (27 +1)(2k+1)Wcifg; ka)W(cijb; kh)x 
x W(gfjb; kd)W(cdhg; 3). 


jk 
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Contraction of 3-j symbols: 


W(abed; ef) => (2c + aie | ? si bie °) x 
d 
5 


summed over all z-components except yy; 


rentals) Bere JUHA D 


Deft (abcd; ef) (> s ia ‘| i4 ; ier 
ale “Nee 
D, 26+12/+1) W (abcd; ef) ke i i) i ; 3) (; 5 * . 
corolh 


> (20 +1)(2e-+1)(2f-+1)W (abed; al? : ‘) E j = < 


cef 
dc e\fab e€ Sig tay er ee 
x(% y We B Jie us 


If a+b-+e is even, a special case is 


d b 
(e } ‘\ = —[(2a+1)(2b +1)]* W (abcd; o)(6 0 Al 


Algebraic formula for the general Racah coefficient is given 
by equation (3.15), and for e = 4 and 1, in Table 4. Other 
special cases are 
W(abcd; a+b, f) = 

(2a)! (2b)! (a+b6+c+d+1)! (a+b+c—d)! 
= eaeoman I(c Tie [= al Coy (a+f—c)! 
(a+b+d—c)! (c+f—a)! (d+f—b)! ; 
“(a-+ot+f+) ( Toa Te d)! eral 
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(= ae” 


Meee P= ea $1) BeF 1) 


6(a, b) d(c, a). 


6-7 symbol, definition: 


‘i : A ae (—)etote+d W(abcd;ef). 


Triangular conditions: the four triangular conditions which 
must be satisfied by the six angular momenta in the 6-7 symbol 
nay be illustrated in the following way: 


(lal ecltile” Nalttlonc’ | 


Symmetries: the 6-7 symbol is invariant for interchange of 
any two columns, and also for interchange of the upper and 
lower arguments in each of any two columns, i.e. 


ena ferd—Geam {ere — te ssp e 


Contraction of 3-7 symbols: 


> (—)A+B+Ctataty E Be Ne vNs eae b i 
ey « —By'NB —y a’) \y —a 8") \a’ 8 ¥4 


abe 
= sn ates ay B of 


Sepnenim(d BYE CAYO 48) 


Sum rules: 


» (—1)24(2k-+1) 4 b sto 


k 


Pe) )at+bt+k/( (ebay t e Ar drof( (2a+1)( )(26+1)}, 
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> ChE NOPE Ge ag = 
Z J ea Ae om fs 


Zeng Nao of “1h o ah 


APPENDIX III 
9-7 SYMBOLS OR X-COEFFICIENTS 
THE X of Fano is defined by the transformation (3.23), with 


((ab)e, (de)f; i[(ad)g, (be)hs 2) , i 


[(2¢+1)(2f+1)(29-+1)(2h+4+1)} g fe 


= X(abe, def, ghi). 


Orthogonality: 


are \ (a bc 
S ve snernes ann e ie e i = bg4 Onk. 
cf dm OCP 1G 


Symmetry: (72 relations [39]). The X is invariant under 
interchange of rows and columns (reflection about a diagonal) 
and is multiplied by (—)? 

(where p = a+b+c+td+e+ft+gth +12) 
upon interchange of two adjacent rows or columns. 


Sum rule: (others are given in the literature [2], [61]) 


@ 0c) (anergy aoe 
Sim napsaneen} d if b i -{? e i} 


ik pkisig ht Goh 


Contraction of Racah coefficients: 


X(abc, def, ghi) = > (2k +1) W(aidh; kg) W(bfha; ke) W(aibf; kc), 
k 

> (2c+1) W(aibf; kc) X (abe, def, ght) = W (aidh; kg) W(bfhd; ke), 

hie, 
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W(ghjk; wl) X (abc, def, ghi) = > (28 +-1)(2t+1) W(efjk; ts) x 
st 
x W(desk; ft)W(belk; ht)X (abc, dts, glj). 
Contraction of 3-7 symbols: 
abe 
abcfbeh eed eee 
(2 1 
{ied wed (TB No cally d+ x 6 p}~ 
d 
6 


ght 


summed over all z-components except «. 


a One . 
CoN -BG IES ICS IG SA)» 


cy, gent gcaied 
ee ale ANS, I 
S e+nec+y(4 } j ie ae 5 yee . i 


be Y 
_f(ad ae Ne h ' 
= (6: ded uci 
Special cases: 


, ? / _ bap bpp =)P4-2* W(abdes og) 


d 
Sa [(2e+ 29+} 


dec a 
pif ele IVR Ng G+) 27+} 
x (—)*+9-2-¢ W (abde; cg). 
With g = 4, c-+d-+e even, and using A.1, this gives 


ale 
[6(2e+ 2d+1)(2e+1)}(5 | nM ; =() ik 
4 g—!}, 


fief a(a-+1) dd +1) —b(b-+1) +e(e+1) 


2 


we oO 
— 
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When c+d-+e is odd, we have two other relations 


ab ¢ 
c+lde 
( 0 0 ‘} ‘i “ah 
44 1 


= 


q a eae ed, (2 b ‘); 
- [6(c +- 1)(2c + 1)(2c +3)(2d +1)(2e +1)]* 


c—lde 7 im, 
Wiacoc Sank 


l 
easel Se ee De en (e >). 
a (Gees ay ( 221) (2a 27 1 )\(2e2 1 4 —2 0 


Algebraic formulae for g = A = } are easily obtained from 
these. 


APPENDIX IV 
SPHERICAL HARMONICS 


TueseE are defined with the same phase as Condon and Shortley 
ar Y= (2k+1/47)'C,,, where 


a 
' *] Pub, if g > 0, 


and 
Cy (86) = (—)*C,,(84)*. 


The P% are the associated Legendre polynomials, with 
P? = P,; their properties are well known [40], [42]. Special 
cases are 


(Oh) =e Che = ORR O14 = e(4)) cin Oe 
Co 


Crp = ome 


3(3 cos’6 —1); Cy. = (3)* cos 6 sin 6 e***; 
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Orthogonality: 
(2k-+1) | C,o(6$)*Cpo(04) sin 0 d0dp = bp Sy47. 
Sum rules: > |C,.(66)|? = 1, 
q 


¥ (2k +1)C,,(6¢) = 26(eos 0—1). 


k 
Addition theorems: >) C,,(0¢)C,,(6'¢')* = P,(cos a) 
q 
if w is the angle between the two directions (64) and (6’¢’). 


C4080, (06) =D c,09i2e+1y-¥ (2% (EG Gl) 


C 


2 C08) 0,9(98)( 4 : =) = c,(08)(—Y (8 a 


| .(68)608)0, (64) sin 6 dédé = i-(“ 5 s) (" : a 
abe 
000 


% 


. { HG) Gee 0) feos Py ctin 2 9 = 2( 


APPENDIX V 
ROTATION MATRIX ELEMENTS 


Our definition of 2?,,.(aBy) for the rotation of axes through 
(afy) is the same as that of Rose [54], Messiah [46]. Bohr and 
Mottelson [13] use rotation matriees whieh are the eomplex 
eonjugate of ours. Wigner [78], Fano and Racah [31], Edmonds 
[22] and Rose [53] use the same notation, but («fy) are then 
Euler angles for a rotation of the system. Associated with this 
is the equivalent eoordinate rotation (y B «)-! = (—x« —B —y), 
so their usage is equivalent to a different sign convention for the 
angles of rotation. 

With our definition (equation (2.17)) 


Bi Aapy) =e di (Bye ty, 
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The formula for d/,, Ae: is given by equation (2.18), and given 
explicitly for 7 = 4, 1, $ and 2 in Table 1. 


Orthogonality : 
(25+1)[ Diu* Dare sin PaBdady = b55 Syae Sqynee 80°. 
Sum rule: ps DAD in”! ay = Onan 


ae 


Closure: 2, Drm %2 2 2) Prem %By1) = Dw XBY), 
where (ay) is the resultant of first («,6,y,) then (a.Bo7o), 
» Y Fone(B2)Epem(Br) = Tane(Pr +B2)- 
Symmetry: 
Bye) = (=) big (B) = FE gy mB) = Evin —B) 
SP ET uted)! (i) 
Deo = (en Pe apy) — Dn = 8 — 
where (—y —f —a) is the rotation inverse to («fy). 
Special cases: 
Zino( By) = Cjm(Bax)* 


a ! $ 
. dig(6) = ("FEE | Pre) if m >, 


di(B) = P,(cos 8). 


di,(B) = (—1)°-™[(29)![(j +m)! (j —m)!]* x 
x (cos $8)7*"(sin 48)'~™. 


Contraction: 


4,(ap/\ DB (aby) =>. (00+1)(4 2°) (4 Vac iapy 
Cc 


a (apy)*(4 a =a i B as (xBy) Dor (XB); 


abe 
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Ip )*= > C0+(4 § oscsnairain(S 2 ), 
aba’b’ 

> Mob DB ep) E(aby)(4 PO) (4 2°), 

abe 


[ 22 of) PAaBy) Dol aBy) sin faBadady 


APPENDIX VI 
TENSORS AND THEIR MATRIX ELEMENTS 


THE commutation rules for spherical tensor components with 
the spherical components J,(u = 0, +1) of J are 


oo ae 

2 k+@tyu-1 

(Fur Peg] = Tags {HE +1)(2k-+1)}( k++ (Ri ny 

Basic tensors are the spherical harmonics C,. When k = 1 

we have vectors, a = aC,(6d) = > (—)*a,e (64) are the 
u 


BY — yt 
polar angles of a, so C, is a unit vector along a. The e, are unit 
spherical vectors, eg) = e,, e,, = ¥ (e,+ie,)/2*, so the vector 
components a, = aC,,(66). Product tensors are defined by 
(4.6), 


MAU, Si) == Baal, by) 


= 7 \2ke+K—Q K ky ke 

= (2K +1)*(—) >(¢ —% —GQ Bo Sisay 
When k, = k, = 1, /2T,(a, b) = iaa b, and T (a, b) is given 
in section 4.5. Other examples are the bipolar harmonics of 
section 4.6, and the spherical harmonic addition theorem. We 
may derive various re-coupling relations for such tensors; for 
example, if tensors k, and k, eemmute 


Tic(ky, ky). Tic(Kg, kg) 
= (26 4-1)(—)9** Y W(kykokgk,; KR) T (ky ks) . Ty(Kgk,). 
ra 
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For example, with K =0,k, =k, =k, =k, = 1, 
(a.c)(b.d) = }(a. b)(c.d)+3(anb). (ca d)+ 

+T,(a, b). T,(c, d). 
When c, d are Pauli matrices, with o.o = 3, T,(¢, 0) = 0, 
[20], 
(a.o)(b.c) = (a. b)+i0.(anb), 

and when a = b = Jj, c = d = J,, we get 


T,(J,Jj) . T (JoJo) a (J, . J2)? +4(J, . J,) —4JiJ5. 


Again, ifa = s,,b=s,,c =d =r, we get the tensor force 
Sie 


T,(S,S2) . T,(rr) = r°S,,. = (s, . r)(s, . r)—3(s, . s.)r? 
where T2(rr) = di (5) r2c2168) This may be recoupled, using 
r(s.r)—}r*s = —(10/9)*r?T,(C,,s) 
to give T,(s,8,).T,(rr) = —(10/9)+r?s,.T,(C,,s,). 


Tensors may be formed by ‘polarizing’ solid harmonics 
oC. with n < k vectors A, B, - - - [27], [75], 


TA, B, 089 H, Ta) _ >; (Bi y, 


Ape 


AB ee VV ee 


ib 


-- V_,1*O,,(64), 


and if B, -- - H commute with V, 
= (A.V)(B.V)---(H. V)"C,,(6¢). 


Each step replaces a vector r by a vector A, B::-, 
leaving unchanged the transformation properties. Other 
polarized harmonics may be formed with the operators L and 
V; vector harmonics were introduced in section 4.10.2 and are 
further discussed by Hill [35] and Edmonds [22]. In the 
product tensor notation above, the vector harmonic is 
py rea 7 OY. = y_yk-ofX 

Ga = Tro(%e €) =>, 2K +MY Q | _ 

gu 
x Y,,(60)e,,. 
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Some useful properties of the various tensors follow from 


L,G,.(9¢) 


k ok 
= er yerro(—rrer( tA 1) Cha) 


auioeneiie | SRA) Spee he k 
psa ra96)(5——)st00 +] | eco 


r 


k(2k—1) 
2k+1 


k-1 k 1 rs F d k+l 
ase =F E Ree yi ait : J go, 
VA" C4) = (—jPtet*[k(2k—NQk+1)F x 
ie ee 
= ee me _,)? Cetato 


V,(r-F1G,,) = (=F tt +1)(2h +1) (2k +3)] x 


k+1 ik 1 \ 
(Ft = alr ‘ "Ch ateta 
So we have [kK(A+1)]}'YZ,, = L( i) 


Pa ee, = VE me) 
et a SI. 


Further, using the relation 


vats ~ ~[ee4-o( 2) 


we get 
Ook 
(2k+1)V A L(fY,,) = ik(e+1) Yhou($—-)f+ 
‘ ; d k-+1 
+1(k+1)k Y¥} k-ll (5 —) ae 


and in particular 
VaL(rO,,) = i(k +1)V(r*C,,). 


Similar tensors arising in radiation theory are discussed in 
section 6.1.1, and others have been applied to B-decay theory 
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[3], [45], [57]. We also have 
d k+1 
(2E+N) (FO) V = THe —P(S +=) for) Ty(Gy_.¥) 


ak 
—[(k-+1)(2k+3)}* (= ==) 400 P gl Gz 41) V)- 


Some of the multipole tensor expansions useful in physics 


are 


oe — (21-+1)i9,(kr)C,(0,¢,) . C,(8,4,) 
6(a—b) = Ea Niet 1) (20+ 1)C,(0,¢,) . C,(8,4,), 


e-v(a-b)” — Se (21+ 1)e-v(27+b") 5,(2iyab)Ci( Gada) . Cr(Ovd0); 


and if p = b—a, with b >a, 
Ip = > (a'/b'*1)C,(0,$4) « C,( 9.45) 
ely = ikh (kp) Sebel (21+1)j,(ka)h{?(kb)C,(0,¢,) « C,(9,4,), 


Ag) — (bel jdiaa)h{? (tab) C(8.b4) - G85); 


ifr =—a-+b 


a 2l! 4 
Ll a nish ok ola — 
r'Cim = ys (saaa—ai) ee 
x C1-a,m-u(Bapa)Cru( Boho) 0 —AAm —ppe| lm ). 


Reduced matrix elements, definition: 
The Wigner-Eckart theorem (4.15) states 
—)2 (JM |I'kM’g) (J||Txl| J’ > 


ee 


(JM|Tiy|\J'M’) = ( 


The reduced matrix element (J||Tx||J’) defined and used by 
Racah [31, 48] and Edmonds [22] is related to ours by the 


equation 
(A|Pel 7’) = (27 +1) IIITa|"). 
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The factor (—)** is included in our definition so that the phases 
of both reduced matrix elements should be the same. It is 
relevant only if k is half-integral (for example if 7'x¢ is a creation 
operator for a spin-} particle). 


Reduction formulae in terms of 6-7 symbols: 

For convenience we collect the principle reduction formulae 
of section 5.3 and rewrite them in terms of 6-7 symbols. In a 
two-component system the tensor Rx, (1) acts only on the first 


part and S,;,(2) only on the second part. If 
Pro(kike) = X Reya,(1) Skya,(2) (hikeqiqe| K@ ) 
Nn 


then formula (5.12) gives 


2 Ja ke 
(2.1) GallRag|li (292 +1)* (Jal Sual 32. 
Special cases of this result are 


S=1 , K=k, =k (equation (5.9)), 


J J’ 
(hide ||Takike)||Ii37") = (er ser sa ji ke x 


a — JJ’ k 
HidoF Re G95") = 6(9295)( 2 +05 j a x 
hi 1 2 


x (—)F ath (29, +-1)8 [RCL 9 
Ree a ee 
py ‘ A 
JAS 2M FGI’) = bj 2I +0 a a 
Js Je Ji 


X (— )e tier 2a(256 +-1)F (Jol[S.(2)I1 95); 
K=0 , kj) =ke =k (equation (5.13)), 
Rx.Se = (—)#(2k+1)! Too(kk) 


hrjoFlIRe- Sill 7557") = HII )(— yore J ; x 
jo Jo J 


X (291-1) jal[Ralldt (252+1)? Gell Sullge ). 
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Basic reduced matrices are 
GJ") = VU +1))'4,,., 
aload 
Eve £9 ’ $ _ yi 
aCe) =v + r(, Fs il 
(uel Ne 
¢ = , : 
avy = r= Kt =e e 
Fi 2 
= (—yran( 3) « 
aaa 
where the last factor is a radial integral and a = —l’ if 


el la lif b= 1 —1, 
The matrices of tensor products and for composite systems 
are discussed in Chapter 5. Application of (5.9) gives 


ee: 
asters) = even -y() FF) 

provided J+l'+k is even, and zero otherwise; 
: ee  k 
TAIT (Cy, oN") = axl -M4gtsyi(3_F 8) 


where Oe = (4 —2')/4/{k(k+ 1}, 

any = —(k-a-e')/V{R2k+ D}, 

ane = (kt 122 (RAN) 2R+)}, 
with aw == (L—j)(27+1), # = U—7')(27' +B. 


Waid ek = 6, 6, (— yiet1-J ha x 


Titi Jede 


xX(~lat (24, + V2T +E W Aids’ 1je), 


(Ade 


and for example, since 6 = 2s for spin 4 
(aylleolt337) = SylF+IG+D—U+DIGG+DI ifs =J', 
= 6,2(0-+1)/(21+-)F ifj =j'-1 
and (44S||o,||34S’) =[S(S+1)]}? if S=S’, 
= —/3 if Ss — 
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Equations (5.12) and (5.13) give 
GIA Sr - Felli ieF) = 45,55, ¢81F J +1) Alii +) Jae +, 
(445llo,.0,||345) = —3 if S=0,and1 if S=1. 
Also ($485||T.(¢,, ¢5)[|$48’) = dg5-65,(72)*. 
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ASYMPTOTIC EXPRESSIONS FOR LARGE 
ANGULAR MOMENTA AND CLASSICAL LIMITS 


The classical ‘vector model’ has already been introduced in 
connection.with the interpretation of the matrices for rotations 
of the coordinate system (p. 27), and of the vector addition, or 
Clebsch—Gordan, coefficients for the coupling of two angular 
momenta (p. 32). We would expect from the Correspondence 
Principle [46] that the various quantum mechanical expressions 
that appear in the treatment of angular momenta would reduce 
to the geometrical forms obtained from classical mechanics, 
when the magnitudes of the angular momenta become very 
large compared to h. It should be remembered, however, that the 
reduction of quantum mechanics to classical mechanics may 
not converge uniformly. The existence of quantal interferences, 
which have no classical analogue, can result in the probability 
of observing a given value for some variable that oscillates as 
that value changes. These oscillations grow ever more rapid as 
the magnitude of the variable increases. Then, in order to see 
convergence to the smooth behaviour predicted classically, one 
must average the quantal oscillations over a small interval. We 
illustrate this property here for the spherical harmonics 
Y,,(9, @) when / is large; see also the discussion of the rotation 
matrices 2,4, on p. 29. 

Another consequence of quantum uncertainties is that a 
variable may assume values that would fall outside the range 
allowable in a classical description. However, the probability of 
finding such a value decreases exponentially as the value moves 
further outside the classical range, and the decrease becomes 
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more rapid as the magnitude of the variable becomes larger and 
more ‘classical’. This phenomenon was mentioned on p. 29 in 
connection with the probability of finding a projection M’ of a 
vector J on a new Z-axis after a rotation of the coordinate axes; 
other examples featuring spherical harmonics are illustrated by 
Zare [98]. Such nonclassical behaviour is closely analogous to 
another well-known quantal effect, the penetration of potential 
barriers [46, 99] that would be insurmountable in classical 
mechanics. 

In this Appendix we give a few examples of asymptotic (‘large 
j) expressions for some angular momentum coefficients and 
some alternate forms they may assume in this limit. Many 
others are listed in the exhaustive compendium of Varsholovich 
et al. [96], who give other referenes and sources. (See also the 
early discussions by Brussard and Tolhoek [15] and Wigner 
(78].) Some of these asymptotic expressions are surprisingly 
good for quite small values of the arguments. Besides illuminat- 
ing the transition from the quantal to classical domains, 
asymptotic, large-) relations may find practical use in compu- 
tations. Often one has the problem of evaluating sums over 
angular momenta that may take on very large values; this often 
occurs, for example, in descriptions of atomic, molecular, and 
heavy-ion collisions [97, 98]. The sums may be convertible into 
integrals after the use of one or more asymptotic relations. 
Sometimes the integrals may be performed analytically, thus 
greatly simplifying the results. 


Spherical harmonics: Consider a particle which, classically, 
would be moving in a circular orbit with angular momentum of 
magnitude [((/+1)]'?h=(/+})h and z-component mh. The orbit 
(Fig. Al) would be inclined to the z-axis at an angle p= 
sin~(m/[l(i+1)]”). The polar angle 6 of the particle is limited to 
the range B<0<x—f. 

In quantum mechanics, the x- and y-components of angular 
momentum become indeterminate, as though the orbit were 
precessing about the z-axis (p. 28), and the orbit is replaced by a 
probability distribution. This distribution is proportional to the 
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Z 


= 7 


Fic. Al. Particle moving in a classical orbit with angular momentum 

J[i+ 1a (1+})h and z-component mh. The plane of the orbit is inclined at 

angle f to the z-axis, where cos B =m/,/[l(1+1)] ~m/(1+ }). In the vector model, 
the vector | is supposed to precess around the z-axis. 


square modulus of a spherical harmonic | VAG o)|”, which 
represents the probability of finding the particle with polar 
angles (0, @). The indeterminacy of the x- and y-components of 
angular momentum is represented by the uniform probability 
d@/2x of finding the azimuthal angle between ¢ and ¢+d¢. 

The m/l=0 case is the easiest to visualize. Then the z-axis 
passes through a diameter of the corresponding classical orbit 
(B=0 in Fig. Al), while the plane of the orbit is distributed 
uniformly around the z-axis, with the probability (2z)~' rad71. 
Classically, it is easy to see that the probability of finding the 
particle with a polar angle in the interval @ to 0+d@ is 
d@/(2n* sin 0). Consequently, we would expect to find the same 
result in quantum mechanics in the limit of large 1. 

Indeed, there is an asymptotic form for the spherical 
harmonic, derived for !— 00 (but 1> m0) which allows us to see 
this correspondence, namely 


emt T 
Y nO, @) & Saw snl +5)0+(2m+1) =| (A.1a) 


/sin 6 
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provided @ is not too close to 0 or x (@>1-! or m—0>I1~). We 
illustrate this for m=0, when (A.1a) becomes 


sn] (+-)0 4. _| 


TOO) = ee sin 0 


(A.1b) 
In fact, the ‘large /’ approximation (A.1b) is remarkably good for 
small / (even J=1!), provided @>1~'. This is shown in Fig. A2 for 
1=4 and /=10, where the approximation (A.1b)is compared with 
the exact values of Y,(0, @). 


! | | (Lele | 
0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 


@ (rad) 


Fic. A2. Illustration of the relation (A.1) when m=0, for /=4 and 10. The solid 

curves represent the values of Y,.(0, @), while the dashed curves are given by the 

approximation (A.1). Note that Y,.(x—0, b)=(—-)'Y, (8, @), and both are 
independent of ¢. 
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When @=7/2, approximation (A.1b) predicts that 


Wal 6) ee) leven 


Tt 


=), l odd. (A.1c) 


This is exact for odd /, and very close for even /. Even when 
1=2, this estimate differs from the exact value by only one per 
cent. 

When m #0, the range of small 0 for which (A.1a) is not a good 
approximation increases relative to that for m=0. For example, 
the limit (A.1a) is not at all accurate for Y,,,(0, 0) until 0>4, and 
for Y,,4(0, 0) until @>1. 


Classical limit for spherical harmonics: The last sine factor in 
the expression (A.1a) oscillates very rapidly with @ as 1 becomes 
large. This is an example of the quantal oscillations mentioned 
above. Consequently, the probability PAG, o)|? must be aver- 
aged over an interval A@ in @ which is small enough that the 
variation of the (sin 0) denominator may be ignored, but which 
contains many oscillations in the numerator; thus A@> 72/1. The 
the square of the oscillating sine factor may be replaced by its 
average value of 3. (Note that the quantal oscillations persist 
however large / becomes, so that the averaging is always 
required, in principle, in order to make the transition from the 
quantum mechanical result to the corresponding smooth 
classical one.) 

Thus, the averaged probability distribution becomes 


L¥,.(0, &)°| m9, b)"| = (2? sin 0), (A.2) 


provided />1, />m, and @ is not too close to 0 or z. Since the 
condition />m implies m/!~0, this approaches the classical 
case discussed above for m//=0. The approach to the correspon- 
dence (A.2) is illustrated in Fig. A3 where two examples of 
| Y, 2 with m=0 and with /=10 and /= 25, respectively, are 
compared to the classical expression (27? sin 0)~!. (See [98] for 
examples where m#0.) The value of / is not very large in the 
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‘i \ x = : rE ; L 
Ue Veh 
0.6 0.8 1.0 ile, 1.4 1.6 


6 (rad) 


Fic. A3. Illustration of the relation (A.2) when m=0. The smooth, solid curve 

represents the Satan (A.2) while the oscillatory curves are values of 

| Y, (8, )|? for /=10 (short dashes) and /=25 (long dashes). Note that 
|¥,(@— )|?=| ¥,o(8, o)|?, and both are independent of ¢. 


former case, and an averaging interval A@0>2//=0.3 rad =18° 
would be required to smooth the oscillations. None the less, the 
average behaviour of | Yio0(9, $) | 2 follows the relation (A.2). The 
larger value /= 25 would require a smaller averaging interval 
A@>0.13 rad7°, and the classical value (A.2) is approached 
more closely. None the less, the oscillations of | ed, ¢)|? about 
the classical value continue, no matter how large / becomes; the 
spherical harmonic passes through zero at intervals of x/l in 0. 
Thus the averaging over @ is always needed in principle before 
we retrieve the classical limit (A.2), even though the averaging 
interval required becomes smaller as / becomes larger. 
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Another spherical harmonic limit: One can relate spherical | 
harmonics to Bessel functions when /> 0 and 1>m230, with: 
l@>1 but @<1; then 


IN 1 


Y,_ (0, prem)  Snl(l+4)0. (4.3): 


When m=0, this reduces to 
Pxcos 6) x J,[(1+ 56]. (A.4) 


The asymptotic limits for the Bessel functions when /@>1 show 
that the relations (A.3) and (A.4) are equivalent to (A.1) when 6 
is small. 

These expressions can be used to show that equivalence of the 
Fraunhofer diffraction and partial-wave expansion descrip- 
tions of the scattering of a particle by a strongly absorbing 
(‘black’) disc or sphere [97, 99]. 

A generalization of the relation (A.3) holds for the reduced 
rotation matrix elements, 


di nO)=J,,_.ld+p0l, (A.5) 
again as 1300, 1>m, n, with m—n20, and /@>1 but 0<1. 


Vector addition (Clebsch—Gordan) coefficients: The interpreta- 
tion of vector addition (Clebsch—Gordan) and other angular 
momentum coefficients as probability amplitudes, and their 
visualization in terms of the vector model, was discussed in 
Chapters II and III. It suggests that there are expressions for 
these quantities in the large-j limit that correspond to the 
geometrical relations that hold for the classical analogues. 

Ponzano and Regge [100] (see also [96]) derived such 
expressions for the Clebsch—Gordan coefficient when all three 
angular momentum vectors are large, and showed that they 
have simple geometrical interpretations. 

Another class of asymptotic expressions for Clebsch—Gordan 
coefficients arises when a small vector b is added to large one a; 
consequently, the resultant c=a+hb is also large. Two such 
expressions were given on p. 33, taken from the detailed 
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discussion of Brussard and Tolhoek [15] (see also [96]). When 
a, c>b, 

<abaB|cy) ~ 5,445 Upe—a(8) (A.6) 


where 


cos 0=7/,/[e(e+ Dl] = y[(e+ DX y/E, 


so that 0 is the angle between the vector c and the z-axis. In the 
special case a=c and f=0, this reduces to 


<eby0|cy> = P,(cos 8). (A.7) 


These approximations are quite accurate even when a and c 
are not very large. To illustrate this, consider the relation (A.7). 
Table 3 (p. 36) may be used to obtain explicit expressions for the 
Clebsch—Gordan coefficient when 6b is small. For example, if 
b=1, the left side of (A.7) becomes 


<cly0|cy> =y//[e(e+ 1), 


which is exactly the right side for all values of c, 


PJ {ele + DI) =7/V [ele +1]. 


(An error of order 1/8c” would be introduced if the semiclassical 
approximation cos 0=y/(c +4) were used.) 
Again, if b=2, the Clebsch-Gordan coefficient is 


3y?—c(c+1) 


(02907) = Fe @e+3) en) 


while the corresponding Legendre polynomial is 


3y?—c(e +1) 


Piy/./le(e +1) = FEEH) 


The numerators are identical, while the denominators differ by 
relative amounts of order (2c+1)~”. 


n-j (n>8) coupling coefficients: Asymptotic expressions for 
coupling coefficients involving more than three angular 
momenta rapidly become more complicated and less easy to 
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visualize as the number of variables increases, although there: 
remain underlying geometrical structures that can be related to) 
classical mechanics in the large-j limits. Various asymptotic: 
relations for 6-j and 9-j symbols are included in the compilation \ 
of Varshalovich et al. [96]. Here we only present two examples | 
for 6-7) symbols (Racah coefficients). 

The first is known as Edmonds’ formula [22], 


abc oo (Se ; 
{ae © (sy Gace Fri O), (A.8) 


where m=d—b, n=e—a and a, b, c>1, f, m, n. The six vectors 
form a tetrahedron (Fig. A4), and @ is the angle between the 
vectors a and b, so 


_ aat+1)+ 66+ I—e(e +1) 


cos 0 
2,/{a(a +1)b(64+1)] 


(A.9) 


In the present case the vector f is short compared to the other 
vectors, and thus the differences m and n are small. In 
particular, if m=n=0, relation (A.8) reduces to 


abc 7 (yee rer Pp a oe 
loot “Jloanerrmy eee AO 


baf 


Fic. A4. Tetrahedron whose edges are formed by the angular momentum 


abc 
ae The length of the 


vectors indicated makes use of the semiclassical approximation Jla(a+1)} = 
a+, ete. 


vectors represented by the arguments of the 6-7 symbol } 
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The other example is a simple expression, known as Wigner’s 
formula [78], and applies when all six vectors which constitute 
the tetrahedron in Fig. A4 are large. If a, b, c, d, e, f>1, 


bc)? 1 
go ae (A.11) 
def 24nV 


where V is the volume of the tetrahedron [78, 96]. However, the 
relation (A.11) is only valid on the average because the 6-j 
symbols oscillate rapidly in value as the arguments are varied. 
An average of the left-hand side of (A.11) over a reasonable 
range of values of at least one of the arguments must be taken 
before it converges to the right-hand side. This presents another 
example of the nonuniform convergence of quantal results to 
the classical expectations. 

Extensions of the relations (A.8) and (A.11) are given by 
Varshalovich et al. [96]. 
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Addition theorems, 28, 55. 

Adjoint, 59, 60. 

Angular momentum, commutation 
relations, 10, 53. 

——, matrix elements of, 17, 58, 
126, 127. 

Asymptotic relations, 33, 154-60. 


Bipolar harmonics, 55. 


Clebsch-Gordan coefficients, see 
Vector addition coefficients. 

Central forces, matrix elements of, 
101-5. 

Classical limits, 27-9, 33, 154-60. 

Cogredient, 28. 

Contragredient, 28, 60, 62. 

Correspondence principle, 8. 


Density matrix, 107-12. 


Electromagnetic multipoles, 71-7, 
88-91, 98-100. 

— —, transition matrix elements, 
91-5. 

Equivalent operators, 99. 

Euler angles, 20. 


Fractional parentage, 83-6. 


Graphs, arrows on, 114, 119. 
—, basic components, 113, 114. 
—, normal form, 120. 

—, rules for constructing, 115. 
—, — — reducing, 126-33. 
Group symmetry, 3, 4. 

— representations, 13-15. 


Hyperfine interactions, 74, 76, 


98-100. 
Isotopic spin, 86. 


jm-coefficients, definition, 126. 
—, graphical representation, 126. 
—, reduction formulae, 127, 131. 


Matrix elements, electromagnetic 
multipoles, 91-5. 
—, factorization of, 57, 80-3. 


—, reduced, 57, 58, 61, 79-83, 85, 
151, 152, 153. 

—, rotation 21—4, 28, 32, 146, 147. 

— of spherical harmonics, 57, 81, 
82, 153. 

— time reversed, 62. 

Multipole expansions of a scalar 
field, 63, 64. 

—— — a spinor field, 70. 

— — — a vector field, 65-70. 

—-—w— a vector plane wave, 69. 

— —, special functions, 127, 131. 

Multipoles, see electromagnetic. 


Parity, 10, 19, 59. 
Pauli spin matrices, 23. 
Projection theorem, 78. 


Racah coefficient, 41-4, 140, 141. 

— —, explicit formulae, 43, 141. 

—w—, graphical representation, 
117. 

— —, sum rules, 44, 140. 

— —, symmetries of, 45, 140. 

Reduced matrix elements, defini- 
tion, 57, 78, 151. 

— — —, phase, 57, 151. 

— — —., reduction formulae, 79, 
81, 82, 152. 

— — —,, special cases, 57, 58, 153. 

Rotation of angular momentum 
eigenfunctions, 28. 

— conventions, 5, 21, 146. 

— of spherical harmonics, 28. 

— of spherical tensors, 51, 53. 

Rotation matrices, 21-4, 146, 147. 

— —, explicit forms, 22, 24. ~ 

——, reduction of product of, 32. 

— —, symmetries, 147. 

Rotation operator, 9, 20, 21, 53, 
63, 65. 


Selection rules, 59. 

Spherical harmonics, 18, 145, 146. 

— —, addition theorems, 28, 54, 
55, 146. 

— —, matrix elements of, 57, 81, 
82, 153. 
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Spherical harmonics, as tensors, 54. 
——, vector, 66, 67, 68, 149. 
Spin Hamiltonian, 100. 

Spin of a vector field, 65. 

Spin orbit forces, 83, 108. 
Statistical tensors, 109. 

Sum rules, construction, 44. 

—  —, Racah coefficient, 44, 140. 
— —, 6-7 symbol, 134, 142, 143. 
——, 9-7 symbol, 143. 

Symmetry and conservation laws, 


—, geometrical, 4. 
—, group, 3, 4. 
—, rotational, 5. 


Tensor Cartesian, 49, 50. 

—, irreducible spherical, 50. 

—, products of, 51, 52, 54, 55, 67, 
79, 148. 

—., transformation of, 51, 53. 

— forces, 106, 107, 149. 

—, operators, commutation rules 
with angular momentum, 53. 

—, adjoint of, 59, 60, 61. 

—, matrix elements of, 57. 

—, time reverse of, 61, 62, 91. 

Time reversal, 11, 61, 91. 


Vector addition coefficients, 30, 39. 
— explicit formulae, 33, 34, 35, 
36, 137, 138, 139. 
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— — —,, graphical representation, 
117. 

— recurrence relations, 34, 114. 

— symmetries, 136, 140. 

— model, for rotation matrices, . 
29. 

—w—, for vector addition coeffi- 
cients, 32. 


Wigner 3-7 symbol, 39, 126. 

— — —, graphical representation, 
113. 

—— —, orthogonality, 115, 116, 
136. 

— — —,, recurrence relations, 137. 

— — —,, special cases, 138, 139. 

———., see also vector addition 
coefficient. 

Wigner 6-7 symbol, definition, 4, 9, 
142. 


— — —,, graphical representation, 


118, 
— — —, sum rules, 134, 142, 143. 
— — —, symmetry, 142. 
———., see also Racah coeffi- 
cient. 


Wigner 9-7 symbol, 45—7, 143, 144. 
— — —,, graphical representation, 
125. 


X-coefficient, see Wigner 9-7 sym- 
bol. 
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Angular momentum introduces the quantum theory of angular 
momentum to students who are unfamiliar with it and develops it to 
a stage useful for research. 

The first part of the book contains the basic theory of rotations 
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